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This issue of AVTOMATIKA I TELEMEKHANIKA contained an article commemorating 
the 90th anniversary of the birth of V. 1. Lenin. Since that article did not.discuss either the 
present state of Soviet research in automation or its historic development, and as a matter of 
fact, contained no information of scientific interest, it was omitted from the translation, 














EDITORIAL BOARD OF 
AVTOMATIKA I 
TELEMEKHANIKA 
M. A. Aizerman 
L. A. Charikhov ‘ 
V. A. I'in | 
ya ea Automation and Remote Contro 
M. P. Kostenko 
S. P. Krasivskii 
ri aga A translation of Avtomatika i Telemekhanika, a publication of the 
(Aasoc. Editor) Academy of Sciences of the USSR 
B. N. Petrov 
Shumilovskii 
“ = en conn Russian Original Dated April, 1960 
B. S. Sotskov 
V. A. Trapeznikov 
(Editor in Chief) 
Ia. Z. Tsypkin 
G. M. Ulanov Vol, 21, No. 4 November, 1960 
(Corresp. Secretary) 
RUSS. 
PAGE PAGE 
Analytical Controller Design, A, M, Letov eeereeeere#ese’ . eoeenerneeeneeeeneeneeeneeeneee 303 436 
Application of the Krylov-Bogolyubov Method of Constructing Asymptotic Approximations to 
the Investigation of Systems with Lags, V. S, Kislyakov.......ececessesceecess 307 442 
On Steady-State Modes in Automatic Control Systems, B, V. Shirokorad. ....seeeesee08 316 456 
~*= Determining Pulse System Optimum Weighting Function Which Ensures the Extreme of a 
Functional, PF FP errrrrrrrrrTTTT TTT. 322 465 
~The Effect of Stationary Random Processes on Automatic Control Systems Containing 
Essentially Nonlinear Elements, G, I, Pyatmitskii..... cc ccecccccceseseseses 328 474 
An Interpolation Method for Analyzing Automatic Control System Accuracy under Random 
Stimuli, We Be CE, fo 6 6 6600 006606606006 060 cheese beens 333 481 
- Theory of Ideal Models of an Extremal Controller, V, V. Kazakevich .....ceeeeeee008 338 489 
Electromagnetic Control-Element Dynamics, L E, Dekabrun, N, P, Kozlov, and I, M, 
nn 66660066 000680000006000000006600000800066 60006046000" 348 506 
‘Differentiation of Slowly Varying Signals, S, N, Diligenskii..... ccc cecescecevceces 353 513 
Choice of Frequency Range for Industrial Pulse-Frequency Telemetry System Devices, 
i. i, PET 6 6000600000000 0206 66bseeeesneenesetesecaneen 361 525 
On the Reliability of Schemes for Connecting Dispersed Objects, V, A, Il"in .....eeeeeee 364 530 
On the Question of Relay-Contact Device Reliability, K, P. Kurdyukov......esceeessees 366 533 
Rectified Voltage Regulator, OO O Bee  ccctocceeseesesesecesseseseceuane 372 542 
Stabilizing the Temperature of Heated Thermistors, G, K, Nechaev,...eseeeeeeeesees 376 548 
Use of Heater Resistances for Current and Voltage Stabilization, V. S, Popov. ....ssseeees 381 555 
eatees tp the Béiecr. A. 1. CREUGNORSY .ccccccccoccceceseseeeeeeoceseeeceees 384 558 

























A. M. Letov 


(Moscow) 


Translated from Avtomatika i Telemekhanika, Vol. 21, No. 4, pp. 436-441,April, 1960 


Original article submitted December 23, 1959 


A solution is presented for the problem of analytical controller design in accordance with a given optimizing 
functional. The solution is given for the open region of definition of the system's differential equations. The cases 
of closed regions are worked out in Parts II and III of the present work, which are to be published in succeeding num- 


bers of this journal. 


1. Two Classes of Problems in Optimal 





System Theory 





The coastruction of optimal systems leads to the solu- 
tion of the allied mathematical problems, which are di- 
vided into two classes. The first of these is comprised of 
those problems which are related to the determination and 


calculation of the undisturbed modes of motion. Here, one 
has to do with those programs of automatic control for which ~~ 
the given undisturbed motion possesses the requisite extremal Problem 
properties. We shall call such systems optimal with respect 


ANALYTICAL CONTROLLER DESIGN. I 





to control mode. 

Systems for which the concept of optimality connotes 
the arrival at a given state in minimal time we shall find 
it convenient to call brachistochrone. 

In the other class of problems fall those coatrollers 
which guarantee the existence of given properties of the 
disturbed motions (transient processes). 

We shall call such systems optimal with respect to 
transient response. The set consisting of the object of con- 
trol, the optimal programming element,and the optimal 
controller comprises an optimal automatic coatrol system. 
In both cases the optimization problem can be treated as 
atwo~point boundary problem, for the solution of which 
one may apply all the methods of the variational calculus 
and, in particular, the latest methods, developed in [1-7]. 
However, certain differences do exist between these two 
Classes of problems. Thus, in the first case, the problem's 
solution is obtained in the form of a known function of 
time, which then becomes the basis for the design of the 
programming elements. In the second case, the analytical 

















form of the control law appears as some function of the ori- 


ginal coordinates of the control system, i.e., the problem 
consists of the construction of the controller's differential 


equation. Of course, the process of analytic controller con- 


struction is nothing else than the synthesis process. In par- 
ticular, the synthesis process is here given the same mean- 
ing as in[1-5]. However, we choose not to use this term 
since many other authors use it with differing interpreta - 
tions. On the other hand, the equivalent terminology "an- 
alytic construction” must suggest itself more strongly to 


engineers, since this terminology implies merely the search, 


by means of the tools of mathematical analysis, for the 
form of the controller's differential equation which will 
answer to the accepted criterion of optimality. In the pre- 
sent paper we consider the problem of the analytical con- 
struction (design) of optimal systems which are defined in 
open regions. In the succeeding papers we shall deal with 
the same problem as formulated for closed regions. 


2. Initial Posing of the Analytical Design 
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We consider a closed controlled system (cf. Fig. 1) 
in which the object's disturbed motion is given by the set 
of equations 


= %— (2 brat, t+ me) = 0 (2-4) 














(A =1,..., m) 

in the generalized coordinates » y, while the equation of 
the controller for the coordinate € remains unknown. For 
definiteness, we shall assume that all the by, and my are 
given constants on which no restraints have as yet been 
imposed. 

Let N be the open region in which (2.1) are given; 
any conditions in N of the form 


M1 (0) = qo, - ++ Ny, (9) = "no §(0) = Co (2.2) 


m (00) =... = 74, (00) =§(00) = 0 

we shall call natural boundary conditions of the problem. 
They mean simply this, that whatever the transient re- 
sponse which arises in N, it must terminate, for t = oo, 
with the system being found at the origin of coordinates. 

As the criterion for system optimality, we choose the 
integral - 

1(@) =\ Vae (2.3) 


0 
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of the positive definite quadratic form 


2.4 
V= >} ayn. -+ cé?. an 
k 


We shall seek such continuous functions €, 1; ...1, 
(class C,) with continuous first derivatives which give the 
integral in(2.3) a minimum value. The existence and u- 
niqueness of the solution of this problem were proven in [7]. 
We now give ourselves the task of determining the form 
of this solution as well as the possibilities of using it for 
the purpose of an analytic design. The integral K€) is 
a functional defined on class C, of functions, and its val- 
ue characterizes the integral of the squared error, weight- 
ed by the constants a, and c, which the system will have 
in the course of the transient response lasting until t * = 
@. The problem consists in writing in analytic form 


(2.5) 


F (&,&, m,-...7,) =0 


the control law which, in conjunction with (2.1), will 
form a stable system and will guarantee the existence of 








be the system's characteristic determinant. It is easily 
proven that, if 4,...,4, are simple roots of the equation 
4(u) =0, then the numbers -/,,..... rH, are also simple 
roots. We assume that jt,....fin are simple roots, Rez}, #0, 
and the roots are numbered in accordance with the in- 
equalities 
| Rep,<0O (k=1,...,2), (3.5) 
This latter step is always possible by virtue of the 
aforementioned property of determinant (3.4). The gen- 
eral solution of the equations of the problem will consist 
of a linear combination of exponential functions of the 
form cye BE and Coa ke” kt (k = 1,...,n), and will contain 
2n arbitrary constants 


Cy veees CoCry tee -rCgn - 
2. By virtue of (2.2)* at infinity we set 
Cuk 29 (kK=1,...,0) 


where the constants c,, can be chosen from the initial con- 
ditions. However, this operation is not necessary. 
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8. Solution of the Problem 





a minimum for the integral in (2.3). Let 
2 
m, myMe 
by = -, ae | Din, “Oe ’ ’ = 
m,m, m2 (3.4) 
A (pu) = Dine <ss4 ban — #2, 2 ae) “Oe 
et a 0, —dby—p,..-, —On 
Pp dees 2an, — bin yeeey —Onn— 


We are dealing with the Lagrange variational pro- 
blem, the procedures for the solution of which are very 
well known. Specifically, we first set up the function 


H=V+)) 8, (3.1) 
: k 
where the A, are arbitrary multipliers. 
Then, 
0H 0H 
= =, x 2am, — >) haben —_— 0, 
on, k a 
(3.2) 
0H 


» a 
Ee, = 2cé —2 Marg. (3.2) 


The equations of the variational problem have the 
form 


Re = — J dychj + Zax, 0 = 2c —D) may. (3g 
; j 


To these, (2.1) should be adjoined. 

The further steps in the procedure for analytic de- 
sign (excluding singular cases) are as follows: 

1. The equations of the variational problem must 
be solved jointly for the functions n,, and A,. 














3. In the remaining 2n formulas which define the 
functions 1 4,...,Mn, Aq,---»-An, We eliminate the n function 
of time cyePk®, as a result of which we find that 


A= Dp Aen, (i= te--- sms 


where the 4; ., are completely defined constants. 
4. We substitute the values found for the A, in the 
last equation of (3.3), thus finding the connatinlt equation 


& = Prt + Pata +--+ - + Pa%n: (3.0 
It remains to verify that (3.7), in conjunction with t 
initial equations (2.1), forms a stablet automatic contr 
system. This verification is easily carried out. Indeed, 
the optimal system found has a general solution which is 
a linear combination of the exponential functions oye? 
(k =1,...,n). Consequently, the numbers py (k = 1,...0) 
are the roots of the characteristic equation of system (2.1), 
(3.7). We obtained ideal controller (3.7) with rigid feed- 
back and with an infinitely fast servomotor [8]. Thus, We 
arrive at the conclusion that functionals like (2.3) lead @ 
controllers which can not be realized. 


(3.6) 





(3,4) 


the 


(3.6) 


1 the 





equation 
(3.7) 


with the 
control 
deed, 
rich is 
o,eF K 


»m (2.1), 
d feed 
hus, we 
lead t 





4. Possible Generalizations of Functional 


——_—_—_ 


(2.3) 


——_$— 


We now consider an indefinite quadratic form such as 


shat 2s Dear; + © (4.1) 





and repeat the formulation of the Lagrange problem for the 
functional 


1) = \ (T +V)dt, 


The value of the functional K€) characterizes the in- 
tegral of the squared position and velocity errors which the 
system will have in the course of the transient response 
lasting until t* = o. 

We remark that the solution of the problem may be 
simplified by eliminating the derivatives 1 ;, in accordance 
with the original equations. It is reasonable, therefore, to 











co 
1@)=\ V+ bar. 
0 
The computation will differ from the foregoing only 
in that now the expression 0H/3E = 2€ and does not equal 
zero. Therefore, instead of the last equation of (3.3), we 
get 


(4.2) 


2F = 28 — D} mr. (4.3) 
j 


This equation, imconjunction with (2,1) and the first 


n equations of (3.3), solves the given problem, The solu- 
tion procedure reduces to the implementation of the fol- 
lowing operations; 


1. The equations of the problem must be solved for 
the functions n;, A, and —. 


consider first the simpler functional of the form Let 
by —P,.--, bin ’ 0 ’ ’ 0, my 
bay ,*** 9 ban — [4 0 ,** 9 0, Mn 
A () = 2a, ,° 999 0 ...> by —p ’ —— bn, 0 (4.4) 
0, re) 2a ’ —— bin ’ ’ = ban —fp, 0 
0, pees 0, ’ — im, , ’ — Mn, 2c — ps? 


be the system's characteristic determinant, As in the case 
of (3.4), it possesses the property that A(y)= Ay -y). 
Therefore, the equation A4(j) = 0 will have the roots 


Hys--+eH,, (which we assume to be simple), with the proper- 


ty enunciated by (3.5), plus the roots -p4,....~Hn4;- The 
system's general solution will consist of a linear combina - 
tion of tees. functions of the form c,pe"k* and 
Crake MK (k = 1,. 
constants, 

2. By virtue of (2.2) at infinity we set Cp, 34k = 0 
(k=1,...,.n+1). The remaining constants can be deter- 
mined from the initial conditions. 

3. We now take the solution found for the function 
§. In the n+2 formulas found for the functions 1 4,...,0 p, 
— and £, we eliminate the functions of time c,eFit, 
Cee nett , as a result of which we find the equatiate 
equation 


é= 2 Pe — 7. 


Here, the p,, and r are completely defined constants, 
determined in the course of implementing the enumerated 
Operations. In conjunction with (2.1), (4.5) forms a stable 
Optimal system. The stability, as in the previous case, is 
4 consequence of the fact that the numbers /,,.. 


(4.5) 


oH na 


are roots of the optimal system's characteristic equation. 
In contradistinction to the solution of Section 3, we 
obtained a linear controller with a limited variable ser- 


, n+1) and will contain 2n+2 arbitrary 








vomotor speed. Such controllers can be realized if the 
region N is sufficiently small. The analytic design pro- 
cedure in the case when the motion of the system is in a 
closed region N, when the optimal soluticn can be found 
on the region's boundary, will be considered in the sub- 
sequent papers. 
5. Examples 
For a first-order system, the equations of the varia - 

tional problem will be 

n= byn+m—, h=—bvA+2ay, 2E = 208 -—- md, (5-1) 


To these there correspond the characteristic equation 
(5.2) 





pt — (b? + c) u* + am® + cb? = 0 
and the solution 
y= cye"* + ce", 
betty me Bab owt, 


Pa (Ha — 
pee) 


g-o— 





(5.3) 


b) etal ries) eget, 


The controller's equation has the form 
i i n 
Hi—b pa—b . 
m m = 0 


#1 (Hi — 5) ps (pa — be 
m m 














or, in expanded form, 





_ (er —>) (Hs — 9) 
m 


ba (5.4) 





4+ (pit pe — 5) E. 


For a second example, we consider the second-order 


system 


Th = burma + diay + mE, Ha = barns + dasa + mag.(5-5) 





bui—w, bis 
ba, bas — BP, 
A (pu) = 2a, 0, 
0, 242, 
0, 0, 





Let f14,/4g, and pg be the three roots for which Rey, <0. 


We then have 


3 3 
m= >) Artue*e, m= Dy Aa(u,d ee,» 


“ = (5.8) 
bd t ° : ust 
= >. As (p,) ‘8 Cy: = > pss (v,)€ : Cs 
s==1 s=1 


Here, 4,, 4,.and Ag are the minors of the first, sec- 
ond,and fifth elements of the first row of the determinant 
in (5.7). It is obvious that the equation sought for the op- 
timal controller will have the form 


A; (1) A: (142) A; (us) m1 

As (v1) As (#2) As (Hs) a 

Mis Adee «tee og |" O* 
His(yi) pwads(s) psAs(s) € 


The author sincerely thanks E. A. Barbashin and N. N. 


Krasovskii for their active participation in the discussions 
of the present problem, and for their helpful suggestions. 
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To obtain the equations of the variational problem, 
we must append 


hy = — dyyAy — bars + 205m, 
Ke = — dyads — Dagha + 2a2%, (5.6) 
QE = 2cE — myAy — mads. 


The characteristic equation has the form 





0, 0, my 
0, 0, ms 
—ba—p, —bn, 0 = 0. 
— bss — p, 0 
— mM, 2¢ — 2p? 
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APPLICATION OF THE KRYLOV - BOGOLYUBOV METHOD 
OF CONSTRUCTING ASYMPTOTIC APPROXIMATIONS 
TO THE INVESTIGATION OF SYSTEMS WITH LAGS* 


V. S. Kislyakov 


(Moscow) 


Translated from Avtomatika i Telemekhanika, Vol. 21, No. 4, pp. 442-455, April, 1960 


Original article submitted July 4, 1959 


The present paper treats of the application of the Krylov-Bogolyubov method of constructing asymptotic ap- 
proximations to the investigation of linear and nonlinear systems with lags. The effectiveness of this application 
of the method is illustrated by the example of the investigation of the Minorsky equation. 


Introduction 





As a second approximation to the solution of (1), we 


Today, as an effective instrument for the approximate take the expression 


investigation of nonlinear instruments, one may use the 
widely disseminated method developed by N. M. Krylov 
and N. N. Bogolyubov [1] for the construction of asympto- 
tic approximations. This method, in a simplified variant 
(when only first approximations are considered), is called 
the harmonic balance method [2]. 

The nub of the asymptotic method of Krylov and 
Bogolyubov as applied, for example, to a nonlinear dif- 
ferential equation of the form 


- . 1 

x -}+-w*x = ef (x, 2), me 
where € is a small positive parameter, reduces to the 
search for the general solution of this equation in the form 
of a formal series expansion in powers of the small para- 
meter €: 


t= Acosd-+ eu, (A, $) -+ e*ug(A, $)+.. 


in which ug{A, ¥), u,(A,¥),.... are periodic functions with 


» (2) 


period 27. An additional condition is imposed on the func- 


tions uy, Up,..., Which is expressed by the absence of the 


xz = Acosd + eu, (A, 9), (7) 


in which A and ¥ are defined by the equations 


A = eD,(A) + #*D;(A), 





: 8) 
§=0+6B,(4)+0B(4),  § 
and uy is defined by the formula 
mn 1 <1 Bn (A) cosny +h, (A) sin ny (9) 
mby Co n=2 ra 


where g,,(A) and h, (A) are coefficients of the Fourier 


series: tie 
g,(A) = = \i(A cos }, — Aw sin ¢) cos nd dd, 
1 (10) 
h,(A) = = \ {(Acos$,— Awsin $) sin nd dd, 
v 


and the coefficients Dy(A) and B,(A) in (8) are defined by 


first harmonics in the Fourier expansions of these functions. the corresponding formulas (cf. [2]). 


As a first approximation to the solution of (1), we take 


the expression 
z= Acosd, (3) 


in which the amplitude A and the phase angle ¥ are de- 
fined by the equations 


A=eD,(A), ¢=o+8eB,(4). ™ 


The coefficients D,(A) and B,(A) in (4) are defined as 
the first (fundamental) harmonics of the expansion of the 
nonlinear function in a Fourier series: 


D,(A) = — ah {(Acos$,— Ao sing) sin ¢ dg, 


v 


on (5) 
B,(A) = — — \i(A cos ?,— Awsin ¢) cos ddd. 
U 
(6) 


As a third approximation to the solution of (1) we take 


the expression 


z= Acos}-+eu,(A, $) + 6%(A, $) ete. (11) 


In automatic control systems (ACS) one is frequently 
most interested in the so-called stationary (steady-state) 
oscillations. This is related to the fact that all nonstationary 
oscillations, after a quite rapid transient response, approxi- 


*Some results of this paper were briefly presented at the 

Sixth Scientific-Engineering Conference of Junior Scien- 

tists on Automatic Control of the IAT AN SSSR (Jan. 19- 

21, 1959) and, synoptically, at a seminar on differential- 
difference equations at Moscow State University (May 28, 
1959). 
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mate to steady-state oscillations and, almost immediately In each controlled system, the argument o expresses the 
after the initiation of an oscillatory process, this process total (sum) pulse (signal) formed in accordance with the 
can be considered as being stationary. A stationary mode, control law adopted for the problem. It is assumed that 
as is well known, is characterized by a constant amplitude, the function f[9(t -T)] lies in a subclass A’ of class 
i.e., such that the derivative of the amplitude with respect A of functions [4], i.c., satisfies conditions of the form 


to time is zero. For the system described by (1), the steady {J{[e(t—-t)]}=0 for ¢ (t—t) = 0; 
state is defined by the condition df |o(t —*)] ] : 

: | dco (t — t) lo (t—v)==0 > k > 0, 

A =eD,(A)+e°D,(A)+...=0. (12) 


o(t—t)o[o(t—r)] > 0 


The method under consideration, i.e., the Krylov- for 9 o(t—t)#0, 
Bogolyubov method of constructing asymptotic approxi- ' 
mations, can be successfully employed, as will be shown where p[o(¢ —*)] = /[9(t—*)] —ko(t—*), andk 


below, for differential equations with lagging arguments. _— js a given constant. This function also satisfies the Dirich- 
With this, the Krylov-Bogolyubov method, when combined jet conditions in the control range 2| o|, and can be pre- 


with the Routh-Hurwitz criterion, turns out to be a simple _ sented in the form of a linear portion and a small nonliney 
and effective means both for finding the steady-state per- _ correction: 

iodic mode and for determining the boundary of the region f[o(t —+t)] =ko(t—1)+ ep[o(t—*)], (12) 
of stability of the controlled system with lags described by 
linear and quasi-linear differential equations with lagging where €(€ “ 1),introduced by means of the substitution 
arguments. The present paper is based on [3], in which x(t) = € x(t), is a small parameter. 

there is contained the foundation for the method of approx- 
imate computations of periodic solutions of differential 
equations with lagging arguments, as well as on the results 
of simulation studies. The foundation for the possibility im the form 
of using the asymptotic methods of Krylov and Boglyubov . (t) = a,(A)o(t) + a,\(A)s (t), (1.3) 
for the investigation of systems with lags will be given in 

the second part of the present paper. where a;(A) and a,(A) are the amplitude-dependent co- 
efficients of the linear element. 

It is required to determine the values of the coeffi 
Oscillatory ACS with Lags by the Krylov- cients a,(A) and a,(A) and to resolve the question of the 
physical interpretation of equivalent formula(1.3). For 
this we can use the standard electrical engineering con- 


We now find a linear element which is equivalent, in 
a first approximation, to the nonlinear element of (1.1). 
We write the differential equation of the linear element 


1. Equivalent Linearization of Nonlinear 











Bogolyubov Method 
































Fig. 1. Block cepts of active and reactive power. The active power de- 
schematic of an veloped by the executive organ during one period of oscil- 
LP “ ACS with lags. lation corresponding, in the given case, to average power, 
LP is the linear will equal the work of the force f{o(t-7 )] during a 
; portion and NE period T of oscillation, divided by the magnitude of T, 
Z NE is the nonlinear 1.€., 
element with , T P 
lags. + \ Me )ie(ae. (1.4) 
We consider an ACS consisting, in the general case, . ; ; 
of a linear portion and a nonlinear element with lags (Fig. The reactive power will equal 
1). The equation of the linear portion may, in general, { 4 ° T (1.5) 
be arbitrary? while the equation of the nonlinear element + \ fis(t—)]¢ (¢ a =) dt. 
has the form 
z(t) = f[o(t —*)}. (1.1) We now show that the value of the coefficient aA) 


can be obtained by equating the reactive powers of the 
nonlinear and the linear elements, and that the value of 
the coefficient ag(A) can be obtained by a corresponding 
equating of the active powers. Thus, in accordance with 
(1.5) and presenting f[o(t-r)] in the form of (1.2), we 


An equation of the form of (1.1), in which the lag 
enters into the argument of the nonlinear function, may 
be treated, in particular cases, as the equation describing 
the executive organ of an ACS. In this case, the function 


f[9(t-r)]° represents the acting generalized lagged obtain 
force engendered by the executive organ in accordance t The Krylov-Bogolyubov method permits the linearization 
with the value of the argument 9 (t-T), wherer is a of a system's nonlinear element independently of the other 


positive constant which characterizes the temporal delay. system parameters. 
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+\ flote— a(t —+)at = 


me 
a 
7 


e(t—t)o(t — 1) dt + 


€ 


+ —- els(t—z)16 (tL) a. (1.6) 


T 
T 

We assume, in a first approximation, that the quantity 
o(t), during the time interval T = 2%/w, varies harmoni- 
cally: _ 

a(t) = Acoswt, o(t) = — Awsinot. (1.7) 

By substituting the quantity in (1.7) into (1.6) for the 
reactive power, and introducing the new independent vari- 
able ¥ = Wt, we get 





T 
+\flet—)le (t —+)ae = a= cos wt + 
0 


on 
\el4 cosw(t—r)]cospdp, (1-8) 


v0 


eAw 
3 2n 





By substituting the quantity defined by (1.3) into (1.5) 
for the reactive power, and by taking (1.7) into account, 
we have 

1 i T 
F \ 4 (A)e(t) a(t —) dt - 
0 


1¢ T A) At%w (1.9) 
7 \ aa(A)o (a(t — 7) atm BE 


From a comparison of (1.8) and (1.9) we obtain the 
value of coefficient a;( A): 


2n 


- “3 \ [A cos w(t — t)] cos ddd. 
0 
In exactly the same fashion, by substituting f[9(t- T)}, 
in the form of (1.2), in (1.4) for the active power, we get 
T 
#\ fle (t— 2) 3 (tat = 


(1.11) 


g T 
=) o(t—s)a(t)de + F \ ele (t— aa (t)ae. 


0 


By substituting the quantity (1.7) in (1.11), we obtain 


T 
F\ flo (t— 916 (t)de = — Aosinet _ 


1 ° . 
— \ ¢[A cos w(t — t)]sinddd. (1.12) 


0 





By substituting in (1.4) for the active power, the quan- 
tity defined by (1.3), and by taking (1.7) into account, we 
have 


T T 
gr \ 44(A) 0 (t) a(t) dt-+ Jr\ ap(A) 3° (t) dt = 
0 0 
: oA. (1.13) 


From a comparison of (1.12) and (1.13) we obtain the 
value of coefficient a(A): 


Qn 


a, (A) = a | ¢1A cos w(t — 2)] sin pag. 


’ (1.14) 


The method of determining the equivalent coeffi- 
cients a;(A) and ag(A) by (1.10) and (1.14) can be called 
the power principle, or the energetic balance principle 
for a nonlinear element with lags. It is easily shown (cf. 
the Appendix) that, when the harmonic balance principle 
is applied in the case when n = 1, one can obtain the val- 
ues themselves of the coefficients. In accordance with 
(A.8) and (A.9), (1.10) and (1.14) can be given in the form 


a,(A) = (Kk -|- sai) cos Wt — =“) sin wt, (1.15) 





a,(A) = — (i +- al) eet _ Sais wt, (1.16) 


A o Ae 





where g,(A) and hy(A) are the coefficients of the nonlinear 
function's Fourier series expansion for the case n = 1. 


2. Determination of the Steady-State Per- 





iodic Mode Parameters in Systems with 





Time Lags 


As the first example of the use of the Krylov-Bogolyu- 
bov method of asymptotic approximations, we consider a 
system described by an ordinary first-order linear differen- 
tial equation with constant coefficients and lagging argu- 
ment: 





z+ kx(t—x) =0. (2.1) 


Equation (2.1) was investigated by Schmidt [5] by 
means of a rigorous method for extracting the roots of 
the so-called transcendental function? . Schmidt showed 
that (2.1) possesses a periodic solution of the form ‘ 

r=, coswt + c,sinot, (2.2) 
where c; and c, are arbitrary constants, with the condi- 
tion that w is defined by means of the equation 


1 


o=k= = (+ + 2n) (2.8) 


- (n=0,1,2,....) 


tE. Schmidt considered the transcendental function 


1 (w) = iw — k exp (— iwt) = — k cos wt + i (@ + k sin wt). 
(2.1a) 


In his notation, instead of k,r, w in (2.1) and (2.1a), there 
standA,h,¥, 
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By the Krylov-Bogolyubov method, the coefficients 
a; and a, of the equivalent system, defined by (1.15) and 
(1.16), will equal for the example under consideration, 
(€ = 0); 
1=keoswt, a, = ———— 
(2.4) 


and the system's equivalent equation has the form 


2.5) 
(14 — see t+ keos wre = 0. é‘ 


A periodic solution of (2.5) occurs with the conditions 
that 
k sin wt 
1— ——__ = 0, 
@ 


kcoswt=0. (2.6) 


As is easily seen, (2.3) follows from this. Thus, the 
first approximation gives a result which coincides quanti- 
tatively with the exact solution of the ordinary linear dif- 
ferential equation with lagged argument which is being 
considered. 


As a second example of the use of the Krylov-Bogolyu: 


bov method, we consider the questions connected with the 
search for the steady-state periodic mode (autooscillation) 
which occurs under certain conditions in nonlinear oscil- 
latory systems with artifically introduced damping. The 
system is described by the so-called Minorsky** differen- 
tial equation 


y(t) + ry (t) + wby (t) + 2gy (t — +) = py (t —2). 
(2.7) 


where r, and q are positive constants and p is a small 
positive constant. By considering p to be sufficiently 
small, we set 
p = Ev 
and present (2.7) in the form of a system of two equations: 
y= 2. z= — 2rz — 2qz(t — t) — wey + evz? (t — t). 
(2.8) 
On the basis of the theorem on existence and unique- 
ness proven by Krasovskii [3], one can assert that, if 
the quasi-linear system (2.8) has a periodic solution for 
€=0 (the so-called generating solution) y(t,0) then, if 
the condition |€|<€, holds, where &€ is some positive 
constant, there exists a unique periodic solution y(t, €). 
We start initially with the periodic solution for € = 0 
(the generating solution), defined by the system of equa- 
tions 
y= 2. 2= — 2rz— 2g2(t —t)— iy. (2-9) 


By defining the coefficients a; and a, of the equivalent 
system by (1.15) and (1.16) (for the case € = 0), we obtain 


the following equality: 


2qz (t — t) = ayz(t) + ay2(t) = (2.10) 


= 2q cos wrz (t) — 4 z(t). 


By substituting the value of 2qz(t-rT ) in(2.9), we get 
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os. 2= —b2—dyy. (2.11) 
where , 
i 2r +- 2q cos wt wes wp 
~~ 4 — 27 Sin ot 2 2qgsinet * 


—— 1 


@ 10) 


The characteristic equation of (2.11) has the form 


—hk +1 
D(d) = ahhh: wala aed = 0 (2.12) 
or 
h® + bh + b2 = 0. (2.13) 


Equation (2.13) is an ordinary algebraic equation, and 
one may apply to it the familiar classical formulas of the 
theory of linear oscillations, in particular, the Routh-Hur- 
witz stability criterion. By this criterion, the system de- 
scribed by characteristic (2.13) will lie on the boundary 
of stability and will manifest undamped periodic oscil- 
lations in case 


R= A,-, = 6, =0, (2.14) 
where Anes is the(n-1)th Hurwitz determinant. 
Condition (2.14) occurs when 
R= 2r + 2qcoswt = 0. (2.15) 


The frequency of the periodic oscillations is defined 





by the formula 2 

® 
w? = oaks (2.16) 

os 2q sin wt 
@ 
By eliminating q from (2.15) and (2.16), we find that 
w? — w (2.17) 
tg or=>-— — . 


It is easily seen that the periodic solution of (2.13) 
occurs under the following conditions: 


tg wt < 0. (2.18) 


cos wt < 0 
or 


> + nn Cot < e+ ann (n =0, 1, 2,.. .)(2-19) 


In case of equality, wr = 7 +2Tn, the frequency 
w of the periodic oscillations either coincides with, or is 
a multiple of, the natural frequency u» of the system, and 
the phenomenon of resonance supervenes. 

The frequency of the periodic oscillations is deter- 
mined by a graphic solution of transcendental (2.17), and 
must satisfy (2.19). The amplitude A is an arbitrary con- 
stant, as is ordinarily the case when linear systems are be- 
ing considered. It should be mentioned that, in the exam- 
ple under consideration, an arbitrarily small delay T > 
may give rise to periodic oscillations in a system described 
by (2.9). In the case wkent = 0, no periodic oscillations 
exist in the system described by (2.9) for any combination 
of the coefficients since, by our assumptions, the coeffi- 





** An equation of the form of (2.7) was investigated by 
Minorsky [6,7], and then by Pinney [9], as related to the 
search for periodic oscillations. 

+ T Independently of Krasovskii, an analogous theorem was 
proven by Khalanai [10]. 
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cients 4, r, q are positive nonzero numbers. 
We now determine the periodic solution of (2.8) for 


«#0. Asa first approximation, we give the periodic solu- 


tion in the form 


y = Acoswi, (2.20) 


where the amplitude A and the frequency w are determin- 
ed from the solution of the equation obtained from (2.8) 
by use of the Krylov-Bogolyubov method of equivalent 
linearization. Determining the coefficients of the equiva- 
lent system from (1.15) and (1.16), we have the following 


equations: E ‘ 
2" (4 —. t) = yz (tl) -+- agz (t) = 


(2.21) 


9 SIN wt 


——~ z(t). 


3 eo q 
j pA? cos G)TZ (t) _ = pA 


By substituting the values of 2qz(t-r) and p z*(t-T) 
in (2.8), we get 


y = Z. . —— bz = bay. (2.22) 


where 3 
2r -+- cos wt (29 _ z 4?) 


by rs sin wt . 3 : 
i— are (24 =< pA?) 











— “9 
_ sin wt /, ryt 
i— - a _- = pA?) 
In this case, the system's characteristic equation will 
be 
~ F 
DQ)=|~"* T* Shae geen 
—by —bs—k 
or 
h® + bsh + bg = 0. (2.24) 


For fixed valves of the amplitude A, characteristic 
equation (2.24) is an ordinary algebraic equation, to which 
we apply the Routh-Hurwitz stability criterion. By this 
criterion, the system described by (2.24) will lie on the 
boundary of stability, and periodic oscillations will occur 
in it if 


R(A) = Ay_, = bs = 0. (2.25) 


which occurs for 


R(A) =2r + cost (2g — + pA?) = 0, (2.26) 


With this, the frequency of the periodic oscillations 
is determined by the formula 
2 
2 is 


- sin wt 3 ; 
i-a—5 (29-7 p*) 


By eliminating 2q-3 p A®/4 from (2.26) and (2.27), we 
find 





® 


(2.27) 








i oe 
witht ee (2.28) 


2ro 


The expression for the amplitude is determined from 
(2.26): 


= | 8 (r +- q cos wrt) }'/: (2.29) 
“Lo Speosar |” 


Equations (2.28) and (2.29) must satisfy (2.19) and, 
moreover, the condition 


q| COS wt | mT (for q> S pAt) a (2.30) 
which, if satisfied, makes (2.29) physically meaningful 
(A is real). 

The joint solution of (2.28) and (2.29) is most easily 
begun with a determination of the frequency w of the steady- 
state periodic oscillations, specifically, with a graphic solu- 
tion of transcendental equation (2.28). For this, we present 
(2.28) in the form of two terms: 

= iS. ae (2.81) 
the first term of which is the equilateral hyperbola €, = 
a /2rw with the— and w axes for asymptotes, and the 
second term is the straight line Ey = w/2r, with slope 
-1/2r. Using algebraic addition of the ordinates of these 
two curves for fixed values of u andr, we obtain the sum 
curve (Fig. 2), whose intersections with the tangent curve 
in the second, sixth, etc. quadrants define the frequencies 
of the steady-state oscillations of which, if the graph is to 
be believed, there is an infinite set. However, it is easily 
shown that not all, but only a strictly finite number n, of 
the frequencies obtained from the graphic solution of (2.28) 
satisfy (2.30) and, consequently, (2.29). Indeed, it is clear 
from Fig. 2 that, as the w, tend to infinity (i = 2,3,..,n,..), 
where we understand by w, the frequencies obtained from 
the intersections of the tangent curves with the curve in the 
sixth, tenth, etc. quadrants, tan wjT- oo. With this, cos 
wr — 0 and, for any frequency w,,;, there can occur a 
violation of (2.30), i.e., the condition of physical realiz- 
ability of (2.29) no longer holds. Thus, as a first approxi- 
mation to the solution of (2.7), we obtain the expression 

n 
y(t) = > A, COS Ol. 
s=1 

in which the amplitudes A, (s = 1,2,...,n)are determined from 
(2.29) and the frequencies corresponding to them, which satis- 
fy (2.28) and (2.29), are determined from the graphic solution 
of (2.28). 

As a second approximation to the solution of (2.7), neg- 
lecting quantities of the order of smallness of p*, we get the 
expression 


y(t) = >) A, cosa,t + p >) uy (Ag, ds). (2-83) 


s=1 =] 


(2.32) 


where u, is defined by a-formula analogous to (9): 


1 Bn (A +h, (A) sin n 
Rost + 8 (A) cosny y . 2.20 


ow ne—i 





+The proof that every solution y(t) can be given in the 
form of a sum of partial solutions is carried out in the sec- 
ond part of this paper. 
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Fig. 2. Graphic solution of transcendental equation (2.28). 


and g,,(A) and f(A) are defined by (A.8) and (A.9). In 


ome yer cos 3w (t — Tt) 


w? d (2.35) 
Then, the second approximation has the form 





1 
Uy = 35 A® 


- n 
AS cos 30, (t — 
y(t) = >) A, COs wgt + FD} 2 cos =e +) 


s=1 s=] 8 


etc. (2.36) 


3. Construction of the Boundaries of the 








Stability Region for the System of Equations 





(2.22) and the Investigation of Stability of 





the Steady-State Periodic Mode (Auto-Oscil- 





lation) *** 


In our example, the boundary of the stability region, 
as was already stated, is defined by the condition 


R(A)= . 
= 2r + cos wt (29 — 5 pA?) =0 (coswt <0). 


(3.1) 

By introducing the notation d(A) = q-3pa*/ 8 and kw) 
= 1/| coswt| , we construct the boundary of the stability 
region in the r, d(A) plane. This boundary is some curve 
whose equation may be written in the form 


d(A)=k(w)r. 


For the fixed values w2=20, r = 0.1 sec, the boun- 
dary of the stability region is shown in Fig. 3. It divides 
the r,d (A) plane into the region of stability, where (A) > 
> 0, and the region of instability, where R(A) < 
<0. Since dA) is a function of amplitude, each fixed val- 
ue of A has its own representative point. Therefore, the 
oscillations corresponding to boundary points of the system's 
patameter values, defined by (3.1), may be stable for cer- 


(3.2) 
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Fig. 3. Boundary of the stability region 
of system of equations (2.22) in ther, 
& A) plane. 


tain values of A and unstable for other amplitude values. 


The oscillations on the boundary of the region will be stable 


if 





(ee 


dA Aeition™ (3.3) 


and unstable if this expression is negative. 

In analogy with the criterion for determining safe and 
dangerous boundaries due to Bautin [11], we may assume 
that the boundary is dangerous “in the small” if the follow- 
ing inequality holds in its neighborhood: 


(en), <0. 


If this expression is positive, then the boundary is safe 
"in the small”. Dangerous and safe portions of the boundary 


(3.4) 


***In investigating the stability of the steady-state oscil- 
latory mode, we take the amplitude A, = A of the funda- 
mental harmonic into account, but ignore the amplitudes 
A; (i = 2,3,...) as being small high-frequency quantities. 
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are separated from each other by bifurcation points which 
are characterized by the condition 


(a) ae 
dA R(A)=0 a 





(3.5) 


Bautin's idea may easily be generalized to the case 
of stability “in the large”, i.e., for Aa. Mathematically, 
this may be determined as follows. The boundary of stabil- 
ity of the system is safe “in the large” if the following in- 
equality holds: 


@ (A) (3.6) 


dA he > 0. 
If the inverse inequality holds, then the boundary is 
dangerous “in the large”. 
Since, for nonlinear systems with lags R.A) depends 
not only on amplitude but also on the frequency of oscil- 
lation, then (3.3) can be rewritten in the form 





(aia . 
dA ais ‘ 
aR (A) 16 
- _(@R (A) dw 
=( OA Dieta: + (“ R(A)mod A > 9. 
For the example under consideration, we have 
(22a) ; 
: - dA JR(A)=0— . 
ee 2 — si ane ne cate 
at Ap cos wt — t sin wrt (2q i pA \5 ’ 
do 2wpA sin wt (3.8) 





dA 
2 | (20 — S pat) (sin wt + wt cos wot) — 20 | 





where (3.9) 
: wo? — w* 
sin wt = = ; cus we = — 
® (29 _ ey pt) 2q Tr pA? 
(3.10) 


By substituting the values of (3.9) and (3.10) in (3.8), 
we get dR A) . 
dA /R(A)=0 


; 3pA [w* (4r2r + 2r) 4- 2rw? + + (wo? — w2)?] (3.11) 


3 
2 (29 —> pA?) [w? (2rt + 1) + 02] 








The point A = 0 is a bifurcation point, since 


a 


The oscillation is stable in the small, since the quan- 
tity dRA)/dA remains positive as A 0. According to 
Bautin [11], the boundary is “safe” in the small. The point 
A~ qa is also a bifurcation point since 


(3.12) 





(3,13) 





(SP), 0 = 9. 


Thus the question of the stability or instability of the 
oscillations for A > oo remains an open one, and requires a 


more detailed consideration, i.e., investigation of higher ap- 
proximations and taking into account of oscillations of the 
form Aj cos wit (i = 2,3,...,n). This investigation does not 
lie within the scope of the present paper. 
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Fig. 4. Oscillogram of the steady-state periodic process 
in a nonlinear system with lags. 


In conclusion, we compare our results with those ob- 
tained by Minorsky [6-8]. This latter author considered 
the appearance of self-excited oscillations in dynamic 
systems which were described both by linear and by non- 
linear differential equations with lagged arguments of 
the form of (2.7), i.e., essentially linear and nonlinear 
problems. In considering linear problems, the author ar- 
rived at the conclusion that there exist an infinite num- 
ber of frequencies of self-excited oscillations. In practice, 
however, in his words, only the basic frequency of the self- 
excited oscillations is ordinarily observed, which he also 
considers in his investigation of nonlinear problems. It 
was shown in the present work thar in an autonomous dy- 
namic system described by (2.7) there may exist, not only 
the fundamental self-excitation frequency, but higher fre- 
quencies as well, although their number is strictly finite. 
The periodic solutions are obtained in the form of sums 
of series of small parameters, and an estimate of the de- 
gree of approximation can always be obtained from the 
remainder term. Use of the Krylov-Bogolyubov method, 
in conjunction with the Routh-Hurwitz criterion, allows 
one to consider, moreover, the boundary of the region of 
stability of (2.7) and to determine its “dangerous” and 
“safe” segments. It should be mentioned that the theoret- 
ical results obtained from the investigation of the Minorsky 
equation by means of the Krylov-Bogolyubov asymptotic 
method were confirmed by simulating the system under 
consideration on an electronic analog computer. Figure 
4 shows the oscillogram of the steady-state periodic pro- 
cess (autooscillation) arising on the boundary of stability 
of the system considered. 


The author wishes to express his deep gratitude to 
A. M. Letov and Ya. Z. Tsypkin for a number of valuable 
remarks which have been incorporated in the present work. 
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APPENDIX 





Equivalent Linearization of a Nonlinear Element with Lags by the Harmonic Balance 





Method 
Let the function of the nonlinear element to be linearized have the form 
F =q@[o(t —7*)], (A.1) 
while its argument o varies harmonically: 
o(t) = Acos wt = A cos. (A.2) 


It is assumed that, in the range of control 2| 0 |, the function F satisfies the Dirichlet conditions and may 
be expanded in trigonometric series without free term, of the form 
F = 9 [A cos w (t —*)] = g, cos + gz cos 2 +++++¢, cos np f-.s- (A.3) 
e tAysiny + hesin 2p +---+h, sinnp+---, 
whose coefficients are determined by the Euler process, in which both members of (A.3) are multiplied by cos 
n¥d¥ to obtain the Bn (n = 1,2,...), and by sin n¥d¥ to obtain the hy (n = 1,2,..), and then integrated from 0 to 
+2. Then, by taking into account the functions’ orthogonality properties, we get 





RR (A.4) 
En -1{ 9 [A cos w (t — t)] cos ny dy. 
“/ (A.5) 
on 
i, =| @ [A cos w(t — t)] sin ny dy. 
0 
We give sin n¥ and cos n¥ in the form 
By substituting the expressions thus obtained in (A.4) and (A.5), we get 
- , y (A.6) 
Bn = tener @ [A cos 4] cosnay — Sinner | 9 [A cos $] sin np dy, 
$ 0 (A.7) 
on on 
ig \ @ [A cos d] sin np dp + SR 7O* \ @ [A cos }] cos np dp 
or bi od 
0 0 
En = 8, 00S nwt — h, sin not. (A.8) 
h,, =h,, cos not + g, sin not. (A.9) 


where 8, and h_ are the ordinary coefficients of the Fourier series for the nonlinear function without lags. We now 
substitute (A.8) and (A.9) in(A.3) and group terms with identical Fourier series coefficients. By using the theorem 
concerning sines and cosines of a doubled angle, we obtain 
F = [A cos (t — t)] = g1 cos w (t — t) + g2 cos 2w (¢ — t) +----+ Ay sin w(t — t) + (A.10) 
+ hesin 20 (t — t)+--- 

In the sequel we shall use either the form in (A.3) or the form in(A.10). Let the nonlinear element be given 
in the form of (1.2) and the equation of the equivalent linear element be written in the form of (1.3). We use the 
harmonic balance principle to determine a,(A) and a,(A). By substituting (A.2) in (1.2) and taking into account only 
the fundamental (first) harmonic of the oscillation, we get 


f [A cos w (t — t)] = kA cos w (t — t) + eg; cos } + eh, sin = 
= Ak (cos sin wt + sin sin wt) + eg; cos} + chy sind, 
where g, and h, are the coefficients of the trigonometric series defined for the case n = 1 by (A.8) and(A.9). We 
now compare the oscillation of (A.11) with the oscillations obtained by substituting (A.2) in(1.3). From a com- 


parison of the coefficients of the sines and cosines we obtain two equations for determining the unknown coeffi- 
cients ay and ag: 


(A.11) 





a; (A) A = Ak cos wt + #g1, — a2 (A) Aw = Ak sin wt + eh. (A.12) 
From whence . ch 
a, (A) = (x +- “e) cos wt — |G —sinwt, (A.13) 
hy\ si h 
a,(A)=— (x Ff 3) — -_ 3 COS WT. (A, 14) 
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The present paper discusses certain basic concepts and methods of automatic control system (ACS) theory which 
are related to the existence of steady-state modes of autonomous systems with the coexistence of several independent 


frequencies. 


We limit ourselves to the consideration of the broad- 
est class of autonomous ACS, the properties of which are 
described by a system of ordinary differential equations, 
given in the Cauchy form: 


w= fx (Ly, La, s--y In) +9 MN) (1) 


where the right member lies in some previously given class 
of functions, defined in a fixed region G, and does not ex- 
plicitly contain the time parameter t. 


(k = 1,2, ee. 


Such systems were called dynamic by Birkhoff [1]. 

Lyapunov [2] called each (undisturbed) solution of (1) 
(following Routh's example) a steady-state solution with 
respect to the initiation of disturbancet. Subsequently, 
however, it turned out that such a general term was in- 
convenient, and it fell into desuetude. It developed that 
it was more convenient to apply the term "steady-state 


modes” to the more specialized motions of (1) which corre- 


spond to stable equilibrium positions (or points of rest). In 
the investigation of an ACS a basic role is played by the 
problem of determining the steady-state modes (for exam- 
ple, resting points, periodic, almost periodic, and the most 
general recurring [1] motions) which are stable with respect 
to other, transient, motions which approximate to these 
steady-state motions (control quality). All other motions 
(for example, generation or autooscillation of missile 
steering, generation of incorrect alignment of the autono- 
mous contour of a power drive or steering mechanism) turn 
out to be undesirable, and they are considered from the 
point of view of the necessity for finding the limitation to 
be imposed on the ACS which will prevent their appear- 
ance in the operating modes (from whence the interest in 
the tuning of the controller's transfer ratio, the determina- 
tion of the possible external noise conditions, etc.). 


The introduction of the concept of the steady-state 
mode was an obvious step for linear or for two-dimension- 
al nonlinear ACS with continuous right members of (1), 
but this was definitely not so for two-dimensional ACS 
with discontinuous right members, and even less so for 
multidi mensional nonlinear ACS. For these, other funda- 
mental concepts also seemed unclear, concepts such as 
establishment (of the steady-state mode), roughness (of 
adjustment), physical realizability, and others. 
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It was known even in A. Poincaré's time that the two- 
dimensional system (1) with continuous right members 
could have, in the phase plane (in the given case, in the 
velocity plane of A. Liénard), only two forms of steady-state 
modes: an equilibrium position with one state, or a period- 
ic mode with many states. 

If the right members of the two-dimensional system 
are discontinuous, then oscillatory motions of a new type, 
not converging to those mentioned above, might appear 
(for example, such motions appear in the systems treated 
in [3,4]). 

In multidimensional systems described by (1), even 
with continuous right members, guaranteeing the existence 
(fort from - @ to+ oo) of a unique and continuous de- 
pendence of the solution on the initial conditions, there 
might appear very complicated oscillatory modes, with 
the coexistence of several, and even infinitely many, in- 
commensurable frequencies. 

No mention was made of the existence of such modes 
even in such fundamental textbooks as [5-7]. Moreover, 
Teodorchik [8] asserted} that, in many cases, one of the 
frequencies must, in the course of time, dominate over the 
others, such that after the lapse of some time the (estab- 
lished) mode will be nonetheless periodic (i.e., in actual 
systems, almost periodic modes do not exist). 

Such circumstances led to an incorrect estimate of 
the value of studying the aforementioned modes in de- 
signing ACS. In particular, the problem of uncovering 
the very simplest of them — the almost periodic ones — 
in an ACS in the case n = 3 was completely undeveloped, 
although the properties of such modes in differential equa- 
tions were described by such authors as Poincaré, d’ Anjou, 
Hayashi, Nemytskii [9], and others. 

The etiology of the longstanding incorrect opinion 
as to the nonexistence of actual steady-state modes with 


Presented at the All-Moscow Seminar on Nonlinear ACS 
Methods of March 18, 1959, and at the First All-Union 
Session of Mechanics of January 29, 1960. 

T Lyapunov assumed that the f, were analytic functions. 
+ However, his proofs suffered from mathematical in- 
correctness. 
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many states of a more complicated nature than periodic 
ones is related to a very essential factor having to do with 
experimental observations, as well as to the widespread 
use of the harmonic balance method without account be- 
ing taken of its domain of permissible applicability. The 
fact of the matter is that, in making our observations, we 
cannot distinguish between, for example, periodic and al- 
most periodic steady-state modes, since the accuracy of 
the instruments and the intervals of observation time are 
quite limited. These individually noncharacterizeable 
steady-state modes present themselves to us phenomeno- 
logically as recurrent motions, i.e., those which were in- 
troduced into consideration as the most general motions 
by Birkhoff [1] in 1928. They are characterized by the 
property that, for each €> 0, there exists its own interval 
of observation T such that any trajectory of duration T of 
the recurrent steady-state mode will contain in its €-en- 


velope all the trajectories (as time varies from ~ a to 
+). Thus, an operator taking an observation during a 


time interval greatly exceeding T will always touch upon 
a “diffuse cycle” (a periodic mode with noise). 

However, from the heuristic point of view, it is not at 
all a matter of indifference whether we have to do with a 
periodic steady-state mode, distorted by random noise or, 
for example, with an almost periodic mode. In the first 
case, an analytic investigation of the ACSis possible; in 
the second case,such an investigation is frequently im- 
possible for the reasons previously cited. 

Not only Birkhoff but, in the 1930°s, Andronov **, 
spoke of the importance of investigating almost periodic 
steady-state modes. Finally, Aizerman [10] came out 
convincingly in defense of the physical reality of these 
steady-state modes and the necessity of studying them 
analytically for designing ACS. 


1. Definition and Internal Properties of 





Dynamic Systems 





We preface the more graphical explication and illus- 
trations by a somewhat abstract definition of the dynamic 
system given by (1) and by a description of certain of its 
properties, borrowed from topological dynamics [11]. This 
digression allows an exact mathematic formulation to be 
given, but does not impair the understanding of the follow - 
ing simpler material. 

Let X be a topological space of certain elementsz € X, 
which we shall call a phase space, and let the x be phases 
(points). Further,let T be a topological group of elements 
t€T (such that every pair t, tT is put into correspon- 
dence with a third element i, 7). We call T the time 
group, with moments t. In X we define a transformation F 
which puts each point = € X into correspondence with 
another point of X which we denote by zF € X. The topo- 
logical product XT is mirrored by the transformation F in 
X, i.e., each pair (x,t), where z € X and ¢ € 7, corresponds 
to a point (x,t)F. Descriptively (axiomatically), we define 
the topological transformation group, denoted briefly by 
the triplet (X,T,F), which satisfies the following axioms; 


1. Identity axiom. In the group T there exists 
an initial moment of time 0 € T (0 need not be a num- 
ber) such that the pair t and 0 corresponds to t , and the 
pair (x,0) is put into correspondence, by transformation F, 
with the point (x,0)F = x. 

2. Homomorphism axiom. If the pair 41. 
ts € 7 corresponds to tg €7,then{ (x, Ty)F,t,)F =(x,t)F. 

3. Continuity axiom. The transformation F 
is one-to-one and continuous ( i.e., topological). 

If Xis a metric space with metric p, and the moments 
of time ¢ €7 are real numbers then, in accordance with 
the continuity axiom, (x,t)F is a continuous function of the 
collection (x,t), i.e., foreach ¢ €'7 and a given pair of 
numbers €>0 and z¢€ X there exists another pair of 
numbers 6 = 6(x,t)> 0 and fT =1(x,t) such that, for 
any y@X and t,@T7 such that p(x,y)< 65 and|t-t,|< 
< T, the following inequality holds: p[(x,t)F, (y,tyF] < 
<€. 

By a trajectory (motion) we mean the geometric locus 
of the images of the fixed point z¢X under the transfor- 
mation F, corresponding to all the ¢¢7'. 

If, forall t€7' the transformation F takes the set 
I © X into itself, then I is called an invariant. For exam- 
ple, a singular point (a rest or equilibrium state), each tra- 
jectory, the set of all trajectories passing through all the 
limit points of a fixed almost periodic trajectory (in the 
case of a metric space X, cf, below), are invariant sets. 

It is easily seen that the triplet (1,T,F) satisfies the 
three axioms given above and, consequently, is itself a 
topological transformation group (with the induced topo- 
logy of X). 

A nonempty closed invariant set I is called minimal 
( lin if it does not contain an invariant set as a proper 
subset. In the example given above, the invariant set of 
almost periodic trajectories is minimal. 

If, on L nin’ the invariant J is defined with respect to 
the points ¢€7' then we say that I,,,,,, possesses J-internal 
stability. 

Each trajectory (or resting point) possesses a natural 
(trivial) invariant: the transformatioa F takes it into itself 
for all t€ 7. 

Turning now to the other side, we consider metric in- 
variants. For this we assume that X is a metric space with 
the metric p, and the moments of time ¢€7' are real 
numbers. The set I,,j;,, has S-internal stability (in the 
sense of Franklin [12)) if, for any pair €>0 and 2 € J min 
there exists 5 = 5( €,x) > 0 such that, for eachy €1 min 
lying in the 6 -neighborhood of x [i.e., forp (x,y) < 6], 
there will hold, for all moments of time ¢¢7' (i.e.,for 

-mM< t<+@ ‘ 
pl(z, t)F, (y, t)Fl<e. (2) 
If, moreover, I ip is bounded, i.e., there exists a phase 
y € X such that, for each x €1 pinep (7, y) < const < 
+oo,then, by means of the Vitali covering theorem, it is 
easily proven that in this case 6 will depend only on €, 


**This latter was communicated to the author by El'shin 
who consulted with Andronov during those years. 
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and I in Will represent a recurrent motion (in the sense of 
G. D. Birkhoff[1}), i-e., for each pair ¢>0 andz€Ipmin 
there will exist a t, = +, (e, xr) >0O such that, for any 
other pair ¢,, ¢@7' one finds a number + = t (t) € (¢; — t,, 
te |-t,). for which 
pl(z, t)F, (x, FI] <e (3) 


Then, a fortiori, I,)jp Will consist of L-internal stable 
trajectories (in the sense of Lagrange), i.e., there will exist 
a y€X such that, for each fixed *€1,5, and all? G7 


pl(z, t) F, y] <covst < + 00, (4) 


Each recurrent trajectory is P-internally stable, mean- 
ing that it possesses the local recurrence property (repeat- 
ability) in contradistinction to the integral property of (3) 
(in the sense of Poisson [1]). i.e., for any pair €, r> 0 and 
any 7 ¢ Imin there exist two moments of time, ty < -T 
and t >T for which 

ol(z, t)F, rz)<e. (5) 


The converse is not true. It does not follow from the 
P-internal stability given by (5) that the trajectory is re- 
current, as given by (3). 

Since a continuous function on a closed set is uniformly 
continuous, it follows from the continuity axiom and Vitali's 
covering theorem, that a closed trajectory (cycle) has §-in- 
ternal stability and, consequently, presents a recurrent 
motion. 

Markov [13] showed that each bounded S-internally 
stable 1.45, Consists of almost periodic trajectories, i.e., 
trajectories such that each of their components is a Bohr 
almost periodic function, i.e., is expanded in a trigono- 


metric series 
co 


> (a; sin @t + bi cos (;t), (6) 


i=1 


where the a; and b; are constants, and the numbers w; need 
not be integers (in contradistinction to Fourier series). 

The following property holds for an almost periodic 
trajectory (x,t)F. Forany e€> 0 there exists a T, =T,- 
“(€) > 0 such that, for each t,€7 in the semi- 
interval [tg,t9+ Te) there exists at least one number T 
(the so-called almost-period with respect to €) for which 


pl(z,t+)F, (2, t)F]<e (7) 
for allt €7’. 

In particular, if (7) holds for € = 0, then the trajectory 
is called periodic (a cycle) with period r. In this case, all 
the wu (i = 1,2,...) in (6) will be integers. 

We note that dynamic ACS (1) may not satisfy the 
three enumerated axioms. For example, let X be the or- 
dinary n-dimensional Euclidean space, T the set of real 
numbers, and x) a continuous decreasing function such 
that (x9) <0 forsome x > 0. Then, the trajectory de- 
fined by the equation x = (x) does not exist (i.e., can 


+00 
not be continued) for t < § dz / q(x) < 0, if this latter 
Xe 


integral exists. 
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It is known that each L-internally stable trajectory (1) 
can be continued. 

A trajectory of the type in (1) that can be continued 
and is uniformly continuous in t is S-internally stable since, 
in this case, foreach €> 0 there exists a 5 = 5(€) > 9 
such that(7) holds for 0< tr < 5 and anyt. 

Local continuity in t is a necessary but not a sufficient 
condition for S-internal stability. For example, the trajec- 
tory y=1/(1-x) of thesystem x= 1, ¥ = -1/(1-x) 
which passes through the point (0,1) for t = 0 does not pos- 
sess S-internal stability. 

We mention still one more important circumstance, 
A trajectory of type (1) which can not be continued is 
necessarily unbounded. The converse is not true. There 
exist trajectories which are simultaneously S-internally 
stable and L-internally unstable. Tricomi [14) first en- 
countered such motions in a rotor (armature) of a syn- 
chronous motor which lagged the stator’s rotating magne- 
tic field by the angle 6, the motion being described by 
the equation 


64 A6+ Bsind = C. (8) 
where A> 0,B>C> 0, The law of this motion is 
defined by the continuous periodic function 9 = $(t) with 

iod T: 
perio 0= 94 21, (9) 














Fig. 1 
Fig. 2 


Fig. 1. Slipping mode of an asynchronous motor'’s rotor. 


Fig. 2. Cylindrical phase plot of an asynchronous motor’s 
rotor slippage. 


On the phase (velocity) Liénard{g ,@ plane, this law 
takes the form of the unbounded curve shown in Fig. 1. 

This law of motion also has another geometric inter- 
pretation. On a phase cylinder (Fig. 2), it is given by a 
closed curve (a cycle) with elliptical form, and is bound- 
ed. 

The dual interpretation just given has a curious mean- 
ing. According to Nemytskii's classification [9], each 
bounded S-internally stable and (simultaneously) asympto- 
tically stable (external, in the sense of A. M. Lyapunov) 
trajectory represents a steady-state mode. Consequently, 
in the sense of Nemytskii, the mode given by cat T, de- 
pending on the choice of the coordinate system, is simul- 
taneously both steady-state and not steady-state. 


FFI is known that the mode of (9) is asymptotically stable 
in the sense of Lyapunov. 
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2. Definition of the Steady-State Mode the domain of attraction. All these trajectories will also be 


almost periodic, with the same frequencies as the first one. 
The trajectory of a periodic steady-state mode can be u- 





We assume a priori that, in topological space X, a 
limiting operation is defined, and that the group T is a 
well-ordered set of elements tt. nique in its domain of attraction, 

Let Imin be embedded in the open region U. The set We note in conclusion that Nemytskii [9], following 
of points x which lie in U but do not lie in Ipjn is called Franklin [12], insists that each steady-state mode must be 
the open envelope of I,,in, and is denoted by U(I,,,). S-internally stable. He starts from the condition that "if 

If there exists a Ul min such that, for each of its at some moment of time two states were neighboring, it 
points, there is an increasing sequence 1; ¢ T (i = must follow that in the future these states must remain 


Oe Se neighboring, since otherwise the mode under investigation 
either would have in the past, or will have in the future, 
fa Glas s cic Khe Kia (10) some tendency to development consisting of a separation 
such that of these two states from each other, which would render it 
(7. tn) F ->T min for n—> 4-00. (11) impossible to speak of a steady-state mode”. These words 


seem essentially unconvincing and unclear, so that the 


I is called orbit 
ee Laas, © Called cruttally seympestionty cette (in the development in time must be understood only from the 


sense of Poincaré) in the positive direction. In th : 
) a point of view of the limitation given in (2). 


U(Imin) is called the domain of attraction of Ijin. Since every motion in time is a development, then each 
Correspondingly, one speaks of a negative direction, temporal invariant may, with equal justification, be taken as 

if the sequence in (10) is a decreasing one. an indication of nondevelopment (in the meaning defined by 
If there exist simultaneously two sequences, one in- him), stationary nature, or establishment in time. The choice 

ceasing aud the other Cocsuasing, one mys Mat 1... & of some temporal invariant or another must be determined, 

simply orbitally asymptotically stable. not by subjective conceptions, but by the objective problem 
If, moreover, there exists a point y€I pig such that, hich supervenes in the concrete investigation of an ACS. 


in addition to ( 11), there — the relationship On this score, one could give many convincing examples 
(2. tn) F—>(y. tn) FP for n—» + oo. (122) fom practice. 


then one says that L nin is asymptotically stable in the sense 9. On Certain Methods of Inv estigating Al- 
of A. M. Lyapunov [in, respectively, the positive, negative, 
and in both directions where, in the latter case, one says 





most Periodic Steady-State Modes in an ACS 





that I. n 28 unconditionally asymptotically stable (in the 
sense of G. Dirichlet) }. 


If the space X is a metric one ,then the definitions just 


given assume theii ordinary meanings***. 
Each mode of motion which is an internally and or- 
bitally asymptotically stable ( externally) I,,;, is called a 


The simplest method of investigation is a modification 


of Minorsky’s approximate “stroboscopic” method [16], well- 
known in nonlinear mechanics, this method being a further 
development of the principle introduced into science in 1883 
by A. Gulden and Linstedt [17)TTT. The limits of appli- 
cability of the “stroboscopic” method were determined by 


steady-state mode. The corresponding process of establish- Urabe [18]. 


ing the steady state is defined by the limiting process of 
(11) in the sense of the topology of X. 

Each 1,,jn Of steady-state modes has a domain of at- 
traction. For example, such domains do not exist for any 


For greater clarity we demonstrate this method by the 


example of a four-dimensional ACS without, however, 
lessening the generality. 


We turn to the example of a “blind” winged missile, 





invariant set of motions defined by the equationx +x =0. circling (without lateral banking) at a constant altitude un- 
Nemytskii [9] showed that the presence of a domainof er the action of residual voltage U (of a hysteresis magne- 

attraction.in an almost periodic mode essentially simplifies tic or "zero drift” vacuum tube amplifier). 

its structure. Thus, for example, each asymptotically stable 

(in the sense of Lyapunov) steady-state mode of a three-di- 

mensional system is either an equilibrium position (with one ¥}The topological transformation group (X,T,F) is a semi- 

state) or a periodic mode or an almost periodic mode (with ordered set if T is well-ordered [11]. 

many states) having aot more than two independent fre- ***We note that, in the monograph of Zubov [15], orbital 

quencies. The images of these motions in phase space are, stability, as given by (11), is called stability in the sense 

respectively, a point, a closed curve,and an irrational to- of Lyapunov, which can lead to misunderstandings. For 

roidal winding. Whether the analogous suppositions hold example, for an isolated orbit of a heavenly body (a Sput- 

for n-dimension systems, for n>3, is not yet known, but is nik moving in accordance with Kepler's law in a conserva - 

probable. tive central force field 4, = —prf|r|r (where r is the vector 
The domain of attraction of an almost periodic steady- of the body and yp is the gravitational constant), the orbit 

state mode can not keep its trajectories isolated. An in- is unstable in the sense of Lyapunov and orbitally stable. 

finite set of trajectories passing through the limit points of ttt Their work also served as the starting point for the de- 

these almost periodic trajectories will also be embedded in velopment of the results of Krylov and Bogolyubov. 
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In the particular case, the missile equations can be put 
in the form 


B+ eB + cB + ¢95 = 0. (13) 
b=b+e8 
(where cy, Cg, Cs,and cy are positive numbers,and 8, ¥ and 
5 are the angles of lateralslipping, course, and rudder de- 


viation from the neutral position),while the autopilot equa- 
tions are 


= f(s), o=i(b+vp—U)—8. (14) 


where f(o) is the piecewise-smooth velocity characteristic 
of the steering mechanism [of(o) > Ofor o #0, f(0) = 0], 
i > 0 is the autopilot's rigid feedback coefficient and v > 0 
is the transfer ratio of the accelerometer for lateral linear 
acceleration. 

We assume a priori that oscillations are observed in the 


system, and that the conditions for applying the “stroboscopic” 6. 


method are satisfied $+. We may then set 
B= 6, -| r(f)sinnt + y(t)cosnt, o = z(t)sinnt, 

(15) 
where 8B, = cz U/(cgv + c2/i), n is angular frequency, and 
x(t), y(t), and z(t) are slowly varying functions of time (in 
relation to the period 21/n). 

On the basis of the previously given theorem of V. V. 
Nemytskii [9], the state equilibrium positions and the per- 
iodic modes of the {x,y,z} phare space of the stroboscopic 
projections of the hyperplanes of the {8 ,8 5} phase space 
of (13)-(14) completely define all the bounded steady-state 
modes (possessing S-internal stability) of (13)-(14) [owing 
to the degeneration of (13), since the second equation does 
not depend on the first]. A study of the three-dimensional 
(induced) system of equations in the { x,y,z} phase space 
can be carried out by the methods presented in [19-25] and 
others. However, there are as yet no effective methods for 
determining the domains of attraction (except for what was 
presented in (23]). 

It so happens that, under certain conditions, the { x,y,z} 
phase space contains a unique stable cycle. It is then easily 
seen that the cylindrical {;¢,8,5, p} phase space of (13)- 
(14) will contain a unique (bounded) steady-state almost pe- 
riodic mode with two coexisting independent frequencies n 
and c,8», so long as the latter are incommensurable.’ - 

This same result is attained by the (exact) method of 
the topological product, the essence of which is as follows. 

By well-known topological methods[ 23] one proves 
the existence of a stable oscillatory mode in {8 ,8 ,5} phase 
space. c 

The cylindrical {8 , 8,5,¥} phase space is considered 
as the topological product of the {8, 8, 6} hyperplane and 
the {¥} circle. By this we arrive at the result given above 

It should be mentioned that the topological product 
method has thus far been applied only to those ACS whose 
equations are either degenerate, or reduce in some fashion 
to those groups of which one, at least, is amenable to the 
application of one of the methods of investigation given in 
{19-25}. 
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DETERMINING PULSE SYSTEM OPTIMUM WEIGHTING FUNCTION 
WHICH ENSURES THE EXTREME OF A FUNCTIONAL 


N. I. Andreev 
(Kharkov) 


Translated from Avtomatika i Telemekhanika, Vol 21, No. 4 pp. 465-473, April, 1960 


Original article submitted October 17, 1959 


A wethod is described for finding the optimum pulse system which ensures an extreme of a sufficiently 
general criterion in the form of a functional I = ©.(1y,.-Iny 4 ), where ® is a known function, and Ij are 


quadratic functionals of the weighting function of the system. 


The methods for the determination of the optimum 
linear pulse systems minimizing the mean-square-error 
at its output are given in [1,2]. The least mean-square- 
error criterion seems a suitable choice in many cases when 
selecting a dynamic system. Inpractice, however, when 
evaluating and comparing various pulse systems one often 
encounters cases where the most suitable criterion is a more 
complex probability characteristic of the accuracy of the 
system operation. A characteristic of this type may be, for 
example, the probability P of the system error not exceed- 
ing by numerical value a certain given quantity c. Ina 
specific case the probability P may usually be expressed 
in terms of mathematical expectation and variance of the 


separate components of the error, the components having dif- 
ferent distributions. (When the components of the system error 


also depend on one another the moments due to their corre - 


lation should also be considered. For example, if the error of 
the system is composed of the sum of mathematical expecta- 


tion a and a random component normally distributed with 
variance o*, then 


P= 








1 c (a+ zx)? 
Tame) |— “pede =F C04, 0%). ay 


In the case when the mathematical expectation of the 
error is a, and the first random component Z, of the error 
is distributed normally with variance o”, and the second 
random component Z, has a rectangular distribution with- 
in the interval from -b to+b, and Z, and Z, are indepen- 
dent, then 


P= vite: | {ar[-"t |} - 


I (a*, o*, b*). 
In the often encountered case when the quantity c is 
small, (1) and (2) assume a simpler form: 


p= V2 exp [— 4], 


V x0 2c? 
b 
ve __ (¢+y)? (2a) 
P= Vine \ exp | — | dy. 


(2) 


(1a) 
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The quantities a’, o?, are quadratic [3,4,6] func- 


mange €. the system weighting function. In a more general 
case a*,o*, F* also depend parametrically on time. Conse- 
quently, even the optimum system which at any given mo- 
ment of time is extreme depends on time and is thus non- 
stationary. 


We shall assume henceforth that the criterion I for 
selecting the optimum system is a given function ¢ of n+! 
quadratic functionals I; (i = 1,..., n + 1) which in tum de- 
pend on the weighting function of the pulse system. The 
problem of determining the optimum linear pulse system 
according to the I criterion is formulated in the same way 
as in [3,5] when determining the optimum linear system 
of a continuous type. 

We assume that there is a method of finding an opti- 
mum weighting function kp{n,m,€]* of a pulse system 
with given time of the transient process and such that the 
function will ensure the extreme of the quadratic functior 
al 


Thy = 0,1; (k[n, m, 


s]) ++... +9n/,(k[n, m, ®}) 


+ Ingi (k(n, m, e}) (3) 


for arbitrary values of the parameters 6;. A method is 
sought to determine the optimum weighting function 
k*{n,m,€] with the same transient process time interval, 
and such that it will ensure the extreme of the functional 
(4) 
IT=I(k[n, m, e])=O(,..., In4). 

The time is here expressed in integral units [1,2], Th 
time in conventional units is obtained from the integral 
units as follows; t = nT + mT + €T at any moment of 
t = mT at the discrete values of time, The quantity T# 
the pulse sampling period, otherwise the control time in- 
terval of the pulse system, and € varies from 0 to 1. 
*The weighting function K{n,m, €] of a pulse system with 
variable parameters of the equivalent (reducible) contin- 
uous part expresses its response at the moment t = n +m 
to a unit impulse applied to the input at the moment 
t=n-0. 
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The problem is solved at an arbitrary moment of time 
regarded as a parameter. 

The difference between the formulations of this pro- 
blem and of [3,5] lies in the latter continuous case being 
the optimum system in the class of all the possible linear 
systems with the same given duration of the transient pro- 
cessiand in our problem the optimum system is selected 
from the class of pulse linear systems with given pulse 
sampling period as well as duration of transient process. 
All the pulse system weighting functions with a given 
pulse sampling period form a linear manifold (which also 
holds good for weighting functions of continuous linear 
systems)t. Therefore, all the derivations and proofs ob- 
tained in [3,5] for linear systems of continuous type could 
be repeated by replacing everywhere “linear system” (of 
continuous type) with “pulse system" (with given pulse 
sampling period), We just state the fundamental result 
without reproducing the derivations, 

The necessary conditions which must be fulfilled by 
the weighting function k*{n,m,¢] ensuring the extreme 
of the functional Iare: ~ 

a) The function k® should ensure the extreme of the 
functional Iyy, i.e., it is required that: 





k® = ky{n, m, &, 9),..., On) lo,0,,3 (5) 
b) The following should be valid: 
él (k,[n, m, €, 6,,...- 8] 
a0, —" (6) 


i Pein 
(i, j mi,...,M). 
or (a special case) the condition below must be fulfilled: 
00 [7, (hk), - Ing (OD 


al, =@ (i=1,. 





Tr oe 


which illustrates the fact that the function $ reaches its 
extreme within the domain of variability of the Ij func- 
tionals. This special case is usually of no practical im- 
portance, 

The condition a) enables one to reduce the complex 
task of determining the function k® at which the extreme 
of I occurs to a less complex one of finding kp at which 
the extreme of Ij; occurs, Determining the function ky= 
= ko(n,m,€,@4, ..., On] is the principal stage in the 

N 


E[n, ¢] = { >) x(n — mjk[n — m, m, ¢] — Hx[n, el} + 


m=0 
N 


(the effect of initial conditions being neglected). 


solution of our problem. Once the function k{n,m,e, 
64, . - - »On] has been obtained explicitly for arbi- 


trary values of the parameters 6j, it is easy to deter- 
mine the values of the parameters 6; = 6;, satisfying 
the condition b)t¢. Where only approximations to the 
function are obtained, the values of the parameters 

6; = Oj9 are also evaluated approximately, for example, 
graphically, As long as the number n of the parameters 
6j remains small, the determination of the values Oj 

is fairly straightforward and easy. This is why the de- 
termination of the weighting function kg[n,m,€,64, 

- ,On] ensuring the extreme of Iq should have pri- 
mary attention. Methods of determining kp with data 
given in other forms have been sufficiently developed([6, 
7 and others). The functionals I, 1j, Iyy obtained when 
investigating the pulse systems differ from the corres- 
ponding functionals I, Ij, Ij; obtained when examining 
the continuous linear systems (the integrals in the lat- 
ter case correspond to sums in the former one). 

In the case of a pulse system, the determination of 
ko can be reduced to the solution of linear algebraic 
equations [1,2]. As is known, in the continuous case 
ko is determined as a solution of a certain linear in- 
tegral equation. The actual form of the equations de- 
pends in both cases on the way in which the data are 
given. 

We shall consider one method of determining the op- 
timum weighting function of a pulse system from the 
criteria (la). Let on the input of the pulse system be im- 
pressed a useful signal X[n, € ] which can be represented 
as a sum of a known nonrandom continuous function 
x(n, €] and a random function X°{n, €] distributed nor- 
mally with xero mean and correlation function Rx{n-m, 
n-#], The interference Z[(n,¢€] also occurs at the input; 
it is a random function normally distributed with zero 
mean and correlation function R,{n-m, n-2}, In order 
to simplify the subsuquent exposition, it is assumed that 
interference and signal are not correlated, that is, 
M[XZ]= 0, The pulse system at the output is given by 
the function HX{n, €], where H is a linear operator. The 
time of the system transient process is equal to N. 


In this case the error is equal to 


(7) 


N 


+ {a X°[n — m| k(n — m, m, ©] — HX° {n, el} + > Z(n— m\k[n — m, m, ¢] 
m=0 


m=0 


In the expression (7) the first term represents the mathematical expectation a of the error, the second and 


the third term being random components of the error. 





TA set of M is called linear manifold when its elements satisfy the following rules; when k, and kg belong to M, 
then a,ky + a9k, also belongs to M, whatever the arbitrary numbers aj, a». 
tit was shown in [3, 5] that in order to determine @ = Ojg one can utilize another condition: 





ol | al 9 


instead of condition b). 


(j wi,..., 8). 


323 








We denote by I, the variance of the error: 


N N 
I, = o?= M[(£E —a)*] = > > R,[n— m,n—lk[n—m, m,e] k(n—1,l,e]+ 
I=0 m=0 , 
NON (8) 
4. > > R,(n —m,n—I]k([n—m, m, se] k(n —1, l, e] — 
l=0 m=0 
N 
—2 > HR, (n—l, n-+-s)k(n—l, l, e] + HHR,[n +e, n+]. 
1=0 
In the last but one terin, the operator H on the correlation function Ry operated with respect to its second 
argument. 
We denote by I, the square of the mathematical expectation of the error E: 
. (9) 
I, =a* = {M [E}}? = { } xz {[n—m)|k[n— m,m,e] — Hx{[n, el}. 
m=0 
The functional 1; in this case takes the form 
(10) 


Ty = 61, +1. 


Bearing in mind what was stated previously, one has to find a weighting function ky [n, m, €,64) which will 
ensure the extreme of the functional (10). The necessary condition satisfied by such a function ky is [6,8] the van- 
ishing of the first variation of the functional Ip 


dl, {k[n, m, e, 0] + Ax[n, m, e}} 


aa 11 
dA "sates - 





for any weighting function x of a pulse system with pulse sampling period T and time N of the transient process. 
We rewrite the expression for the quadratic functional I,,[k + A«].in the expanded form 


Thy (ke + Mx] = Ty [k] + 24D, + A*Dy. (12) 
where 
N N 
D, = > x{n —l, l, el > k[n — m, m, ¢, 9,)(9,R,[n — m, m —1] + 
i=0 m=0 
+ 6,2 [n — m, n—1) + 2[n —1) x[n — m]) —2z[n— 1] Hz[n, e] — 
— 0,HR,[n —1, n+ €}}, 
N N 
D,= > ys x(n —I, l, e]x[n — m, m, ¢} (6, R,[n — m, n — 1] + 
l=9 m=0 
+ 6,R,[n —l, n— m] + r[n — m] x[n —]}). 
It follows from (11) and (12) that 
N N 
x{n—l, l, €]} k{n — m, m, ¢, 9,)(6,R,[n — m, n—1} + 
a {2 ue (as) 


+ 0,R,(n— m, n—1) + 2[n—l1)x[n— m}) —9,HR,[n—/1, n+ e] — 
—2a(n—l] Hx{[n, e]} = 0. 
In order that (13) be valid for an arbitrary function « from the admissible class, it is necessary that the 
relation given below be valid: 


N 
> k(n — m, m, ¢, 9,) (0,2, [mn — m, n— 1] + 9,R,[n —m, n—1] + 
mand (14) 
+ a2[n—l|z[n— mj) — 9,HR,[n—1, n+ ¢]—2[n—l] Hz[n, e] = 0 
(}=0,1,...,N). 
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To prove it, let us assume that the relation (14) is 
not valid for certain! =, 0 51, * N; the parameters 
n, €, 64 are assumed constant, When / = i, the left-hand 
side of (14) is equal to some S # 0, Knowing that « is an 
arbitrary function, we choose « in the following way: 


x[n—1,,1,,e]=1, x[n—1,l,e]=0 (1h). 


But if x is chosen in this way, (13) is not tue, Con- 
sequently, our assumption has been proved to be wrong. 

The relation (14) is also a sufficient condition for 
the functionk to minimize I, when 6, = 0. It follows, 
namely, from (12) that, whenk =k» satisfies (14), then 
Ti [ky + *) = Ty [k,] + De,\where D, = 0. 

It follows that J); [ky + x] > 11; [ko] for an ar- 
bitrary « which is what we had to prove. 

The values 6 < 0 are, as a rule, of no practical im- 
portance, 

The solution kp[n,m, € ,64] of the simultaneous linear 
algebraic equations (14) must be substituted in (8) and 
(9) to obtain the relations a? = a*(0,), o? = o?(8,) ; 
subsequently, by using (1a), it is easy to obtain P as a 
function of 6;. Then the value 64 = 6g) corresponding to 
the maximum of P is found from dP(@3)/d6@,=0. The 
optimum weighting function of a pulse system is in our 
case: k* = ko[n,m,€,649]; the greatest value of proba- 
bility - which points to the optimum system — is Pp = 
= P(649). 

The conditions which ensure the extreme of the func- 
tional I actually occurring at k° are not considered in the 
present paper. For pulse systems these conditions can be 
formulated analogously to those in [5] for continuous sys- 
tems. 

When dealing with a complex problem the main dif- 
ficulty lies in solving the simultaneous equations (14) or 
similar equations. In a number of special cases of a sim- 
ilar character, the simultaneous equations can be solved 
exactly [1]. In the general case, the solving of such or 
similar equations can be accomplished by means of various 
approximate methods. The method of the steepest descent 
[4,6] appears to be particularly appropriate. We shall also 
apply it directly to determine the extreme of the quadratic 
functional lyj- (A description is given of the main features 
of this method when it is utilized to find the minimum of 
the functional ly. 

Let the selection of the first approximation of k{n—m, 
m, €] = ky{n—m, m, ¢], based on arbitrary considerations. 
We find next the direction of the gradient in the normed 
space of functions k{n—m, m, ¢], that is, the direction 
along which the functional Ij; decreases most steeply. The 
norm of a function k is defined by 


N 
= on WO Bite ae 
|k(n— m, m, ¢}|| = >} k? [mn — m, m, ¢]. (15) 
m=0 


The function x}{n—m, m, €], which corresponds to 
the direction of the gradient,should ensure the extreme 
of the quantity 


dl i, (ki[n— m, m, e] + Ax, [n—m, m, ¢}) 





qa a—y (46) 
under the condition 
N 
>) ln —m, m, ee) = 1. (17) 
m=0 


It follows that this function should also ensure the un- 
conditional extreme of the expression 


dl 5; (ky [nm — m, m, €]+4+-Ax, [n—m, m, e]) 





ss dA oe, (18) 
+ u > x2(n —m, m, 8}, 
m=0 
with » being the Lagrange multiplier whose value can be 
found from the constraint (17) (in our problem the multi- 
plier does not appear to be real and therefore is not found 
from the condition (17)). Bearing in mind that Ip is a 
quadratic functional one is able to write down 
In (A, [n —m,m, e] + Ax, [n —m,m, &}) = (19) 
= I [ky] + 244, + 4°A,, 
where Ay, A, depend only on K and Kg but do not depend 
on A, Ay being of the form 
N 
A, = >} B,([n — m, m, ¢] Lx, ([n — m, m, €}. 
m=0 

The expression By[n—m, m, €] here depends only on 
ky, but not on Ky; L is a linear transformation (on «,) 
operator, It is easy to see that in the case often encoun- 
tered in practice, when Lx, = xy, the necessary condition 
for the extreme of the functional (18) is, with an accuracy 
of up to multiplier y, of the form 


%,[n— m, m,e] = B,[n—m,m,e], (20) 


As the second approximation of the function k one 
takes 


ky = ky + Apa, 


where x, is given by (20) and Ag is obtained from the con- 
dition that when A= Ay then the function (A) = 

= J\;(k, + Ax,) reaches its minimum (its extreme), Then 
ky = By It follows from (9) that Ay= ~(Ay / Ag). 

The third and subsequent approximations are deter- 
mined by applying analogous considerations. The evalu- 
ation of successive approximations can easily be pro- 
grammed for a digital computer. 

The described method can also be applied to evaluate 
the extreme of the I functional. Now the determination 
of Ky {n - m, m, €) is somewhat more involved (in the 

ntl 57, dl; 


expression aE | = ) one has to 
ab |,» 1, dB \, 
i=0 
evaluate also the derivatives ir when k = k;) and the 
i 
procedure of finding the consecutive values 4, 4y,..., 
becomes appreciably more difficult because of the dif- 
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ficulty of determining the extreme of the generally com- 


plex function —(A) = (k; + A«;). However, in this 
case it is not indispensable to determine the values of the 
parameters 6; = Oj9 which were introduced only as an arti- 
fice. In each individual case one has to decide whether 
one should try to find directly the extreme of I or whether 


N N 


to find first the weighting function which will ensure the 
extreme of Ij; at various values of the parameters 6; and 
subsequently determine the values 0; = Gig- 

In the specific case when the functional I); is given 


by (8), (9) and (10), Ay and A, are expressed by means of 
ky and x, as follows: 


A; = > >) k(n — m, m, ¢}(0,R,[n —m, n—1} + 0,2, [n — m,n —tj + 


l=0 m=0 


+ 2[n+l)x[n—m)) —0,HR,[n —1, n +e] —2[n—l] Han, el} x [n—l, 1, €], 


N N 


A, = )) > (6,4. [2 — m,n —1]-+0,Re[n — m, n—l]+2 [n—m] z[n—l]} x 


l=o9 m=(0 


Xx, [n —l, 1, e]x,[n — m, m, s], 


N 
% = B,= Sa [n—m, m, ¢] (0,, R,[n — m, m, ¢] + 6,2, [n — m, m, 6] + 


m=0 


+ z[n—l]x[n — m]) —0,HR,[n —1, n, &] —z[n—1] Hxf[n, ¢], 


N N 
° YY} 4 la—1L 1, e] >! B3 [n—1, l, e] 
l=0 


l=0 
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As a rule, the method of steepest descent implies a 
rapid convergence. In each case the speed of convergence 
depends on the starting element k, as well as on the form 
of the functional taken as the norm of the function k . The 
starting element k, should be taken such that its norm be 
finite. 

Example. Let useful signal X and interference Z 
be impressed on the input of a linear pulse system. The 
useful signal may be represented in the form X = vt + X°, 
where X° is a random variable with a large variance, v a 
given quantity, and t the time. 

The interference behaves as a stationary random pro- 
cess of normal distribution with the following characteris - 


tics: its mathematical expectation M[Z] = 0, its correlation 


function HK, (t) = e—*'*!, The variables X and Z are inde- 
pendent. 

It is required to determine the weighting function of 
the pulse system with duration N = 1 of the transient pro- 
cess (the pulse sampling period is T) which would ensure 
maximum probability that the extrapolation error 

E= X[n+e+eJj—(X [rn] + Z[n]) «[a, 0, e] — (21) 
— (X [(n—1] + ZIn—1)) k [xn — 1,1, €] 
shall not exceed by numerical value a small quantity c. 
The probability of | E| sc, with c small, is found from 
the formula (1a). 

The weighting function k selected is subjected to 

the restriction 


k[n, 0, e] +k [n—1, 1, e] —1, (22) 


which is feasible, since variable X° must have a large 
variance, It is known that the condition (22) produces 
nonstability of the xero order. Taking into account the 
initial data in our example and (22), we shall pick k in 
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A; 


the class of stationary weighting functions, Consequently, 
k = k{m, €] and (22) becomes 


k{0, e] +k (1, e] = 1. (22a) 
Expression (21) can be rewritten in this case as 
E = X°+0T (n+ + eg) — X° —vTnk (0, e] — 
— wT (n—1)k 1, €] — Z[n] & [0, e] — 
—Z[n— 1] k[A, e] = oT (e + ¢,) + Tk (A, e] — 
— Zn} k 0, e] —Z[n— 1] k (1, e], 


where €eT is the extrapolation time, 


We determine a®, o, and Iq: 


a? = (M [tj)? = oP THe +e, +4 k[A, €])*, (23) 


ot = M [(E — a)*] = k* (0, e] + 


My a (24) 
+2 T R[O, e]k[1,e] +f. el], 
Ty; = 9,0? + a* = 6, (k* [0, ©] + 

+2 T R[O, e]R[4, ©] +47 (1, €]) + (25) 


+ v*T?* (e+ e,+k[1, €])*. 
It is necessary to find function ky minimizing Ij; under 
the coridition (22a). This function should ensure the uncon- 
ditional minimum of the functional 


=I, + 27, (k (0, ©] + #1, €}), (26) 


where is a Lagrange multiplier which can be determined 
from (22a). To determine ky we utilize (14) which in our 
case can be written as 





5) 


inder 
1con- 


6) 
nined 





@ 
0, [0, €] + Oe r k[i, e|] =— %, 











, (27) 
6,¢ 7 k[0, e] + (0, +27) k[1, e] = — To — #*T? (€ +-€,). 
Solving the simultaneous equations (27) and (22a) we obtain 
Oe 7 [yo + v2 T2(e + 8)] — to (01 + oT?) 
k,|0, ¢] = . a ’ 
0, (6, —O,e 7 +23 r) 
a 
1 0, (¢ od _j-erece 
bef, 0] = — : (28) 
—2 —_— 
0; (6, — 6,e r -{- v? 1) 
all) * r) 
6, (9; (, —~@ ) + v®T?] +- v®T? (e+- aC —O,- * 
To = bt . - 
20, (, J _ ier, 
Formulas (28) enable one to determine the function ky 3. N. 1. Andreev, "The determination of optimum dynam- 


when the values of € and@, are given. By keeping constant 
the value of the standardized time €, one should find ky for 
several values of 6; = 0 and determine subsequently the func- 


tion P = P (&); it is easy then to evaluate the value 6 = O49 4. 


which corresponds to the maximum of P. 


CONCLUSION 
The presented method of finding the optimum pulse sys- 
tem ensuring the extreme of I seems to be more laborious 6. 
when compared with the method of finding the pulse system 
ensuring the minimum of the mean-square error of the sys- 1. 


tem but basically both methods are equally complex. The 
criterion I is sufficiently general; it can be utilized when 


selecting automatic control systems designed for various pur- 8. 


poses. This method can also be used in the selection of op- 


timum nonlinear systems [9]. 9. 
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A method of calculating nonlinear automatic control systems under the action of stationary random pro- 
cesses is proposed. The method presented is a generalization of the method of statistical linearization and 
allows one to take into account the distortion of the spectrum of a random process by nonlinear elements. 


At the present time the statistical calculation of auto- 
matic control systems containing essentially nonlinear ele- 
ments is accomplished by the method of statistical linear- 
ization [1]. In essence this method consists in replacing a 
nonlinear element by an amplifier equivalent to it in some 
sense. The method is fairly simple and therefore is widely 
applicable in practice. 

However this method may never be considered com- 
pletely justified since the replacement of a nonlinear ele- 
ment by a linear amplifier does not take into consideration 
the distortion of the signal spectrum by the nonlinear ele- 
ment, and this distortion may in some cases prove to have 
a significant influence on the action of the system. 

In this paper a more systematic approach to the pro- 
blem of the statistical calculation of nonlinear systems is 
presented, in which the method of statistical linearization 
appears as the first approximation to the solution of the 
problem. As a consequence,one may estimate the appli- 
cability of the method of statistical linearization to a de- 
finite class of automatic control systems. The method 
presented may thus be considered a generalization and de- 
velopment of the ordinary method of statistical lineariz4- 
tion. 


1. Formulation of the Problem 





We consider a system with feedback shown schemati- 
cally in the figure. It contains one nonlinear element 1 
with so-called standard characteristics [2]. All the remain- 


ing links are assumed to be linear with constant parameters. 


Any system containing one nonlinear element may be 
brought into such a form. 

Suppose link 2 has weight function g(t) and frequency 
characteristic G{w), and the link consisting of 2 and 3 in 
series has weight function k(t) and frequency characteristic 
K(w). 








I zt; vl: J. 
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We suppose the input to the system to be a Gaussian 
stationary random process x(t); let m, be the mathemati- 
cal expectation, and o2 and RA ) the variance and corre- 
lation function of the input process. 

The problem consists in calculating the variance of 
and mathematical expectation my of the system's output 
process y(t). 

The process 2 t) represents the difference between the 
input process x(t) and the process w(t), which appears asa 
result of z(t) acting through the nonlinear element and the 
succeeding linear elements. The presence of the nonlinear 
element leads to the fact that the distribution law for the 
random process at its output will necessarily not be Gaus- 
sian. However, the process is “normalized” by passing 
through the linear elements, so that the distribution law 
for the process w (t) will be nearer to normal than the 
distribution law for the process at the output of the non- 
linear element. Considering that the process z (t) consists 
of the normal process x (t) and the “normalized” process 
w (t), one may assume that the distribution law z (t) will 
be near normal. Proceeding from this, we make the 
assumption that the process z (t) possesses a normal two- 
dimensional! distribution law. 

This assumption provides the possibility of solving the 
problem posed comparatively simply, since it is well know 
how normal.fluctuations are transformed by a nonlinear elt 
ment [2]. A similar assumption is made also in the method 
of statistical linearization. 


2. Formation of the Basic Equation 





Let R(T) and Ry (tT) be correlation functions for piv 
cesses z(t) and w(t) respectively, and R,,(T) and R,,{T), 
the mutual correlation functions for these processes. 

We have the following equation: 


Rx (t) = Rz(t) + Rw (t) + Rew (t) + Rwz (*)- (1) 


We shall find expressions for R(T), R(T) and Ryd 
in terms of RAT). 
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The normal two-dimensional distribution law for the 
process z(t) at m, = 0 may be put into the form of a series 
[3,4] 

- (2) 
) a bh ) one (= o"—» a 
f(z, Z2) = rv) y —h 
z 


z 


n=0 


where pz(rT)=Rz(r)/oz, 
gral defined by 


(x) is the probability inte- 


1 - me . 
(2) = = |e 2 dt, (2°) 


and #) (x) is the nth derivative of this probability integral 
with respect to x. 


Let the characteristic of the nonlinear element, that is, 


the connection between the output v(t) and input z(t) of the 


nonlinear element be given by the relation 


v = F(z). (3) 


Characteristics for a different type of nonlinear element 


are introduced in [1]. 

Using (2) and (3), one may calculate the mathematical 
expectation and correlation function of the process v(t) at 
the output of the nonlinear element [2]: 





my = \ F (z + m,)d® (+), (4) 
Rite Sint . oe *) (5) 


n=1 


where 


Yn = \ F (z + m,)d@™ (=-). (6) 


oo 


The correlation function R(t) may be written in the 
form 


R,(t) = 


FS 


-\{ k (a) k (8) Ry (« — @ + 8) da dB 


or, considering (5), 


mee 


Ry (t) =) = wt \ (a) (B) 92 (* —a + B) da dB. (7) 


n=1 00 


3. Solution of the Integral Equation 





Analogously, for Rz,{T) and R(T) we have the 
expressions 


Rrw (*) = — 2m | k (a) p,(t + a) da (8) 
0 


and 


Rus (t) = — 92m | k(a)e,(¢—2) da. (9) 
7a 


Substituting (7), (8), and (9) into (1), we obtain an 
integral equation for p,(rT): 


o2p (t) — 7,0 JRO, (t + a) 4-p,(t —a)] da + 


+ mp\ | &(a)k (6), (« — a + B) dadp + 
2 ee 
+ Dor) kk Bore —a + B) da dg = otp, (2). 


(10) 


The equation obtained is nonlinear, since it contains 
a power of the functionp,(r) higher than first. However, 
the nonlifiear terms contain the factor n* fnt which, as 
an immediate calculation shows, is small in comparison 
to unity and is less that ni. This circumstance allows us 
to apply the known method of expansion in powers of small 
parameters to the solution of (10), Such parameters in this 
case are n*, /n!. But it is not convenient to expand in 
powers of these parameters since the number of them is 
infinite. Hence, it is more suitable by far to peocead in 
the following manner. We introduce the factor y "in 
front of the terms of the sums in the left portion of (10) 
and expand in powers of 1. From here on the solution 
proceeds according to the usual rules. Finally, one sets yu 
equal to one in the resulting expression. This may be 
done, because it is not necessary to consider p to be a 
small quantity (in comparison with unity), since nt /n! 
are the small parameters. 


Taking the above into consideration, we rewrite (10) in a somewhat modified form, We transfer the sum with 
the nonlinear terms to the right side and adjoin the factor  " before the terms of this sum: 


oo 


ot, () —9,n, \k(a) [p, (¢ + 2) +p, (t —a)I da + 


aa i 


ou" 48 


\k (a) k (B) p, (« — 2 + 8) dadp = 


(11) 


= o8p, (t) — Syun | Qe carecere (« —a +B) dads, 


n=2 
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We represent the solution to (11) in the form of a series: 


co 
P2(t) = > wp, (2). (12) 
m=0 
Substituting (12) into (11), we obtain 


Sum fot, (2) — 142 NiO log ( $2) + Pq (t — a} doe + 


m=(0 0 


e. © 


ng \ }&(a)k(B)p,,(* —a + 8) dads} = a8 9, (2) — (13) 


co 2 co co 
_¥ 5 urtmetmat ... +My n— | \ Prn,Pm, Pim * * + Pon, (a) k (B) da d8. 
00 


n=2 ™,, ..., My=0 


Now it is necessary to equate terms with equal powers of y on the right and left sides. As a result we ob- 
tain a system of equations allowing us to determine successively the functions p af? ), expressing each in terms 


of the preceding. For concise notation we shall consider the left side of (11) as some linear integral operator 
L acting on the function p, (rT): 


Lp, (t) = 2p (t) — 2,7, | k (a) [p(t +a) + p,(t— a)|da + 


(14) 
+ mp | &(@)k() on (© — a + B) daa. 
Considering this, we obtain the following system of equations: 

Lp, (*) = 029, (t), (15) 
Lo, (t) = 0, (16) 

2 co CO 
Lpa(e) = —F-\ Se (a) k(B) 98 (« — a + B) dap, an 

, o> 

Loa (*) = — ar \ | (a) k (B) ef (& — a + B) dads, (18) 


Lp, (2) = — 34 k (a) k(B) pf(« — 2 + B)daedp— 


798 


| 2pspok (a) k(B)dadp, (29) 


Los (t) = — —2(( e(ayki@yetdadp — 
— 2B ( k(ayk(p) 3980, ded p— TEU | ka) @) 2ee pada. (20) 


t 


This system of equations may be easily solved by Fourier analysis. If we introduce the notation 
¢ 21 
Sm (w) = \ Pn (t) e-torde, ( ) 
—oo 
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then, taking the Fourier transform of (15)-(20), and performing simple transformations, we obtain the following 


expressions for S_,, (w): 




















‘ob 4 2 
500) = Te aR or = (> (22) 
S,(o) =0 (23) 
K (iw) |? f (24) 
5,(w) = — = Laat 3 . = = \ Sq (v) Sy (wm — v) dy, 
Kiio)\? 2 : 
S;(o) = 5, rh, =r ey a So (vz) Sq (v2) S, (@ — vy — ve) dvydve, — 
K(io)I* 7% 4 ¢ 
S,(o) = — 7 rua ++ \ 2S, (v) S, (@ — v) dv + 
—co (26) 
" 
TS a oa \ So (v1) Sq (v2) So (vs) So ( — vy — v2 — ¥%) dvydygd |, 
K(io)|? [7% 
54(0) = — 75 | ee { 25, (v) Ss (@ — v) dv + 
7 so 
+ Hr caaye J) 350 (4) Sy (v2) Sa (@ — 4 — ve) dradvy + (27) 





+3 J Solos) So (02) So (on) 5 (04) So ( — 1 — 9 — v0) dndraed | 


: 2S Oe Beh or ee eles “«- s 


Here K(iw) is the frequency characteristic of the 
serially joined links 2 and 3 (see the figure), correspond- 
ing to the weight function Kt). 

Considering that 

co 
i 
Pm (t) = =— 


—o 


Sm (w) e@* dw. (28) 


and also bearing (12) in mind, we obtain for pz (T ) the 


expression 
foe) 


P= | Y Sm(wetordo, (29) 


—oo m=0 


The functions Smif w) entering into (29) depend on the 
values m, and on, which are unknown beforehand. It is 
necessary to calculate them. We obtain the first equation 
for the determination of m, and 0, from (29): 

co co 
t= | of 5 Sm(w)do. 


—o m= 


(30) 


Since the mathematical expectation m, of the input 
process is constant, one may write the following equation 
for the mathematical expectation of the process z(t): 

m 


= 31 
m, 1 + kx (m,, ¢,) . ( ) 





oe ¢ &€ Ue 6 FP £8 £.e + 2S oe 


where ky is the amplification coefficient of the system, 
and X =my/™, is determined by (4). Equation (31) is 
a second equation for determining m, and o,. Similar 
equations are also obtained by the method of statistical 
linearization. For systems whose linear parts are 
described by differential equations of higher order, the 
following method of solution of (30) and (31) is usually 
employed. 

For the initial values of m, and o, one may take any 
quantities, let us say m, and O,. Setting them into the 
right sides of (30) and (31), we obtain m,; and Oz, in 
general different from mzg and Og. Then the values m,; 
and oz, obtained should again be substituted into the right 
sides of (30) and (31), after which the new mzgy and Og are 
obtained. This process is to be continued until upon a sub- 
stitution, results are obtained which coincide with zero- 
order accuracy with the results of the preceding substitu- 
tion. We note that the procedure of substitution consists 
in calculating the functions S,,(w) from (22) - (27) for 
given values of and o,, thereafter calculating (30), 
as well as(31). Practical calculation by the method of 
statistical linearization shows that a similar procedure of 
solving the equation allows a sufficiently rapid calculation 
of m, and o,. 
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After m, and oz are calculated, and also ny(mz, %z), 
N2(Mz,%z) and so on, the variance of the process y(t) at 
the output is easily calculated, by using (7) for r = 0, and 


oo 


substituting the weight function g( « of link 2 in the 
straight circuit, for k( o. Taking the Fourier transform, we 
obtain the following relation: 


oo 


+2 . 
oy= te | (Gio) |? {a8S, (0) +E ge | 52) 52(@—v) dy + 


—oo 
2 


—oo 


co 
where S, = > Sm(@) for the calculated values m, 
m=0 


and o.. 

We return to (22). It represents an equation for the 
calculation of the spectral density of the process at the 
input to the nonlinear element by the method of statistical 
linearization. Hence,the method of statistical lineariza- 
tion is obtained asthe first approximation to a more com- 
plete solution to the problem. The functions Sp{w) for 
m > 0 are further approximations. Since they contain the 
parameters n?/ 21, n?/ 3! etc.,which under definite condi- 
tions are negligible, the further approximations under these 
conditions carry but small contributions in forming the spec- 
tral density, and the method of statistical linearization gives 
the correct solution. However,there may occur conditions 
under which one may not neglect the further approxima - 
tions, and then the method of statistical linearization yields 
an invalid solution. This may occur, for example, when the 
spectrum of the signal x(t) is essentially concentrated in a 
neighborhood of the frequency u, and the linear link has 
a resonance at the frequency 2u or 3ug and so on. In or- 
der to clarify these conditions, it is necessary to carry out 
numerical calculations for different specific cases. Such 
a calculation may be accomplished simply with contem- 
porary computing machines. 


CONCLUSION 


In this paper it is shown that under the assumption on 
normal distribution of a random process at the input to a 
nonlinear element,the problem of calculating the disper- 
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+ 2 sae \ S, (v1) Sz (ve) Sz (@ — vy — v2) dv, dvg + .. +. 


—~0o 


(32) 


sion at the output of a system with feedback containing a 
nonlinear element may be solved more accurately than 
may be done using only the method of statistical lineariza- 
tion. 

The method of statistical linearization is obtained as 
the first approximation in the <pproach to the solution of 
the problem presented, and consequently, the possibility 
arises of being able to test its applicability in different 
cases. 

The calculations by the formulas given may be carried 


out with the aid of contemporary computational techniques. 


In conclusion,I take this opportunity to thank R. L. 
Stratonovich deeply for valuable counsel and direction 
given by him in connection with this paper. 
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AN INTERPOLATION METHOD FOR ANALYZING AUTOMATIC 
CONTROL SYSTEM ACCURACY UNDER RANDOM STIMULI 
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Translated from Avtomatika i Telemekhanika, Vol. 21, No. 4, pp. 481-488, April, 1960 


Original article submitted December 23, 1959 


An approximate method is suggested for constructing the solution of a system of ordinary differential equations 
containing random functions and random variables, and formulas are derived for the computation of the moments. 
The method is a development of certainideas of Chebyshev [1]. 


INTRODUCTION 


In designing or, when necessary, in improving an auto- 
matic control system, the problem arises of analyzing the 
accuracy of control when one takes into account random 
disturbances acting on the system. 

adY; 

a = F; (t, oe ae 

t=0, Y,(0)=y,, 
where the Y; are parameters characterizing the control pro- 
cess, X(t), X_(t),....X7(t) are random functions of time 
whose probability characteristics are assumed to be given, 
and Ay, Ag,-..,Ag are random variables whose probability 
characteristics are also assumed to be given. 

The initial conditions (t = 0, ¥;(0) = yj9) are also ran- 
dom variables with known probability characteristics. The 
functions K (i = 1,2,..,n) are such that there exists a unique 
solution of system (1) for each choice of random initial con- 
ditions, random variables A,,...,A, and random functions 
X1(t), Xg(t),..-.X7 (0). 

Ordinarily, to analyze accuracy of a control system, 
the nonlinear differential equations are replaced, with cer- 
tain assumptions, by linear ones, and then one uses the 
quite well developed methods for analyzing the accuracy 
of linear systems with random stimuli present [2,3]. Of all 
the existing linearization methods, the most practical in 
the given case is Kazakov's statistical linearization meth- 
od [4,5]. 

In [6,7], for the analysis of accuracy of nonlinear sys- 
tems with random stimuli, there were presented quite gen- 
eral theoretical methods which could serve as starting points 
for the construction of engineering methods and, in the sim- 
plest cases, for engineering calculations as well. Thus, in 
[7], Dostupov put forward an engineering method for an- 
alyzing the accuracy of nonlinear control systems with ran- 
dom stimuli, this method being a development of Pugachev's 
method which assumes the use of the latest computing tech- 
nology. 

In the present work we present a new approximate meth- 
od for analyziag the accuracy of an automatic control sys- 
tem with random stimuli, wherein the use of electronic com- 
puters is also assumed, but where the method is based on 


» Yn, Xz (t), Xe(t), . . +, Xr (t)s Ag, Ay, - - 


The systems of differential equations which describe 
the operation of an automatic control system with random 
stimuli are, in the general case, nonlinear: 


en hs). (1) 
(i az 1, 2,..2, 8). 

other principles than those of [6,7] and does not require that 

one solve an auxiliary system of algebraic equations, as is 

the case in [7]. 


1. Posing of the Problem 





If we use the canonical expansions of the random func- 
tions X;(t), Xq(t),.--.X, (t) and limit ourselves to a finite 
number of terms in these expansions, we may approximately 
replace system (1) of differential equations by a system of 
the form 

ay; 
= _ F(t, \,, Y2, *. 


t=0, Yi(0)= yw 


*? Ta A, he, ** * Am). (2) 


(i am 4, 2, . oe &)- 


which contains only random variables in the right members, 
thus significantly simplifying the investigation. 

All the assumptions made relative to (1) remain in force 
for this new system. Moreover, without reducing the gen- 
erality, one may assume that the initial conditions are not 
random, since otherwise one could attain such a position by 
a linear change of variables, which would only entail an in- 
crease in the number of random variables entering the right 
members of (2). 


The random variables Ay, Ag,...A,, will be assumed to 
be independent, which is also not a very rigid requirement, 
since in many practical cases one may attain this state by 
means of a linear transformation of these quantities. As the 
probability characteristics of the random variables, there are 
given the distribution densities Py(A4), Pp(Ag),--.Pm(A 1), for 
which there exist moments of every order, i.e., 


\N P5 (yj) diy = m", where the m{¥) are certain con- 


stants for any k and j (k = 0,1,2,..;j = 1,2,...,m). 
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Since, by assumption, there exists a unique solution of 
(2) for each choice of the random variables A,, AgseeA re 
which are independent parameters of this system, then the 
solution of system (2) can be given in the form of certain 


functions of time and the random variables A,, Ags-sA rp! 
Y,=Y;i(t,? » » Am) (i= 1,2,...,n). (9%) 


For the set of functions in (3), which are solutions of 
(2) for any values of the random variables A,,...,Am, We 
shall assume the existence of the second moments, i.e., 


M [Y?] - 


. rod [P10 


9 


phen cs oT 


“Ly he, ee 


\-.-\ YFG) hess a j)dh; = (4) 
A 


= K; 


(i=4 


where the K; are constants and A is the range of values of 
the random variables A,,. dn and is an m-dimensional 
finite or infinite parallelopiped. 
We choose some set of numbers { Hiky> Hake» +-ssmkm } 
(ky = 0,1,2,..., 5 Kp =0,1,2,...,q; Kp = 9, 1,2,..., q) 
which represents (q+ 1)™ points in region A. 


Let a known solution of (1) be computed at each of the 
aforementioned points: 


Yj (t, Min? Hox, = fore nk) 


(i= 4,2,...,0) (5) 


(By =O, 1, 2,....¢; Bee 0,1, 2,....G5 -2- 


k,, =0,1,2,...,4). 


If the probability in (7) tends to zero as q increases, one 
has convergence in probability, which is ordinarily written 
in the form of a limiting equality, 


lim P [Vi (t, Ad, ha, eeey hm) Pa 
qo (8) 
~Yi(t, Ay, de, -++yhm)| >8] = 0 (i=1,2,..., 0). 


Instead of the minimum probability principle, one may 
use the computationally simpler principle of the minimum 
of the second moment of the difference between the exact 
and the approximate solutions: 


M [(Yig — Yi)*] = 





m 
-{...( Pa —¥i)? [] 250s) 2; = min (9) 
A j=1 
(t==4,...,n). 
Indeed, according to Chebyshev's inequality, 
5 M ((¥,,—Y, 
P[\Yig—Yil>el< a 1 ao) 
It follows from Chebyshev's inequality that, if 
lim M [(¥4, — Y;)*] = 0, (11) 
q->-co 


then there will be convergence in probability. 

If the function Tig is chosen from the class of all pos- 
sible m-dimensional polynomials of degree q, then the pro- 
blem of minimizing the second moment is, in theory, easily 
solved. 


Indeed, by virtueof the assumptions as to the weight 
functions, for each function Pi(A;) (j = 1,2,..,m) there exists 


It is required, from the given set (5) of solutions of (2), [9] a system, unique to within a common factor, of ortho- 


to construct the functions 
Yig (i, Ad, he, ee Xm) (i 


which are sufficiently accurate approximations to the solu- 
tion and, from the functions in (6), to construct the simplest 
possible formulas for computing the solution's probability 
characteristics. 


(6) 


1,2).. 


oo R)> 


© 
“. 


_Probability Principle for Constructing 





Anproximate Solutions 





One of the most general principles for construction of 
the functions in (6) is the principle of minimum probability 
for each fixed q. 

In the case under consideration, this principle may be 
formu'ated as follows: from among the set of all the func- 
tions Pig (t,Ag, Ag, ++» Am)e constructed on the basis 
of the functions yj(t, W1ky» H2kg,..-» Hmkm), we choose 
a function such that 


P{\Yig—Yi(t, m1, > hm) | > ©] = min 


292) 


Qs. 


(7) 


(t=1,2,.. 


gonal polynomials Hjg,Hj1,Hp,...., i.e., such that 
b; 
j 
0 for isk 
H:(h;) He (h:) P: 0: — , 
) ji (Aj) Aj (Aj) P5 (5) dd; ‘ct for i=k, 
j 


where the numbers aj and bj can be infinite. 
Let 





m 
(F) 
Y= 2 Ginn, ... Ny, (t) I] H in; (Aj). (12) 
TMHNz ... Ny j=1 
where 
m 
“as ¥, Ua gn; (Ay) P5 (Ay) dh, 
Ainn, ... Nm (t) = P ‘i=! (13) 
- -\ lla in, (A;) P;(A;) dd; 
A j=1 
(ny, Na, ++. »n,, =0, 1, 2,...). 


From the theory of Fourier series we know that the 
minimum of the integral in (9) in the class of all possible 


for each fixed q, where € > 0 is any arbitrarily small given m-dimensional polynomials of degree q is furnished by seg- 


number; the symbol P[...] denotes the probability that the 
inequality within the square brackets actually holds. 
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- ments of the Fourier series for the functions ¥j(i = 1,2,.-0M) 
e., the polynomials 





one 


nay 
m 
ct 


0) 


12) 


(13) 


he 
ssible 
by seg 
2. . 0), 





én nm, =0,1,2,...,9¢ (14) 


\ny,= 0, 1,2,....9 

However, the problem we have posed can be considered 
solved only theoretically since, for the computation of the 
m-fold integrals entering into (14), there are only (q+ 1)™ 
points at which the values of the integrands may be assumed 
to be known. 

If the solution of (2) can be presented exactly in the 
form of m-dimensional polynomials in the random variables 
Ay, Age--Arn of degrees not exceeding the number q, then 
the problem as posed can be solved exactly by means of 
Lagrange interpolating polynomials. in 

In fact, in this case, from the given(q+ 1) interpo- 
lation nodes, one can construct, using the Lagrange inter- 
polation formula for the m-dimensional case [8], an m- 
dimensional polynomial of degree not exceeding q: 

Yiq = 


>) Yi (t, ikys okss sallintlhad » mkt) (15) 





Kyky,.. Rm 
mh 
[ 5. g-41 (A5) 
: : a9 


(= 1, 2,..., m) (ky = 0, 1, 2,..., 93 


ko = 0,1, 2,..., 93... 5 km=0,1,2,...,9). 


coinciding identically with the solution of (2), where 
Wj,q41 (Aj) os 

= (hj — pans) (Aj — Pos) + © (Aj — Hangs) j,a41 (H45x:) 

is the derivative, computed for Aj = H fky: 


However, the case more frequently arises in the appli- 
cations when the solution can be given in the form of m- 
dimensional polynomials of degree q only approximately, 
and the interpolation polynomials of (15), therefore, will 
in general be only approximations whose closeness, in some 
sense or other, to the solution of (2) depends on the choice 
of the interpolation nodes. 

Since the interpolation nodes can be chosen arbitrarily, 
it is then necessary toso manage this choice that the sec- 
ond moments of the difference between the solution and 
the approximation to it in the form of interpolation poly- 
nomials be, for each fixed q, as close as possible to 
their minimum values attainable in the class of all 
possible m-dimensional polynomials of degree q, and so 
that the approximation, in the form of interpolation poly- 
nomials, converge in probability to the exact solution of (2), 

If one chooses as the (q + 1)™ interpolation nodes all 
possible combinations of roots of the orthogonal polyno- 
mials Hig 41 (A.Xj=1,2,:--6M) of degree q + 1, which 
may always be done since the roots of these polynomials 
are real, different, and lie in the intervals of orthogonality, 


then these requirements can, to a known degree, be satisfied [9]. 
in this case, the interpolation polynomials which are 
the approximate solution may be written in the form 


Yig = >) Yi (t, Pikis Pokss ++» Lemke) 
kik... Km 
™m 


I Hy, q41 5) (16) 
jt A j,q-41 (jug) A; — Hye) 
(i=1,2,...,m) (ky =0,1,2,...,93 
ke = 0, 1, 2,00 Gee ee tky, =U, 1, 2,..-,9)§ 





Fi 5,941 (hj) = (Aj — pans) (Aj — pears) « » « (Aj — mes) 


H j,q41 (inj) is the derivative, computed for A, = By? 
the numbers y iky satisfy the algebraic equations 


A 5,941 (Hjx5) =0 (j =1, 2,...,m). 

It can be shown that the construction of the interpola- 
tion polynomial by (16) with the thus chosen interpolation 
nodes, which we shall for brevity henceforth call Chebyshev 
type nodes, is equivalent to the construction of the segment 
of the Fourier series from the (q + 1)™ first terms if, in com- 
puting the coefficients of the expansion, one uses a Gauss 
type formula for mechanical quadrature corresponding to 
weights equal to the given probability density functions. 
Since Gauss type mechanical quadrature formulas are the 
quadrature formulas with the highest degree of algebraic 
accuracy then, in general, it is impossible to give a set 
of (q+ 1)™ interpolation nodes from which one could con- 
struct an m-dimensional polynomial of degree q lying closer 
to the segment of Fourier series (14), from the point of view 
of minimizing the integral in (9), than the polynomial con- 
structed from the Chebyshev type nodes. 

As was established earlierfcf, (11)], to prove conver- 
gence in probability it suffices to establish a limiting equal- 
ity to zero of the second moment of the difference between 
the exact and the approximate solutions. 


In the given case, 


lim)... %—Yu)* il P;(,) d,; = 


— j=1 
(2. (V5 T] Pi) as — 

4 j=1 

—2him |. J ¥tig 
J] Ps Qa + lim |...) PE, [] Psd as = 0, 
j=1 ne - j=1 (17) 
if 
lim \---\ VAT] P; (hj) d; = J---) PET] Ps Qs drs 
Te a j=1 A j=1 = (18) 
and 
lim |... | ¥i¥e J] P53 (dj) ddy=].. | YET] Ps Qi) as. 
rr eo j=1 A j=l 

(19) 
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By substituting, in (18), the expanded expression for 
the interpolation polynomial, and by taking into account 
the orthogonality of the polynomials Hj q 4 4 (Ay), we 
obtain, after a number of transformations, 


lim \ “y \ Y?, [] 7A) @ = 
weer @ j= (20) 
= Tim DS) yi ts take Pokey +++ stm) TD 5x; 
tains ee a int 


(i= 1, 2,...,m) (ky = 0, 1, 2,....9; 





ka = 0,1, 2,000) Q3--03 hy =O, 1, 2,...,9)5 
© Heo. (hz) (P; (A) dd 
where Pj, = \ jeata (Aj) (P35 (A;) 2A; se 


a; As 9-41 (H5jx;) (Aj — jx) 


Christoffel numbers. 
Formula(20) is a Gauss type mechanical quadrature 
m 
formula with weights lI P;(d;) for computing the in- 
j=1 


j.- (YP [PP Oars, 


A j=1 


tegral 


It is well known®* that the mechanical quadrature 
process converges for any Riemann- integrated function if 
the domain of integration is bounded. If the region is un- 
bounded the same state of affairs can be attained if, in 
constructing the interpolation polynomial, one does not 
take into consideration nodes which lie at a sufficiently 
great distance from the origin of coordinates, a process 
which can only have a beneficial effect in calculations 
with electronic computers (since cases of overflow are 
avoided). 

Since (19) follows from (18), it can be asserted that, 
when the solution of (2) is such that the second moments 
in (4) exist, then the interpolation process constructed on 
Chebyshev type nodes tends in probability to the solution 
of this system. 


3. Calculated Formulas for Computing the 





Moments 


The approximate formulas for computing the moments 
are obtained from the ordinary formulas expressing these 
moments in general form if, in these formulas, the exact 
solution of (2) is replaced by approximate ones in the form 
of interpolation polynomials constructed from Chebyshev 
type nodes. 


We now obtain the formulas for computing the first- 
and second-order moments: 


M [Yi (t, a, do, - 2+ 5 dm) =[...(% ][ Pi Qs) AO = 
ay j=1 
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oy}. 


=~. +. [Via] [Pi 0a (i= 1, 2,.. 
A j=1 


M (Y: (t;, hy, he, eves Xm) Y; (ts, A, he, She, hm)] = 


=| is | Fig (ta ba, Be, » © opAm) ) (ta, Ria, he, . « » hm) 
. “< (21) 
j=1 
(i= 1, 2,...,0; f= 4, 2,0005n). 


By substituting in these formulas the expanded ex- 
pressions for Yjg and Yjg and carrying out a number of 
transformations, we get 


= 
’ Umitm ) ] | P jk; 


j=1 
’ hm)] = 


M [Yi] = >) Yi (t, Pars Poms ++ 


Kykz...kyn 


M [Yi (ty, Aa, Ae, ee » Am) Yx (te, A, ha, ee 


nt >» Yi (t, Prrkys Poke, » + - 


Keyke...km 


, mk) 


m™m 
Yk (to, ik, Poke» ees mk mn) I] Pjk; 


j=1 


ken, =0, 1, 2,--+59)- 
where the p jk, 2 Christoffel numbers, 
1 


The approximate formulas for computing higher-or- 
der moments have an analogous form. 

In computing the moments of the solution of (2) by 
(21), it is most convenient to determine the functions 
yi (te Hakys Hekgr--+s H mkm) (i = 1,2,...,n) by numerical 
methods, carrying out the computations on electronic di- 
gital computers, since the accuracy of analog computes 
is, in many cases, inadequate; moreover, the final com- 
putations by (21) equally require discrete (digital) tech- 
nology. The Chebyshev type nodes and the Christoffel 
numbers Pjk,° Pik,’ Pikm: are computed from the 
given probability density functions. However, in many 
practical cases this is not necessary, since one can use 
prepared tables [10] of roots of orthogonal polynomials 
and the Christoffel numbers corresponding to these poly- 
nomials. 

For example, let the random variable Aj be distributed 
(A; —%,)? 


normally, P;(4;) = = exp & _ We set 
A: 
j 


202 
Aj 


*From the Steklov-Polya theorem. 
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(21) 
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stributed 


We set 





= hj 4 %5,;Y2, The distribution density of the ran- 


dom variable x will be P(x) =e “* "SG. The Cheby- 
snev-Hermite polynomials 
1\" ., da” . 
H, (2) = - 3) ex* —_ [e—>*] 
z 
form a system of polynomials which are orthogonal on the 
interval (- 0, ©). 


+1 
Let x : oa (k= 0,1,2,...,q) be the roots of the 


Chebyshev- Hermite polynomial of degree q+ 1. In this 
case, the Chebyshev type interpolation nodes and the Chris- 
toffel numbers are computed by the formulas 


hin; on ¥ rs ato, V2 





(k =0, 4, 2,...,9), 
(22) 
: ( H,,, (2) P(x) dz tote) 
‘-, = - == 
my) Hy, (@) @— ae) Pk 


(k =k). 

For uniform distributions on the intervals (ay, b,], the 
Chebyshev type nodes and the Christoffel numbers are com- 
puted by the formulas 


a, a Me b, +a; 
pe ee ee Be B 
hie; = 5 Xk = 5 


(k= 0, 1, 2,..-,9), (23) 


1 


a ee 
Pik; = Pk = \ 
—1 


v4 Voss (x) dx 
4 (q+1) (q+1), ’ 
Uo+s (z°t )@—a2e™) 





iq+4)! att 
(2q + 2)! qxati 
the Legendre polynomial of degree q+ 1 and the 





where U44,(x)= (z?—1)2+1 iis 


(k= 0, 1, 2,..., q) are its roots. 


Since the method proposed here, for q = 1, gives a 
linear approximation to the solution of (2), it can be used 
in that case as one of the possible linearization methods 
which does not require that the equations of system (2) be 
transformed. 

CONCLUSIONS 


The method presented here can be used in the investi- 
gation of the accuracy of operation of an automatic control 
system described by a system of ordinary differential equa- 
tions for which, with given initial conditions and for any 
choice of the random variables entering into the system's 
right members, there exists a unique solution with finite 
second moments. Thus, the method can also be used for 
the investigation of systems with essentially nonlinear 
characteristics. 

In conclusion, the author wishes to express his deep 
appreciation for the friendly aid, valuable advice, and com- 
ments, which were advanced during the discussion of certain 
questions touched upon in the present paper, to F. P. Popov, 
V. 1. Zubov, I. Ya. Diner, A. A. Sveshnikov and A. D. Mak- 
simov. 


Example. An automatic control system is described by 
equations of the form 


+ Vt = (A+ 4), t=0, ¥(0)=0, 


where A is a random variable with zero mathematical ex- 
pectation, and dispersion Mi A*} = 0.01; 
a 


P(\)= e 008 





V 2x 0,1 

It is required to determine the mathematical expecta- 
tion and dispersion of the output coordinate Y. 

By using (21) and (22) and making use of tables [10], 
we find that, for q= 1, M[Y] = 0.949 and D[Y] = 0.0164, 
while for q *= 2, M[Y] = 0.964 and D[Y] = 0.0123. 

The exact solution, obtained by solving the given equa- 
tion with the given initial conditions, is that 


M[Y] = 0.963, D[Y] = 0.0123. 


Thus, in the given example, for q = 2 the approximate 
values of the mathematical expectation and the dispersion, 
computed by the interpolation method, coincide with the 
exact values to within one unit in the third decimal place. 
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By means of the theory of point transformations, an analysis is made of the operation of an extremal control 
system with a commutator. In the discussion, no essential limitations are imposed on the form of the nonlinearities. 
The existence of a limiting cycle is proven, and the character of its establishment is determined; the possibility of 
multiple oscillations, and the conditions for their generation, are investigated. The effect of asymmetry on the 
extremal function is determined, and an estimate of the accuracy of controller operation is made. 


Existing automatic control systems solve the problem 
of maintaining a controlled quantity either constant or 
else changing in accordance with a previously prescribed 
law (in the case of so-called “programmed” control). 

After a disturbing factor has acted, the controlled 
quantity changes its value, and the automatic controller 
must, by means of its action on the controlling parameter, 
take it to some given level in some interval of time. 

If the problem to be solved by the automatic controller 
consists in maintaining a given program, then the control 
system just sketched is theoretically valid for its solution. 
Thus, for example, we have the problems of maintaining 
the number of rotations of a machine's principal shaft, 
steam or gas pressure in a trunk, voltage or current in a cil 
cuit, the problem of maintaining a given ratio of compo- 
nents, etc. 

However, as technology developed, a number of pro- 
blems arose which differed in principle from those just 
cited, these being the problems of the automatic search 
for, and maintaining of, extremal values of the controlled 
quantity, either at its upper or its lower boundary. 

Such problems arise in the most diverse domains of 
technology (aviation, motor design, etc.). 

Naturally, such problems can be solved only approxi- 
mately by the ordinary control systems. 

Several years ago the author of the present paper de- 
veloped the automatic control principle which solves these 
problems, the principle which one may call that of extre- 
mal control. 

Some of the questions in the theory of operation of 
extremal controllers were already considered previously 
by the author [1,2]. 

In the present paper, using the methods of point trans- 
formation theory, we investigate the operation of one of 
the extremal control schemes. 

1. Principle of Operation of the Extremal 
Controller 

In view of the fact that the form of the functional re- 
lationship between the controlled and the controlling quan- 
tities is not known and, moreover, varies in the course of 
time by a completely unknown law, to give this relation- 
ship in advance is not possible. 
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{t is in relation to this that the idea of introducing an 
artificial disturbance is made basic to the principle of 
operation of the extremal controller. 

All the motions arising in the system as a result of 
this disturbance are divided into two classes: a) the class 
of admissible motions, b) the class of inadmissible mo- 
tions. 

To the class of admissible motions are assigned all 
those which move the controlled quantity closer to the 
required extremum or boundary, while the motions which 
take the controlled quantity further away are assigned to 
the class of inadmissible motions. 

By means of a special device, the direction of the 
process is verified, and if this direction causes the con- 
trolled quantity to approach the required extremum, then 
the disturbance is accepted; if the quantity moves away 
from the extremum, then the sign of the disturbance is 
reversed, and motion in the required direction then com- 
mences. 

To provide stability of control, the variations of the 
controlled quantity during the establishment process are 
carried out nonmonotonically. 

Basic to the scheme is a sr 2cial device which provides 
selectivity in choosing an admissible disturbance and which 
transfers the motion from the class of inadmissible motions 
to that of admissible ones. This device is a relay which 
takes into account the direction (sign) of the velocity of 
the disturbed motion. It is because of this that we shall 
henceforth call it a signum-relay. 

The signum-relay is made up of lever 2b (Fig. 1) 
friction-coupled to the axis of indicating instrument 1a 
which measures the controlled quantity p, so that the 
lever is dragged along by the pointer as p varies. 

However, lever 2b can move only within the limits 
of a certain zone between contacts 2c and 2d. This will 
be called the insensitivity zone. As the pointer of the in- 
strument moves to one side or the other, lever 2b will touch 
one of the contacts (2c or 2d) and will stop, while the point- 
er will continue its motion. 


* This work was first published in the Trudy TsIAM 


in 1949; it is reprinted here without change. 
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In addition to the signum-relay, the following basic 
elements enter into the extremal controller scheme (Fig. 
1): 

1. The instrument which measures the controlled quan 
tity. 

2. Reversible device 10 which is connected to the or- 
gan which acts on the controlling parameter. 

3. A contact device (a commutator) which provides 
the periodic character of the artificial disturbance. 

If the controlled quantity, whose extremum the con- 
troller seeks, depends on several parameters, then these 
latter are switched in by turn by means of device 4-5-6 
in such manner that, during the control time for one of the 
parameters, the remaining ones either remain free, or are 
kept constant, or are varied in accordance with a given pro- 
gram. 

Control is effected in the following manner. We assum 
that the controller is not switched in; then, the position of 
pointer 1a will not be changed (under the assumption that 
natural disturbances are absent), i.e., the parameter vy, act- 
ing on the quantity p, does not change. With this, let it be 
assumed that the quantity p does not attain a maximum. 


We now switch in knife switch 12. This starts the ro- 
tation of motor 3 (to one side or the other) and commuta- 
tor disk 7. 


With this, let the change of parameter v give rise to a 
decrease of p. With this, instrument pointer 1a, in rotating 
counterclockwise, drags along lever 2b which, after traver- 
sing part of the insensitivity zone, closes contacts 2d and 2 
which causes the operation of electromagnet 8 and stepping 
mechanism 9 which reverse motor 3 and change the direc- 
tion of parameter v. Thanks to this, the quantity p increas: 
es and, if the rotational speed of the disk of contact devic: 
7 is sufficiently small, the quantity p attains an extremum 
and begins to move about a maximum. 


To eliminate the possibility that the stimuli of disturb- 
ing forces would disrupt controller operation, contact de- 
vice 7(a commutator) was introduced. 





Let the permissible deviation of parameter v from the 
value vg which corresponds to a maximum be some Av. 
Then, the number of rotations of disk 7 is so chosen that, 
during one rotation, the change in the quantity v due to the 
action of motor 3 does not exceed Av. 

If, with such a choice for the number of rotations of 
disk 7, there now begins the action of disturbing forces, 
then two cases are possible: 

a) The effect of the disturbance on either side is great- 
er than the effect of variations in v. In this case, changes 
in v from increasing p cease when 43 ~ Ov = vSv,+ Av, 
where vy corresponds to the beginning of the disturbance. 

b) The effect of the disturbance on either side is less 
than the effect of variations in v. In this case, control pro- 
ceeds normally until one reaches the value of v* + Av (for 
disturbances which increase p) or v’ - Av (for the opposite 
case), where v® is the quantity attainable in the absence 
of disturbing forces. 

If the rotational speed of the commutators can not be 
considered small, and if a number of other design require- 
ments do not hold, then the value of v* attainable by the 
controller may differ from the optimal value vg. 


2. Investigation of the Establishment Pro- 





cess in the Extremal Controller 





An extremal controller is a special nonlinear oscillatory 
system with one, two,or more degrees of freedom as a func- 
tion of the character of the idealization. 

With a number of far-reaching simplifications, in many 
cases the motion of the controlled system in the absence of 
external disturbances, i.e., for fixed values of v{ in other 
words, with the extremal controller switched out), may be 
described by a linear differential equation 


d (1) 
= f(x) —p. 


We shall assume that the function fv) is continuous in 
the interval —a@ < y< +o and 
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/’ (vo) _ 0, e (v) < 0. 
We introduce the artificial disturbance by the law 


v= kt. 


(2) 
Equation (1) may then be written in the form 
dp 1 3 
a =~ l/)— Pl. = 
It is clear from (3) that, by choosing k, i.e., the velo- 
city of the artificial disturbance, sufficiently small, one 
can make the deviation of the controlled quantity from its 
Static value as small as desired. We assume thatk « 1; 
with this,we may consider that 


p= f(v). (4) 


With this condition, the extremal controller takes the 
form of a system with one degree of freedom. We denote 
by 6 the magnitude of the signum-relay’s insensitivity zone, 
and by At the commutator time (the size of the artificial 
disturbance's period). 

Then, Av = kAt is the quantity by which the argument 
v can change during the commutator time. We shall as- 
sume that, for the value v = vy, the commutator operates 
and reversal occurs, after which the argument begins chang- 
ing in one direction or the other. 

We now show that the entire line of initial values is 
divided into a number of regions, where this decomposition 
depends on the relationship of 6 and Av, and also on the 
sign of k. The character of the control process is essentially 
different for each of these regions. 

We shall later have occasion to use the following, al- 
most self-evident lemma: if the continuous function p=f(v) 
satisfies the condition f*(vg) = f**(v) < 0 for v = vg then, 
for each pair of values v, and vp such that v, <v}, =Vo, 
the following inequalities hold: 

a) if we are given two numbers C, and G (0 <C,=QG), 
then 


} (v») “a / (Y 7 C;) < } (a) x i (Va T= C2), (5) 
b) if we are given numbers d, and d, such that 6) 
f (ve) —f(ta— 4s) = f(%»)—flo—as), 
then 
d,<d, (6a) 
and 
Vg — dy < My — dy. (6b) 


The proofs of these inequalities are given in the Appen- 


" We shall consider the initial value v = vg, assuming 
that reversal to a decreasing argument occurs for the value 
v= v.. We assume that there exists a v' < vp such that 


f (v’) — fv’ — Av) = 8. (4) 


This can always be achieved by choosing a sufficient- 
ly small Av or a sufficiently large 6. 

We now show that, if (A) holds, the semiline of initial 
values vy = Vg can be partitioned into the two regions: 


340 





v <<, —oo<y<v’. 


In the region v’ < vy = vg the motion is stable, the 
amplitude being stationary for each initial value (i.e., 
there is a continuum of periodic motions). 

In the region vy < v', independently of which point of 
the region was chosen as the initial value, there occurs ap 
approximation to a limiting periodic motion which is car- 
ried out between the boundaries of v* - Av/2<y< y . 
where v® is the root of the equation 


fo)—f("' — & 

It was assumed earlier that, at the value v = vq, there 
occurred operation of the commutator and reverser, after 
which a decrease of v ensued. With this, the decrease of 
the argument can be no greater than Av: if the signum- 
relay operates, then the variation of the argument will be 
less than Av, while if it does not operate the variation of 
the argument will equal Av. 

As a consequence of the decrease of the argument, the 
quantity p begins to decrease. 

We now show that, if v,< v’ then, with a decrease of 
v by a certain definite amount a(0 < a < Av) given by(7) 
down to a value Ve = V2(V), the signum-relay must operate; 
the controlled quantity will thereby decrease by the amout 
6. 


= 6. (8) 


Indeed, the quantity a is defined by the condition that 
with a decrease of the argument v by a, the insensitivity 
zone be traversed, i.e., that p be decreased by 6: 


1 (v1) —/(%1—a) =8 (1) 


or 


f(%) —f (v2) = 8, (1) 


where Vz = vy ~ a. 
We use the lemma at this point. By comparing (7) 


with (A), we note that (6a) holds if we take 


%1=Ua, v=, a=d,, Av—d,; 
consequently, a < Av. 
At the moment of operation of the signum-relay, 


V = Vg = ¥y—a. 


Immediately after this, the argument begins to in- 
crease, continuing to do so until the next operation of the 
commutator. 

The increase of the argument is by the quantity 
Av -a = Av -(vq -¥g). At the moment of the next opera- 
tion af the commutator, the argument has reached the val 
ue 

Vg = V_ + [Av — (v, — vq)] = 2v,+ Av—v,. 

After this, the control cycle starts repeating. 

Thus, as the result of one complete control cycle, the 
point v, on the axis of abscissas is translated to the point 
and, correspondingly, the point p, of the curve p =f (v), 
with abscissa vy, is translated to point ps, with abscissa ¥- 
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It will be shown below that the image vg of point vy cannot 
leave the region v< v’*. 

Consequently, (8), where vg is defined by (7a), gives a 
transformation of the semiline v, < v’ into itself. We call 
this the transformation T. 

The invariant points, which remain unmoved under 
this transformation, are defined by the condition 


Vs = Vy. (9) 
On the basis of (9) and (8) we obtain 
Ve = Vy — ~ (10) 
or ae Av 
2 


If we eliminate a from (7) by means of this last equa- 
tion, and denote by v* an invariant point, we obtain the 
following formula for determining v*: 


j(e°)—1(v — 3) = 


8. The Establishment Process in the FExtre- 
mal Controller 


(B) 








We now discuss the question as to the character of the 
establishment process, as well as the number of invariant 
points which satisfy Eq. (B). 

After commutator operation, there begins the decrease 
of v down to the value of vg, in accordance with the equa- 
om Vg = Vy -a, (8a) 
where a is defined by (7). 

With the reaching of the value p = f (vg), the signum- 
relay operates, and the argument begins increasing up to 
the value 


Vy = 0, ++(Av—a) =v, —a+(Av—a). (8b) 


It follows from (8b) that the value of vs will exceed vy 
if the difference Av-a is greater than a; for the entire 
range of values of v, in which this is the case, transforma- 
tion T will translate every point of this region to the right. 

If, in some region of values of vy, the difference Av-a 


is less than a, then transformation T will here translate each . 


point of this region to the left. 


The condition Av-a > a leads to the inequality Av/2>a, 





a>>. 


Thus, transformation T shifts all the points of the semi- 
line v; < v’ which lie to the left of v* tothe right, and all 
points which lie to the right of v® to the left. It thus fol- 
lows that, on the semiline v<v* there is just the one im- 
movable point v*. 

In order to determine the character of the establish- 
ment process, we must Clear up still two more questions. 

1. Will the distances p; of the successive images vs, 
Vg,--- Of point vy from fixed point v* decrease, increase, 
or remain invariable ? 

2. Are there conditions under which transformation T 
will translate points of the semiline v < v’ lying to the 
left of v® into the region lying to the right of v* and, con- 
versely, can it translate points lying to the right of v* into 
the region lying to the left of v*? 

We first prove that, no matter what point v\< v* we 
choose, its image vg will not leave the region v < v’. 

Indeed, if the point lies to the right of v* and to the 
left of v’ then, as was proven earlier, transformation T 
translates it to the left; consequently, it can not move from 
the region v*< v,< Vv’ into the region v> v’ 


If, now, the point lies to the left of v* then, by com- 
paring (7) and (B), we see that, on the basis of (6b) of the 
lemma, vy -a< v" - -Av/2; consequently,’ 


vg = 0; —a-+ (Av —a) < v'— 4 + (Av —a) = (8) 


=v +2 —a.- 


But a > 0 and, therefore, by strengthening the inequality 


we get . 
vg <v +3. (12) 


We now compare (A) and (B). On the basis of (6b) of 
the lemma, we conclude that 


. A , . A , 
v— > <v'—Av or v+ ><’. (13) 


Thus, we always have vs < v’. 
We now prove that, with successive applications of 
transformation T, the successive images Vg, Vs,... always 


while the condition A v-a < a leads to the inequality Av/2 <a. #Pproach closer to v*. Let the distances of points 


We consider the two cases vy < v* and v*< vy < Vv’. 


The quantity a is defined by (7). We compare it with (B) 


and note that the right members of (7) and (B) are identical. 
Consequently, in the case vy < v*, the conditions of the 
lemma are met, if we take vy = Var V* = Vp, a = dy and 
Av/2 = 42 but, in this case, we have, according to (6), 


e<?. (11a) 


If we now consider the case when v “« vy < v’, then the 
conditions of the lemma hold under the assumptions that 


Av _ 


v=%, W=% a=d, =F d;. (11b) 


In this case, 


vy(i = 1,2,...) from v* be 


pi = v° — H. 


(14) 
We form the function a’: 
a’ =v — Fv, +4. (15) 
It is easily seen that 
p1—a’ = 4, (16) 


To determine the character of the establishment pro- 
cess, it is necessary to consider the two cases: p> 0 and 
Py <0. 





We first consider the case p; >0. In this case vy < v®. 
It was determined earlier that, for vy, < v*, Av/2 =a. 
Consequently, from (16), 


fi—a’>O0 (p>0). 
Moreover, it was shown earlier that, if 4 < ie’ then 
v* - Av/2 > v-a and, consequently, a’ > 0. 
Therefore, a’ can vary within the limits 
(17a) 


>a >0 


We now rewrite (8), taking (16) into account: 
Vs = 2; +2(F°—a) = 


=v°—p, + 2(p, —a’) =v +p, — 2a’. 


(8d) 


By bearing in mind that, on the basis of (14),for i = 3 
Vs ond v* — Ps» 
and by comparing (14a) and (8d), we obtain 
P3 = 2a’ — py. (C) 
By taking (17a) into account, it is simple to conclude 
from (C) that ps can vary within the limits -p< pg < pq. 
Whence 


lpsi<f1 (p1>0). 


We turn now to the case when p, < 0, whereby 
v' > vy> v*. In this case, Av/2-a < 0 and v* -Av/2< 
<-Vy -a. Consequently, a‘ < 0 and py -a’ <0. 

It follows from this that a* can vary within the limits 


’ 17b 
ai<a<0. —_ 
From whence, on the basis of (C), we conclude that 
18b 
lps|<ler| (<0). aun 


It follows from (18a) and (18b) that, for any values 

y4< Vv’, 
|ps|<| pr}. (18) 

Itis thus shown that no matter where the point v, lies 
on the semiline v <v’, we always have {p,| <|p,| i-e., 
repeated iteration of transformation T shifts its image ever 
closer to v*. In other words, for any value of vy < v’', the 
extremal control process approximates to the limiting mo- 
tion of period rT = 4 v/k, in which the argument changes 
periodically, with the same period T, within the limits 
v* - Av/2 = vsv*, and the corresponding controlled 
quantity p varies within the limits 


i(° — F)<p<f(’). 


Having (C), we can easily find the condition for which 
image vg of point vy turns out to be on the same side of the 
invariant point as vj. 

For this condition to be met, pg must keep the same 
sign as py. 

We consider, as previously, the two cases: p; > 0 and 
py <0. Letp;> 0. Since, with this, we also have a'>0 
then, from (C), it follows that pP,> 0 if 
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(14a) 


(18a) 


2a’ > 9- (19a) 


Now let p; < 0, for which a' < 0; consequently, ps<0 


'24|>|p1|. 
Both conditions (19a) and (19b) reduce to the one: 
=| 2a"|—jp,|>0. 7 


Thus, if the function ¥> 0, transformation T moves 
the image of point V, to the invariant point v* while keep- 
ing it on the same side of v®; if 


b<0, (20) 


then the image of point v, is shifted by transformation T to 


v* while simultaneously transferring over to the other side 
of v*. If 


if 
(19) 


¥ =0 (200) 


then image vs coincides with invariant point v®, and the 
further successive images remain immovable, 

Condition (20a) gives that property of the function 
p = f(v) for which py retains the sign of py. One may 
also give the boundaries for the point vy such that its 
image may not under any conditions coincide with v®, 
These boundaries are defined by the condition that ps# 0, 


Letp,; > 0, We have from (16) that 
(16a) 
, Av 
a=>p— > +a, 


On the basis of (C), the condition being sought reduces 
to the inequality 2a°> Py, Le., 


A ‘ 
2(0 > = +a)> A, 


from whence 
(21a) 
P1 > Av on 2a. 


But if vy < v*, then a can vary within the limits 
0<a< Av/2. Therefore, if py > Av, then (21a) always 
holds, With this, v;= v* — py=v*— Av. 

Now let p; < 0. With this, v'> vy > v* and the 
number a can vary within the limits 


a <a< Ao. 


The condition sought reduces to the inequality 2a°< 
< P, from whence,on the basis of (16a), 


py < Av — 2a, (21b) 


Condition (21b) will be met if P, < Av-2 Av=- Av. 

With this, 41=V° ~ Pi >v + Ay. Thus, if y 
is at a distance from v® whose modulus is greater than 
Av, then the image of point v, cannot, under any condi- 
tions, coincide with v*, 

We now consider the case when y, > v’. As before, 
let a decrease of v commence after commutator operation; 
the decrease of v can continue either until the signum- 
relay operates or, failing that, until the commutator opét- 
ates, 
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In the case under consideration (v, >v’ )the signum- 
relay does not operate, since the controlled quantity p 
can change, with a decrease of the argument by Av, only 
by the amount m, which {fs less than the insensitivity zone 
width 5, 

This quickly follows from the lemma, Indeed, 

m = f (vy) -f (vy -A v). 

From a comparison with (A) we note that v! < vy <V» 
and the variations of the argument are identical; conse- 
quently, on the basis of part "a" of the lemma, m < 4, 
Therefore, the reversal from an increase in v will occur, 
not by operation of the signum-relay, but by operation of 
the commutator, which will occur at the value v, = v,- Av, 

After commutator operation, the argument increases 
up to the value vs =(vy — Av)+ Av=y,. After 
this, there is again a reversal to decreasing v, operation 
of the signum-relay, but by operation of the commutator, 
since the signum-relay does not operate, reversal, etc, 

Consequently, one full control cycle translates each 
initial value which lies in the region v'< v < vp into it- 
self; there is a continuum of stationary amplitudes, 

In other words, if the control process starts from a 
value of v lying in the region v’ < v < vo, then there will 
occur periodic oscillations of period t = 2 Av/k for 
periodic variations of the argument within the limits 
y-ASvzevey and a corresponding variation of 
the controlled quantity within the limits of f (vy - Av) 
= p= f (vj). 

Until now we have considered all the points of the 
semiline v = vp except the point v = v*, defined by (A). 


All the values considered for v = vg correspond to 
definite actual oscillatory processes which occur in the 
extremal controller, This cannot be said of the initial 
value vy = v", 

In fact, according to (A), after a reversal to decreas- 
ing v occurs for the value v*, the controlled quantity 
changes by 5 for a change of argument of exactly Av, 


But with this there must simultaneously occur two 
actions; reversal due to signum-relay operation, and re- 
versal due to commutator operation, 


Such simultaneous operations have no real meaning, 
One might think that to this there would correspond a 
switching out of the entire reversing device; this is not 


admissible, despite the small probability that (A) holds, 
In practice, even if the initial value turns out to be equal 


to v’, one must provide operation either of the commuta- 
tor or of the signum-relay, Consequently, it would be pos- 
sible to assign the value of v’ to one of the regions, 


We now consider controller operation if the initial 
value is v> vo, If we assume that, after commutator 
operation initially, there occurs an increase of the argu- 
ment then, by assuming that there exists a v'* > vg such 
that 


{(v") — f (v" + Av) = 8, 


we can transfer over all the results obtained for v S vy to 
the case when v 2 vo, Instead of transformation T we 
shall have transformation T*, given by the formula 


v3 = 2v, — Av — »,, 
where v, is defined by the formula f (v, ) -f (v,)= 6. 

Consequently, in this case one can partition the 
semiline of initial values v = vp into two regions: vy = vy, < 
<v"* and v" < vy, The value of v"’ is defined as the root 
of (A‘), 

In the region vy= vy< v** the motion is stable, there 
being a continuum of periodic motions, 

In the region v’* < v, there occurs an approximation 
to a limiting perlodic motion, in which the argument 
varies periodically with period t = A v/k between the 
limits v’* = v Ss vy’ * +Av/2, where v’® is the root of 
the equation 


f(v")—f (0 + ZF) = 8, 

In order to exhaust completely all the possibilities 
which can ensue for values of Av and 5 for which there 
exist v= v" and v= v"*, it is necessary to discard the 
supposition that, at the moment when the controller is 
switched in, the commutator necessarily is operated by 
a decrease of the controlled quantity p, 

We note first of all, that if the initial value of vy, for 
which the commutator operates on an increase in the con- 
trolled quantity is at a distance from v» equal to or greater 
than Av, then this case leads to that already considered if 
one takes account, not of the value vy, but of the value 
vj = v, + Ov ("+ for vy < vp and *-" for v¥y> yp»), 

In fact, with this the following operation of the 
commutator will occur for a decrease of p=f(v), which 
is a case already considered; the behavior of the con- 
trolled system will thus depend on whether vj_ lies in 
the region v' < vy < v" or in the region vy< v’ and 
vyj> v". 

In the first case ,there is a continuum of periodic 
motions; in the second case, there are either periodic 
motions about v® and v*® if, respectively, vj = v® or 
vj = v**, or an approximation to these limiting motions, 

There still remains for consideration the case when the 
initial value of v, for which the commutator operates on 
an increase of the controlled quantity is separated from vp 
by a distance less than Av, Let 


vp = V9 t+b>V 


It is easily seen that, with a decrease of p from the 
maximum value p = f(v»), | which will occur with de- 
creasing v, down to the value p = f(v, - Sv), the sig- 
num-relay does not operate, since the change in p, dp = 
= f (vg) —f (ve — Av) will be less than 5, Con- 
sequently, the following commutator operation occurs for 
the value: 


V2 = % + b— Av< tp, 


(0<6< Ar). 


Vo — Vg = Av—b< Av. 
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After commutator operation, there begins an increase 
of v from the value vy to the value 


V3 = V2 aa Av = Vo -}- b —_- Av -} Av = U9 -}- b = VU. 


With the attainment of vs = vy, there occurs a new 
operation of the commutator, and the cycle begins to re- 
peat, Consequently, one full control cycle will transfer 
the initial value of vy into itself, i,e,, there is a continuum 
of periodic motions, The case when v; = vo-b, 0 < b < Av, 
reduces to what has just been considered if one begins the — 


investigation with the following operation of the commutator, 


Thus, the partitioning of the line of initial values, 


~a@ < Vy < +, into regions depends on which side the argu- 
ment begins to vary on at the moment when the controller is 


switched in, 
Case A, If commutator operation gives rise to a de- 
crease in v, then the vy line is divided into three regions: 


region ly 


m1 < 0’, (22*) 
region L, vo << v" + Av. (22"") 
region Ls v” -t Av < v1. (22"**) 


Case B, If commutator operation gives rise to an in- 
crease in v, the regions are shifted to the left, and are the 
following: 


region L{ v,<v' — Av, (23") 
region L} v—Av< vy < v". (23"*) 
region Lj v” <2. (23""*) 


In the first and third regions there occur approximations 
to the limiting periodic motions while, in the second region, 
there is a continuum of periodic motions, 


4. Effect of Variations of the Parameters 


Av and 6 





We have thus far assumed that 4 v and 5 have val- 
ues such that there exist numbers v'< v) and v""> vo, 

We now dropthis limitation and consider the in- 
fluence of changes in the parameters Av and 5, 

Let there exist real values of v' and v"", We shall fix 
Av = Ayy and we shall vary 5, 

The decomposition into regions Ly, Lg, and Ly depends 
on the character of the real roots of (A) and (A‘), 

It is easily seen that, as 5, increases, the points v’ and 
v'* move further away from vp while, for decreasing 5,, the 
points v’ and v’* move closer to vp. For some definite val- 
ues 6 = 6, and § = 5,, the points v’ and v"" coincide with 
Voi the values of these are determined, respectively, from 
the expressions 


1 (¥o) — f (Yo — Av) = 4. 
1 (Yo) — f (Yo + Av) = 8. 


(24) 


Simultaneously with the change of position of the 


points v’ and v"* on the line there is a change of position 
of the invariant points v* and v4, defined by (B) and (B*), 
These points move away from vp with increasing 5 and 
approach vp for decreasing 5, Suppose that they coincide 
with vy for the values 5 = 5, and 5 = 5,, defined by (25); 


‘I (0) —7 (vo a >) = 65. (25) 


1 (%o)— 7 (X + >) == 8. 


It is obvious that 5,< 5, and 5, <5,, We consider the 
case when the curve p = f(v) is symmetric about the line 
v=v> Then, 5, = 5, and 5, = 5,, Let 5,<5 < 5), Con- 
sequently, there exists no v’ which satisfies (A). Let (vp) - 
-f (vo -a,)= 5 and consequently, a, < Av, There: 
fore, operation of the signum-relay occurs to the left of v, 
if the initial value lies to the right of vp in the region ly 
[Condition (26) ]. 

With this, the regions of (22) are deformed, and assume 
the forms: 


region ly 
Vv; << v9 + Av—a. (26") 
region J 
. oy pe Pe ed (26"*) 
region Lg (26 


V+ a < 7%. 


We shall continue to decrease 5, If 5 = 5g, i.e,, if 
v* coincides with vp, then a, = Av /2; then, 

, mtann te. 

Consequently, region L, shrinks to a point, There re- 
main two regions, Region Ly : vy < Vo + 4 v/2; region Ly 
vp tA v/2<vy. 

With still greater decrease of 5, the motion loses its 
simple character and becomes a multiple motion, 

It may be shown that the mode of multiple oscillations 
is not effective in practice here, since oscillations are pos- 
sible with this that have greater amplitudes than the simple 
oscillations, 

It is therefore meaningful to decrease 6 and to increas 
Av until v® gets close to vp, but not beyond that; this is the 
more advantageous since, generally speaking, {t is quite dif- 
ficult to obtain a small 6, but a decrease in Ay is very de- 
sirable from the point of view of increasing the stability of 
the control mode with respect to external disturbances, 

We have been changing § with Av held fixed, It is 
easily seen that an increase in Av acts analogously to & 
decrease in 5, 

We can consider the case of asymmetry analogously, 


Vo + Av — a, = % + 


5. Consideration of a Concrete Problem 





Let the function p = f(v) have the form 


= —p?, (27) 
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We shall consider initially a value of vy > Vo = 0. 
If Ay and 6 are given, then the magnitude of v® is defined 


by the equation } ‘ Av 
v= vI—_—- —)|) = $ 
j(v°) — f(x" — >) ~ 


from which, based on (27), 
* Av 7) ° 


Equation (28) is meaningful when 6 = Av*/4. Region 
(A) is isolated by the boundaries v" and vp = 0, where v’, 
in accordance with (A), is defined by the formula 
as gts Orage: 


en 29 
an 2Av (28) 


We now determine the character of the establishment 
process in accordance with (20a) and (20b), The character 
of the process we are seeking depends on the sign of the 
expression » = | 2a*| -| P, |. We now construct this func- 
tion, According to (14), 


, . Av 
a= (» — =)—G—e), 
where a is the root of (7), 
From (7) we obtain the quadratic equation in a: 


from whence 





a=1,+YV +8. 


(The plus sign {s chosen so that a > 0,) 





Consequently, 
($8) fo. Vw 
Further, 


. Av ti 
. Tending —=(F—Z-n). 


Inequality (20a) takes the form 
3 3 - . # 
p=+4[—fav—L4442Vei+ a] >0. 


(The plus sign is taken for v, < v ®, the minus sign for 


¥, > v® ; this expression is meaningful for 5 = 4 v*/4. 
It is easily shown that this expression cannot be satis- 


fied for any arbitrary vy. It follows from this that it is im- 
possible to choose values of Av and 6 such that the approxi- 
mation to v® occurs for any vy, on one side only, However, 
for certain regions of values of Av and 5, one can cite such 
& point vy, = vy # v® whose image vgs coincides with the in- 
variant point v®, 

In fact, it follows from (20c) that this point is defined 
by the condition that # = 0, 

Whence, 





The minus sign defines the invariant point v®, The 
plus sign gives the point v,, sought: 


4 


». (30) 
3 Av’ 


3 
Vu - » Av -{- 


while it foliows from (28) that v* has a minus sign [ the 
condition for applicability of (28)] if 45 >Av’, 


region 


regions of initial values given by (23) and (26), So long 
as there exist values v' < vp = 0, Le., so long as 5> AY 
[this follows from (29)], the line of initial values is par- 
titioned into the following three regions (22): 


then the line of initial values is partitioned into the regions 
of (26) [where the quantity a, is defined by the condition 
f( v9) — f(vg — ay) = 5, but vp = 0, and, consequently, a, = 


- With further decreases in 5, the motion becomes multiple, 








































Formula (30) gives a negative value of v,, for b< 2,25Ay', 


Consequently, the values of v,, are meaningful in the 


0.25Av? < 8 < 2.25Av?. (31) 
We now consider the entire v, line, and construct the 


1) region ly 
Av = 
ake 
2) region L, 
Ar 8 Av La 8 Av 


7 
Ime << Tt rKae tT 


3) region Ls 7 As 
Zao + 3 S72 


If 5 is included in the interval 0.25 Av? <5 <A y*, 


= V6}: 
1) region ly 
Vv) << ¥+ Av— p33; 
2) region Ly 
vo + Av—Vb << m+ V3; 


3) region Ls yi 
Vo + <7. 


As 5 decreases from 5=6,=Ayv to 6 = 0,25Av', 
region L, decreases and, for § = 0,25 Av, shrinks to a point, 


Returning to (31), we note that, in the region of velues 
5> 2,25 Ay’, there occurs an asymptotic approach to per- 
iodic motions (if v, is chosen in regions Ly or Lg); in the 
region of (31) the establishment process can be either a- 
symptotic or finite; in the region of 5<0,25AV, there are 
multiple oscillations, 

The character of the possible motions in the system 
under consideration depends on the two essential parameters 
Av and 4, 

Figure 2 shows the A v, 5 parameter plane, 

In the region bounded by the arc of the parabola Av*= 6 
and the positive 6 semiaxis, there is a rigid excitation 
mode with a decomposition of the imitial values in accord- 
ance with the domains of (22). 

In the region between the arcs of the parabolas AV = 
=46 and Av = 6, there are rigid modes with decomposi- 
tion in accordance with the regions of (26). 
































eee 
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Fig. 2. 1 is the multiple oscillatory 
mode, 2 is the rigid mode, with de- 
composition into the regions of (26), 3 is 
the asymptotic approach to a periodic 
motion (the decomposition into the re- 
gions of (22)), 4 is the finite and asymp- 
totic establishment processes (decompo- 
sition into the regions of (22)). 








Fig. 3. 


The asymptotic approximation to periodic motion oc- 
curs in the region bounded by the arc of the parabola 
9AV = 46 and the positive 5 semiaxis. 


The region Ay? > 46 is characterized by the presence 


of a multiple oscillatory mode. 

Figure 3 shows the phase diagram of the multiple per- 
iodic oscillations for Av = 3 and 6 = 1. 

Figure 4 shows the development of the multiple per- 
fodic oscillations in time (curve I), where the axis of ab- 
scissas is laid out with time in the scale of 1/k. This 
same figure gives the development (curve IJ) in time of 
the simple periodic oscillations for the limiting value of 
Av = 2 for the chosen value of 6 = 1. It is clear, from a 
comparison of curves I and II, that the multiple oscilla- 
tions are less suitable, since their amplitudes are greater 
than those of the simple oscillations. 

In conclusion, we consider the question of the in- 
fluence of the system's parameters on the inaccuracy of 
the extremal control operation. By “inaccuracy of opera- 
tion” we shall understand the difference between the val- 
ue vg of the argument corresponding to the true extremal 
value of the controlled quantity and the limiting value 
v* which can be attained as a result of the extremal con- 
troller’s action. 
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We denote this difference by x: 


- — Av 32 
xam—o ab Ae, (32) 

Expression (32) is meaningful so long as x = 0, since 
the case x < 0 corresponds to the case when the motion 
becomes multiple,and when (28) can therefore not be ap- 
plied. 

Figure 5 gives graphs of the quantity x as a function 
of Av for different fixed values of 6. It is clear from this 
figure that each value of 6 has its corresponding value of 
Av for which the inaccuracy x vanishes. With decreasing 
Av the inaccuracy increases, this increase not being very 
large for values of Av close to Av; for significant de- 
creases in Av the magnitude of x begins to increase more 
rapidly since, in (32), the first, hyperbolic-type, term be- 
gins to play a predominant role. 


Figure 6 gives graphs of the quantity x as a function 
of 5 for different fixed values of Av. ‘This function is 
given by a family of straight lines with slopes greater than 
zero. 

A comparison of Figs. 5 and 6 shows that the increases 
in 6 have a lesser effect on the increase of the inaccuracy 
than do decreases in Av. 
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Appendix 
Lemma. If the continuous function p = f(v) satisfies 
the conditions f*(vg) = 0, f'*(v) < 0 for v = vg then, for 
each pair of values vg and vp (V,< Vp = Vo), the following 
inequalities hold: 
a) if there are given two numbers c, and o,(0< Cy SC»), 


f (v%) —f (% — 1) < f(%q)— 1 %a—er)s (5) 


b) if there are given two numbers d, and d, such that 


then 


f(%,)—fl%e— a) =f (eP—F(—ar) 
then dy < dy (6a) 
and % _ d < _ ds. (6b) 


Proof. We first note that, since f'"(v) < 0 then, for 
decreasing vy, the function f*(v)increases monotonically. 

First, we prove (5). 

a) We assume initially that vg < v},~-C,; in this case, 
the intervals v,-Cy Sv =v, and v,-c, = v = Vp are non- 


intersecting. 
With this, 
[7 ()lr<v,> [f’ ()losv,— ee (33) 
Further, D 
f (%) — f (% — 1) = APy = ( P(v)dv, — (gaay 
MpD—-C1 
f (%q) — 1% — 2) = OP, = ode I (v) dv. (34b) 


On the basis of (33), and bearing in mind that c, = cy, 
we conclude that AP, < APy. 





b) We now assume that vg > Vp-Cg. In thiscase, the 
intervals vg-Cg S v Sv, and v)-cy S v S vp, partially in- 
tersect each other; with this, the segment vp-c; = v 5 v, 
is common to both intervals. Consequently, on this seg- 
ment, the changes in the function f(v) will also be identi- 
cal. 


It remains to compare the changes of the function f(v) 
on the intervals vg -Cg S v S v)~-cy and vg Sv Sv. But, 
for these intervals, we have case"a". Consequently, here 
also AP,< AP,. Thus, (5) is proven. 


Now we prove (6a). 


We fix the quantity c, in the left member of (5), intro- 
ducing the notation cy = dg. We now, in the right member 
of (5), change cy (denoting it by d,) in such fashion that the 
inequality becomes an equality. 


Since, in (34b), f"(v) retained its sign over the entire 
interval of integration, the quantity defined by the integral 
in (34b) can be decreased only by decreasing the length of 
this interval, whereby, on the basis of (33), equality can be 
reached only for dy < dg. 

As for (6b), it is obvious. 
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The equation of motion is derived for a rotary type electromagnetic control element, and a comparison is 


made with experimental results. 


In automatic control and regulation systems, extensive 


The impedance torque is made up of the load torque 


use has been made of electric controlling elements in which M_,, the friction torque in the bearings Mp and the dampiy 


electrical quantities are transformed to mechanical ones. 
Such elements (electromagnetic and electrodynamic) are 
used in remote control systems, in electrohydraulic servo- 


mechanisms and in other devices of automation technology. 


Despire the widespread usage of electric controlling 


elements, their dynamics are still insufficiently explicated. 
The existing works [1] in this area treat only individual ques- 


tions, and lack a general character. 
In the present paper we consider the dynamics of an 


electromagnetic polarized controlling element, since it is 


the most widely used type. 


2 
i 























Fig. 1 


A schematic drawing of the electromaguetic control 
element is shown in Fig. 1, where 1 is the excitation coil 


(superposed magnetization), 2 is the control coil, 3 is the 
magnetic circuit, 4 is the armature,and 5 is the output shaft. 
When a controlling signal is applied to coil 2, output shaft 


5 is rotated by a definite angle. 


We now write the general form of the equation of mo- 


tion for the electromagnetic control element's armature: 


J§= My—M,. @) 


Here, J is the armature’s moment of inertia, 6 is the 
armature’s angle of rotation, Mg is the driving torque,and 


M; is the impedance torque of the armatures motion. 
The driving torque, equal to 
Mg = PR (2) 


(where P is the force engendering the armature rotation and 
R is the radius of application of the force), depends on the 


current flowing in the control coil. 
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torque M,. 

In considering the operation of an electromagnetic « 
trol element in any actual control scheme, it is necessary 
to take into account both the form of the load and the fri 
tional force in the bearings arising when an external fore 
is applied to the armature axis [2]. For investigating the 
dynamics of an inherently electromagnetic control elemez 
removed from any control scheme (and these may be ven 
diverse), there is the possibility of not taking the load on 
the control element's armature into account. In this cas, 
the impedance torque will be determined by the friction’ 
the bearings under the action of the armature's own weigh 
and by the damping which arises when the armature move 
within the surrounding air, i.e., 


M ;= Qrf sign § + kq6, (3) 


where Q is the armature'’s weight, r is the radius of the 


armature shaft, f is the coefficient of friction in the bea 


ings, k is a proportionality factor,and q is the coefficient 
of viscous friction. 7 
. By substituting (2) and (3) in(1), we get 
J# = PR— Qr/ sign _ kq@. (4) 

The work expended in moving the control element's 
armature, as in every polarized electromagnetic system, 
is carried out at the expense of two energy sources: the o 
trolling (working) and the polarizing sources. In the gene 
case, the force acting on the armature of the system unde 
consideration can be presented in the form of a sum of ti 
Soeces (3) P= P,+ P+ Ps, 
where P, is the component due to the polarizing source(# 
sources), which exists whether or not there is current int 
control winding, P, is the component due to the controllil 
(working) source (or sources), Ps is the component arising 
from the interaction of the magnetic fluxes of the contr! 
ling and polarizing sources which actually determines tht 
direction of armature motion. 

For the case when the polarizing and working source 
are coils, these components equal 
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where ie and i, are the currents flowing in the excitation 
and control coils, respectively, We and W, are the number 

of turns in the excitation and control coils, Ge is the cir- 
cuit’s conductance with respect to the magnetic flux of the 
excitation coil, Ge is the circuit conductance with respect 
to the control coil’s magnetic flux, Ge¢ is the circuit con- 
ductance determining the magnitude of the portion of the 
excitation coil’s magnetic flux in the control winding, y 

is the linear deviation of the armature from the neutral 
position *. 

As is obvious from these expressions, only the sign of 
the third component is determined by the polarity of the 
current in the control winding. The summation signs in all 
these expressions are indicative of the fact that polarized 
systems, as a rule, have complicated configurations of mag- 


netic circuits with several controlling and polarizing sources. 


The system under consideration is a differential one. 
In it (Fig. 1) there are two excitation coils on the arms of 


the magnetic circuit, and one control coil on the cross piece. 


With no controlling signal present, the electromagnet's 
armature is acted upon by the forces generated by the first 
and second excitation coils. In this case, the component 
P, is the sum of these two forces, and equals 
yo (iyw1)? dG, 4 (igwe)? dG, ; 
2 dy 2 dy 


(5) 





where G, and Gy are the circuit's conductances with respect 
to the magnetic fluxes of the first and second excitation 
coils, respectively, iy, ig and wy, Wg are the currents and 
numbers of turns in the first and second excitation coils, 
respectively. 

Since the coils are identical (wy = Wg = We) and are 
connected in series (iy =i, =i,), the expression may be re- 
written in the form 

i ave)" /dG 
p,=ee (+). 

The second component P, generated by the control 
coil, can be expressed by the formula 
__ (iw)? dG 

2 dy’ 
where i and w are the current and number of turns in the 
control coil, G is the circuit conductance with respect to 
the magnetic flux of the control coil. 

The third component, created by the interaction of the 
two counter-directed magnetic fluxes of the first and sec- 


ond excitation coils with the magnetic flux of the control 
coil, equals 


(6) 


(7) 





Ps 


_ = 


iwi QW » ( dGoe = dG, c) 
dy dy /’ 


(8) 


where Gy, and Gg, are the conductances which determine 
the magnitude of the magnetic fluxes in the control coil 
which correspond to the first and second excitation coils. 
In order to express Py, Pg, and Ps in terms of the con- 
tolling element's parameters, we use its replacement cir- 
cuit(Fig. 2). Figures 3 and 4 show the nomenclature of the 
basic design parameters of the control element. With this 












nomenclature, we can write that the impedances of the left 
and right air gaps equal, respectively, 


yey °) 
Rn=Gpye 











Fig. 2. 1 is the excitation winding’s mmf, 2 is the control 
winding’s mmf. 








Fig. 3. Basic design parameters of the electromagnetic 
control element. 


| 





Fig. 4. Design parameters of the armature’s bearing 
portion. 


and that the air gap between the bearing and the armature 


equals 


al 


where c = b- a, 2 = hk,, ky = 6/4, d is the armature thick- 
ness. 
We take R, = op, and calculate a number of auxiliary 


quantities: 2cd 
Rr + Res = ya: 

_ __ 2y% 
Rn — Bra = ayia’ 


* We consider small angular deviations of the armature, 
within the limits of + 1 degree. 
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7,2 





Rp hyo = 7s 


(c? — y?) a? : 


Ry Ryo + Rs (Rp, + Rps) = ae 


ty 4 1\ 

mnt = F\ tT): 
é 1 1 

Mim + B= F\sy57 7) 











R at Ry, Rpg + Wty (A; 1 tro) 
Se a Rp, + fs 
‘ (11) 
1 Rp, Rpg + Rg (Rp, + Ryo) 
R, = —- aed Ti, : 
The conductances will be, correspondingly, 
G Rp, + Rs _ (l+c¢-)-y)(c—y)d 
41 —at ea = ww) ee SF tee Ce die: 8 ro a 
FE. Me, 4 Maton & Bin (I 4- zc)? 
Pi ** pe 8 (4fp. Ps) (12) 
5 leat NE Te _(l+e—y)(e+yd 
~ Rp Rpg + Rg (Rpy + Rpe) (I+ 2c) 
and their derivatives are 
dG, __ I 4- 2y d dG, _ '—2y d (13) 
dys bc 3° dy Il-—-zei’* 
Consequently, 
% y - 2 (ie w e) @ 
P, = -—- ——-—_ . (14) 
(t+ 2c) 6 * 


The minus sign shows that component P, is directed 
opposite to the translation y . 

Component P, = 0 by virtue of the fact that the conduc- 
tance is want since 
— Ry + Pitre 


fae» > (15) 
sf Rp, + Rpg ¢ 


It 2ed 


o~ 27 - Rs " 
and its derivative dG/dy equals zero. 

To determine component Ps we should find the ex- 
pressions for Gy, and Ggc. For this, we sketch the replace- 
ment circuit (Fig. 5) and determine which portion of the 
fluxes engendered by each excitation coil will flow through 
the control coil. The action of each excitation coil will be 
considered separately. Based on Ohm's laws and Kirchhoff's 
laws, we get 





G are 
16 Rp, Rps + Rg (Rp, + pe)’ 
R (16) 
Goe =T nt. 
cep, 4 Rp + Rs. 5 (tp, + nf + Rpg)” 
Consequently, 
Rp, — Rp, 





Rp, Bre + Bs (Ry + Apa)’ 


After substituting of the impedance values, we get 


2dly (18) 


B(L-+ Be) ’ 


Go a Gio = 














<C) . 


whence 
“Gc Ge) _ _2ld _ (19) 
dy B (LU -- 2c) ° 
Consequently, 
- ld 
) , . * a ere ‘ 
Ps = iWigte ae (29) 


Thus, the general expression for the force will have the 
form 


P= Ps,+ P, (21) 


or 
ld (é = - 
a (l- + 2c) o (l- {~ % 


as (22) 





iw iewe ; 


As is clear from this expression, the component P, of 
the force on the armature acts as a spring, rotating the arm 
ture after cessation of the action of the controlling signal to 
the original position. It may therefore be called an electri- 
cal spring. 

Figure 6 shows the family of static characteristics 
M = £6) for various currents i in the control winding, these 
curves having been obtained experimentally. These graphs 
substantiate the correctness of our computations. 
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Fig. 6. Static characteristics of the electromagnetic con- 
trol element. 


For small angles of rotor deviation, we may assu.ne tha! 





y = OR. (23) 
Then, 
p. — iwigwld = 2 (i,w9* Rd (24) 
=F + 2)” ~8(L+ 2) 
or 

P = Ai— Bo, (25) 

where wigueld 2 (igwe)*® Rd 

. —~ 8(. + 2c)” ~ 8(L + 2c) 
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in (4), we obtain 


or 


8m and ip, are the basic magnitudes of the angle of devia- 
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an’ +aAr 


kq 


m= 


+s 


Be, Qir 
A ' AR 


J) = (Ai — BO) R — kh —OQrf sign 9 
sign 6 = i. 
J% (Ll 4- 2c) 8,, 


igupdlw hi, 
Or fb (1 -+- 2c) 0, 
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By substituting the expression obtained for the force P 
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form: 


me + y+ 8p +asigng =/, 


where 
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Fig. 7. Experimental amplitude-phase characteristics. 


tion and the control current. 
In (26), the term a sign ¢, characterizing the friction 


in the bearings, is a negligibly small quantity, since the 
armature axis is supported in ball bearings and operates with 
flood lubrication. In view of this, the motion of the control- 
ling device's armature can be described by an equation of the 


form 


where 


Expression (27) is the usual form for writing the equation 


Tz=m/B, Ty=%/8. 


(27) 


of a second-order oscillatory link. Its transfer function can 
be given in the following form: 


The form of the frequency characteristic of such a link 
is widely known, but the data which characterize the fre- 


1 





W (p) = 
(Pp) 3 


(T2p* 4- Tip + 1) 


(28) 
and the quantity ¥ is the lag phase. 





N= zt 


6 





i 


™ 





By going over to dimensionless symbols, we get the 
equation describing the armature motion in the following 


quency characteristic of an actual electromagnetic control- 
ling device are of definite interest. 

The amplitude and phase frequency characteristics of a 
type REP electromagnetic controlling device were recorded 
for a controlling current of i = 2 ma and for different inertial 


loads on the armature axis (for different moments of iner- 
tia)t. The armature's inherent moment of inertia] equalled 
60.5 g -cm*. The table presents the data obtained for 

a control current of i = 2 ma; the basic magnitude of arma- 
ture angular rotation was taken to be 6m =0.4°. The basic 
value of control current was taken to be i,, = 2 ima. 

The quantity N in the table is the ratio of the dimen- 
sionless amplitude of armature angle of rotation to the di- 
mensionless control current, i.e., 


*The experiments were carried out by L. P. Levin, G. Yu. 
Chubarova, and G. M. Val'kova. 









































Moment of inertia g - cm* 
Fre - 60,5 117 952 
aoe Sadeg N | vy, deg | %adeg | N | deg] %, deg) wn | y, deg. 
16 0.376 | 0.94 13 0.44 1.11 15 0.64 1.6 24 
- 20 0.372 | 0.93 20 0.48 1.2 30 0.32 0.8 35 
30 0.376 | 0.94 30 0.53 1.33 42 0.124 |0.31 55 
40 0.40 1.0 46 0.53 1.33 75 0.05 0.126 85 
50 0.42 1.06 60 0.32 |0.8 90 0.04 {0.1 108 
60 0.44 1.41 70 0.16 0.4 110 0.028 | 0.07 120 
80 0.396 | 0.99 90 0.08 0.2 140 0.02 0.05 150 
100 0.20 0.5 110 0.05 0.126 160 0.014 |0.035 180 
150 0.12 0.3 140 0.025 |0.063 




















The amplitude-phase characteristics given in Fig. 7 


were constructed from the tabular data. It is clear from 
these characteristics that the application of an external 


2. 
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load to the shaft of the controlling element's armature de- 
forms the amplitude-phase characteristic by way of in- 
creasing the phase of the lag. 
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The paper considers methods of differentiating slowly varying electrical signals with respect to time, taking 
into account the nonlinearities inherent in actual devices, as well as methods of replacing the differentiation 
operation by a determination of the increments of the input quantity during fixed intervals of time. 


posing of the Problem 





In the present work we consider methods for differentia - 
ting slowly varying signals with spectra from zero up to sev- 
eral cps. In such a frequency spectrum the signal derivative 
will be small, and the receiving devices will have finite 
sensitivities. Therefore, the differentiator must provide am- 
plification of the signal derivative. On the other hand, for 
slowly varying signals, the inertia of the differentiator it- 
self has a lower value. 

It is convenient to open the investigation of possibilities 
of differentiating slowly varying signals by comparing the 
frequency characteristics of the corresponding differentia - 
ting device with the frequency characteristics of an ideal 
differentiator 


Dig GQ) = Rid jQ, 


where Tg jq is the derivative coefficient (time constant) of 
the ideal differentiator and Q is circular frequency. 

On the basis of such a comparison one can establish, 
for the real differentiator, a working range of frequency 
characteristic in which its distortion will not exceed some 
established quantity. We adopt the convention in this work 
that the difference of the phase frequency characteristics 
of the ideal and the real differentiators 4 ¢ shall not ex- 
ceed 20°, and that the ratio of the ordinates of the ampli- 
tude frequency characteristics lie within the limits of 0.8 
to1.2. The derivative coefficient T, of the real differen- 
tiator is bounded. Therefore, in addition to the working 
range of the frequency characteristic, the magnitude of Ty 
for the given differentiator is determined. 


All differentiating devices can be divided into two types: 
differentiators in which the operation of differentiation is 
modeled (analog differentiators) and differentiators in which 
the operation of differentiation is replaced by the determi- 
nation of the increment of the input quantity during fir ‘te 
intervals of time (we agree to call such devices discrete 
differentiators). For the differentiation of slowly varying 
signals, discrete differentiators turn out to be competitive 
with analog differentiators. 


Analog Differentiators 





In this section we shall consider two basic types of elec- 
tonic differentiators, designed for the differentiation of dc 





signals, and an electric differentiator in the form of a servo- 
system with proportional velocity control. The basic ele- 
ments of electronic de differentiators are RC circuits* and 
amplifiers. After determining the working range of the fre- 
quency characteristics in the linear treatment of these dif- 
ferentiators, one must take into account the effect of non- 
linearities in the static characteristics of amplifiers and 
motors on the corresponding frequency characteristics. 

Differentiation of ac signal envelopes by means of pas- 
sive circuits is, due to large errors, efficacious only for rap- 
idly varying signals (cf., for example, [1]). 


A. Differentiator in the Form of a Differen- 





tiating RC Circuit with an Amplifier Con- 





nected at Its Output? 





The transfer function of an inherently differentiating 
circuit with account taken of the impedance of the leakage 
condenser is written in the form 

Tour (p) 2. TP+ RI Me 

Uin (p) Tpe+i+k/Ry 
where T = RC is the time constant (coinciding with the deri- 
vative coefficient when the amplifier is not connected in), 
R/ Rie is the ratio of the circuit's active impedance to the 
impedance of the leakage condenser, U;, (p) and Up. p) 
are the transforms of the input and output voltages respec- 
tively. 

In Fig. 1 there are constructed the relative amplitude 
and the phase frequency characteri tics. The upper bound- 
ary of the working range of the frequency characteristics is 
due to the presence in (1) of the aperiodic link (in accord- 
ance with the conditions adopted earlier for the phase fre- 
quency characteristic, (QT);,9, = 1). The lower boundary 
of the working range is due to the static term in the ampli- 


(1) 


*In differentiating slowly varying signals, an RC circuit 

is more convenient than an RL circuit since, by means of 
an RC circuit, one can easily obtain a larger time constant 
for a lower sensitivity to noise. 

t An RC differentiator is described, for example, in [2]. 
tTo increase the time constant one uses a condenser with 


a large capacity, in which the leakage impedance can be 
,commensurable with R. 
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Fig. 1. Phase frequency and relative ampli- 
tude frequency characteristics of a differen- 
tiating RC link. 
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fier’s transfer function. With this, the magnitude of 
(2T) min» a8 is clear from the phase frequency characteristics 
constructed for various R/R;., decreases proportionately to 
the decrease in R/ Rie. However, for each type of condenser 
there exists an absolute magnitude of 2 .;,, which depends 
only on the coefficient of impedance of the leakage con- 
denser kje = RieC. 

By substituting C = kj, / Rye in (1) we obtain 


R/R te'“te ? +- 1) 


ale 
RR kp +1 + RR 


Tout (yp) __ 
Vin (P) 





(2) 


le 


It is clear from (2) that, for low frequencies, the effect 
of the static term on the transfer function depends only on 
the magnitude of k,,. For example, for A ¢ = 20° and 
kj, = 3000 wf-meg, Qmin = 9.001 cps. 

It is clear from1) that the differentiating circuit's 
time constant coincides with the derivative coefficient. To 
increase this latter, we connect an amplifier to the differen- 
tiating circuit's output. 

If we assume that the amplifier has a linear static char- 
acteristic, then the derivative coefficient is increased by a 
factor of k , where k is the amplifier’s gain**. 

The size of k is limited by the increase, at the ampli- 
fier’s output, of parasitic signals due to internal noise, and 
also by the decrease in range of variation of the input signal 
corresponding to the linear portion of the amplifier’s static 
characteristic. 

The effect of internal noise on the working range of the 
frequency characteristic can be estimated by the ratio Yp 
of the parasitic signal u TT, as a fraction of the input signal, 
to the maximum value of the useful signal during one period 
(for Win t) = B sin Nt): 
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y= (3) 
4p 
where B is the amplitude of the input signal. 

On the other hand, for operation on the linear portion 
of the amplifiers static characteristic, it is necessary that 


the following inequality hold for u;_, , 
BQT < ig 1. (4) 


The aforementioned peculiarities of amplifier operation 
will not affect differentiator operation if, in the given fre- 
quency range, as B varies within given limits, y, remains 
much less than unity and satisfies(4). Figure 2 shows the 
graphs constructed for (QT) pin =f{B /uy) and (QT) ax : 
f (8/u,) for yp = 0.05 and ujp }/Ulp = 1000 (for k = 1000, 
It is clear from these graphs that, because of the limitation 
on the magnitude of (QT) min and the limits of variation 
of the input signal's amplitude, one can increase the deri- 
vative coefficient only by a factor of 1000. 


B. Differentiator in the Form of an Ampli- 





fier with Aperiodic Feedback 





This differentiating device consists of an amplifier, 
shunted by an aperiodic feedback path, and a subtraction 
block in which the input signal is subtracted from the am- 
plifier’s output signal (Fig. 3). The subtraction block is 
included so as to eliminate dc components from the am- 
plifier’s output signal. If the amplifier's characteristics 











Yout 
































Fig. 3. Block schematic of a differentiator 
in the form of an amplifier with aperiodic 
feedback. 1 is the amplifier, 2 is the aper- 
iodic feedback link,and 3 is the subtraction 
block. | 


**It is assumed that the amplifier’s input impedance is 
much greater than the RC circuit's output impedance. 

tt We understand by Up an averaged value of the parasitic 
signal, the magnitude of which depends on the amplifier 
circuit. 
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are assumed to be linear, then the transfer function of such 
a system will have the form 


Vout(p) _ 
Vin (p) 








(1 +- ki) k ( T 


i i 
ky P+ 1) Tp .a 1, (5) 


k + ky + 


~ (k-+-hi +k) M1 4 <q+! 


where ky = R/Rj, g and Rin g is the amplifiers input im- 
pedance. 

From (5) one may show that, to eliminate the dc com- 
ponent from Up. +, k and ky must satisfy the relationship 


ky =. 7 ° (6) 


By substituting ky from (6) in (5), we get 
Upyr (?) (k— 1)? » Tp y" 
Your (P) _ (K— FV 7 ( 1 
To) EH’) 
Figure 4 gives the amplitude and phase frequency 
characteristics of the system. It is clear from these graphs 
that, as k increases, the magnitude of ( QT) increases. 
However, in the given case, the derivative coefficient 
equals the time constant of the RC circuit. From the ex- 
pressions ky = R/ Rin a, ky = 1/(k-1), and T = RC fork » 1, 
we can cbtain 








R 
T= —i8¢. (8) 
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Fig. 4. Phase frequency and relative 
amplitude frequency characteristics 
of the differentiator in the form of 
an amplifier with aperiodic feedback. 


Thus, as k increases, T decreases. This decrease can 
be compensated by increasing the capacity of the conden- 
ser until Rj. becomes commensurable with R;,, ,. In prac- 
tice, even for an amplifier with R, , = 10 meg, for k = 100 
and C = 500 uf, T = 50 sec. Thus, in a linear system there 
is a limitation on the size of Tp. From whence,the mag- 
nitude of ( QT) max in practice equals 40. 

In a linear system no limitation on (QT) min “28 ob- 
tained, since the dc component due to the condenser leak- 
age was compensated here in the subtraction block. How- 
ever, in practice the subtraction block has a finite accuracy 
in implementing the operation of subtraction. We assume 








that the error with this comprises 0.1% of the dc compo- 


nent of u,** . As in the previous case, we define the 
ratio - 
0.001 
age thn (9) 






If, as before, we set y, = 0.05, we get that ( QT) min®™ 
=0.02. 

The upper boundary of the working range does not de- 
pend on the amplitude of the input signal if u in,}/ 8 =1, 
where this holds for any frequency***. 


C. Servosystem with Proportional Velocity 





Control as a Differentiating Device 





The simplified schematic of the servosystem with 
proportional velocity control is shown in Fig. 5a. The 
system consists of electronic amplifier 1, low-power mo- 
tor 2 controlled by the electronic amplifier, reducer 3, 
and feedback position transducer 4. In a first approxima- 
tion, we assume that all links of the system have linear 
characteristics, and that the time constants of the ampli- 
fier and motor may be neglected TTT. We write the 
equations of such a system in the form: 


Um= “ak (10) 
v=Kim (11) 
t= a ' (12) 

i Vr (Pp) 
fh (p) = 72 — (13) 
s= Win “fb (14) 


where € is the error sabing + is the voltage on the con- 
trol winding, v is the motor speed prior to the reducer, 

v, is the motor speed after the reducer, Vi p) is the trans- 
tens of Vie Up (p) is the transform of the feedback posi- 
tion voltage, Ujn is the input signal voltage, k, is the am- 
plifier gain, k, is the reduction factor, kp), is the position 
feedback gain, showing the change in ug, for one rotation. 

From (10)-{14) we obtain 


Vip) me a O 
kk mn TP +1 Usn (P), (15) 
é u o(u 
; rh had 
! 
| 





Fig. 5a. 





ttt is obvious that, for k = 100, the amplifier’s background 
noise will be less than the error from subtraction. 

*** This can be verified by defining the signal at the am- 
plifier’s input as the difference between u;,, = Bsin Ot and 
the feedback signal. 

ttt The permissibility of this last assumption will be ex- 
plained below. 
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where U (p) is the transform of us (ds V(p) is the trans- 
form of v(t) and T = k/ kakmkp- 

Expression (15) coincides in form with (1) for the case 
when R/R), = 0. For the measurement of slowly varying 
signals, as was established earlier, the value of T should 
be chosen large (for R/ Rie = 0). Therefore, for low -pow- 
er motors, the coefficient k,/kak,,kfp will be larger than 
the electrical and electromechanical time constants; the 
working range of frequency characteristics can be found 
from Fig. 1 (for R/Rye = 0,(2T)min = 9,(9T). 12, * 2). 

As the servosystem's output signal, one may use e€, 
Un’ and v (a voltage proportional to v). Measurement of 
€ is equivalent to the case of a passive RC network with an 
ideal capacitor, while measurement of up and v is 
equivalent to the case when an amplifier is connected at 
the output of the RC network, 

The dynamics of the actual servosystem differ from 
those of the linear system principally because of the non- 
linearities of the static motor characteristics and free play 
in the reducer. It was shown in [3] that the function 
v = F(u,,) for low-power asynchronous niotors has the form 





2 Mm “st _ 
v __ __“m. max. (16) 
? max + in — “ st . 
“mn. max 


where V jax is the maximum motor speed, Up max is 
the maximum voltage on the control winding, u,, is the 
voltage corresponding to starting. 

For Up K Upp max» (16) can be approximated by 
(11). As shown by experiments, (16) is valid as u,, varies 
within the limits of up, min = Um = Um.max: BY Um.min 
we understand a quantity somewhat larger than u,,. The 
starting voltage is determined by the starting torque, which 
depends on the load on the motor shaft. The load torque, 
in its turn, depends on the friction in the bearings, which 
is not a strictly constant quantity. Therefore, the motor 
shaft can turn continuously, starting with some velocity 
Vin Which we agree to call the minimum velocity of 
continuous rotation. 

If the voltage on the motor winding lies within the 
limits of u,, S Up = Up, min there then occurs a discon- 
tinuous mode of motor operation. For ujp = B sin QT, a 
definite portion of the period will be comprised of dis- 
continuous system operation. We may reflect this idio- 
syncracy of actual servosystem operation by the ratio ym 
(similar to the ratio Yp): 

u 
ym= aon (7) 
where Ujn min 18 the voltage at the amplifier input which 
corresponds to Vim, min’ (17) is analogous in form to (3) 
but, in practice, ujn min is significantly larger than Up: 

For input signals with shallow slopes, the free play in 
the reducer manifests itself in the servosystem. Because of 
the presence of the free play, as the sign of the slope of 
Ujn(t) changes,the speed of the motor will not correspond 
to the slope of the input signal during the time that the 


356 


free play is being taken up. It can be shown(cf. Appen- 
dix 1) that the ratio of the motor speed at the moment 
when the free play has been taken up to the maximum 
speed during the period (for Und = B sin t) equals 


o 


(mee)? 4 
Ye = 4.65 (3,2) ; ° (18) 


VOT 

where mg is the free play in the first stage of the reducer, 

The upper limit to the working range can be given ip 
the form of (4), where by u;,, , we shall understand u,, 
corresponding to the boundary of approximation of non- 
linear relationship (16) by the linear static characteristic 
of (11). 

LT 





10 
02 











Fig. 5b, 


Figure 5b shows the curves constructed for coTytit 
and (QT);y4x 48 functions of B/uj,, pin fOr Ym = yf = 9.05, 
Uin.1/ Yn min = 2% Key = 500uUjn min 2nd mg = 0.0001. 
It is clear from these graphs that, in the range 


20 = B/Uin. min = 50,(QT)min © 0-5 and (QT) A 4, *1. 


Discrete Differentiators 





A discrete differentiator measures the increment of the 
input signal during a finite interval of time: 


Ag(t+ At) _ g(t+ At)—g(t) 
At ie At 4 (19) 





where Ag(t +Art) is the increment of the input signal at 
time t+ At, g(t+At) is the input signal at time t +At, 
gt) is the input signal at time t, and At is the time inter- 
val between measurements of the input signal. 

A discrete differentiator may be constructed either on 
the principle of quantized input signal g(t) or on the prin- 
ciple of quantized time. In measuring the speed of change 
of a slowly varying input signal, Ag will be small, so that 
it is more convenient in practice to use the second prin- 
ciple (the value of At can be chosen arbitrarily). 

The dynamic characteristics of discrete differentiaton 
depend on the method of measuring the increments of the 
input signal. One can cite two basic methods of measur- 
ing Ag: 


ttt In correspondence with (17) and (18), two curves are 
given for(Q T),nin- For a given value of B/ujn min» the 
lower boundary of the frequency characteristic corresponds 
to the larger ordinate. 

‘It is understood that we are speaking of discrete differen 
tiators with quantized time. 
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1. Comparison of the instantaneous values of of t+T ) 





and g(") , WhererT is a fixed interval between two measure- 


ments{cf. (19)]. 

In (19) we replace the continuous argument t by the 
discrete argument nr (n corresponds to the moments at 
which Ag is measured). Then, letting g(t) =B sin Qt, we 
obtain 


Agi 1 = [sin Q(n +- 1)t — sin Qnt]. (20) 


By transforming (20), we find that 


Sates) = = sin = cos (Qne a ea) - (21) 

Expression (21) can be considered as the system's fre- 
quency characteristic for a discrete argument. This means 
physically that (21) shows the difference 4g,/Br from the 
ideal derivative only for the moments of time nr. The 
relative amplitude and phase frequency characteristics are 


given in Fig. 6(curve 1). As the argument of the frequency 


characteristics, we have here chosen Qr. It is clear from 
these graphs that the working range of relative values of 
Nr has just the upper boundary of (Q rT), © 0.6. 


2. Comparison of the results of integrating the input 
signal over the intervals of time [t,t+r ] and [t+T, t+ 2r], 


respectively : 








Aa P t+er t+t 
Be =|) 6ae— J g(a). (22) 


As in the previous case, it is assumed that g(t) varies 
sinusoidally. We obtain the frequency characteristic in 
the form 


2 , 3 
Ags (et 207) __ a. sin? 27 cos(@ + Qt), (23) 


where ¢ is the current value of the argument (analogous 
tonr in number 1 above). 
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Fig. 6. Phase frequency and relative ampli- 
tude frequency characteristics of discrete dif- 
ferentiators. 1 is for a discrete differentiator 
of the first type, 2 is for a discrete differentia - 
tor of the second type. 





The relative amplitude and phase frequency character- 
istics are given in Fig. 6 (curve 2)? It is clear from them 
that the upper boundary of the working range has been low- 
ered, and that(Qr)inyax * 0.45. 

It is clear from (21) and (23) that the derivative co- 
efficient (for small QT) equals, respectively, 1 and Tr. < 
In the first case, T q can be increased only by connecting 
an amplifier at the output of the discrete differentiator 
and, in the second case,at thecost of increasing r. How- 
ever, connecting an amplifier reduces the working range 
(cf. above), and increasing r decreases 0 may (Fig. 6). 
Moreover, implementing the necessary operations to de- 
termine increments on the basis of analog technology leads 
to large errors. For example, in the discrete differentiators 
described in [4-6], there will be large dc components at 
the output due to the nonideal nature of the elements. The 
transition to digital technology, as will be shown below, 
retains the working range of the frequency characteristics 
with an arbitrarily large derivative coefficient. 

We consider a discrete digital differentiator. To mea- 
sure the results of integrations over the interval r, it is 
convenient to present the input signal g(t) in the form of 
frequencies of trains of pulses (pulse-frequency modulation): 


1(t) =f, + lg (t), (24) 

where f, is the minimum frequency, corresponding to 
f(t) =0, and 2 is the slope of the analog-digital trans- 
formation‘. 

Figure 7 shows the block schematic of a discrete 
digital differentiator. The transformation corresponding 
to (24) is implemented in device 1. As a reversible coun- 
ter for the differentiation of slowly varying signals, one 
may use the low-power pulsed step-by-step motor 2. In 
switching phase, the rotor of such a motor is rotated by 
one and the same angle 4g. The use of a pulsed motor, 
in addition to simplifying the counter scheme, permits the 
implementation of a digital-analog transformation by 
means of inductive transducer 3. The intervals of addition 
and subtraction are given by synchronous motor 4. The con- 
trol scheme in 5 provides: a) switching of motor phase 
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Fig. 7. 























"In Fig. 6, the relative amplitude frequency characteris- 
tic was constructed in the dimensionless units Ag, / Br? 
instead of Ag, /Bro. 

One can convince oneself of this by expanding sin(Qr /2) 
in a series and retaining just the first term of the expansion. 
“such a transformation can be implemented either by 
means of a special device or by the use of transducers with 
digital outputs. 
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with frequency f(t); b) reversal of the direction of phase 
switching upon transition from addition to subtraction; c) 
transmission to the differentiator’s output of the measured 
increment of the input signal; d) reset of the system to 
its initial state. 

We now consider what the conditions are which must 
be satisfied by the digital differentiator's parameters in 
order to provide the necessary working range of the fre- 
quency characteristic (2,1, t0 2 .,,): 

The integration time rt can be found, for a given 
upper boundary of the working range of the frequency 
characteristic, from the phase frequency characteristic 
of an ideal discrete differentiator in which (2T) max 
equals 0.45 (Fig. 6, curve 2): 
0.45 
Q 


- 
. 


(25) 





max 

In a digital system, Ag will be expressed in the form 
of some number m of pulses. The general expression for 
m, as well as the expressions for Mpjjn 42Nd Mya for a 
given frequency and amplitude B in the case of a sinu- 
soidal input signal, are derived in Appendix 2. The quan- 
tity Min Characterizes the minimum accuracy in meas- 
uring 4g over one period of signal variation. By consid- 
ering Mmin to be a given quantity, we obtain from (A.10), 
after substituting the value of r from (25), 2 = Qmin and 
B= Buin the expression for 





8m Q 2 
Soe min ma 26 
Fa max (% : (26) 


Brin 
With this, the maximum number of pulses necessary 
to make the transformation from digital form to analog 
will equal, on the basis of (A.11), m max ¢2Bmax/ Bmin) x 


X Mmin (Q max/2 min)*®. 


Finally, by taking into account the optimal relation- 
ship between f and IB(cf-, Appendix 2), we find from 
(A.12) the expression for the maximum capacity of the 
counter: 

8Bmax (aes ) 


Nmax= pa "™nin\|\ pa 
Brin Onin 


Of all the parameters of the digital differentiator, it 
is the quantity Mp4, Which may be limited in practice 
since, with increasing Mpa x for a given type of digital 
differentiator,there is a decrease in the magnitude of Ag 
which, in practice, has a lower limit. Therefore, for given 
Mmin and Mmax: 4 digital differentiator has a bounded 
working range. For example, 2.4, / Qmin = 22 for 
Mmin = 3+ Max = 10,000 and Bmax/ Bmin = 5. 


In contradistinction to the second type of discrete 
differentiator which was considered earlier, in the given 
digital differentiator Ty depends, not only on Tr, but also 
on the proportionality factor between the angle of rotation 
of the transducer’s shaft and the change in voltage at the 
output of the induction transducer. Here, the proportional- 
ity factor does not depend on the working range, and Ty 
may in fact be made as large as desired. 


Comparison of Dynamic Characteristics of 





Differentiating Devices 





The results of a comparison of the three types of ana- 
log differentiating devices and the digital discrete dif- 
ferentiator are presented in the following table 




















Ampl, with Sefer. 
Ro network aperiodic feed - Servo sain | Digital differ 
plus ampl. |pback- entiator 
0.4 Not limited 
Qin i/sec 0,001 0.001 > 
1 1 
Qnax 4/sec ar 0.8 r 20 Q nin 
T/T 1000 100 1000 Not limited 
B “in. | 
max Te. “inl “in.d i?” ° 
ir 50 50 2.5 5 
min 














It is clear from the table that only the servosystem 
and discrete digital type differentiators have no limit on 
the smallness of the quantity 2,,;,. In the digital dif- 
ferentiator, moreover, the derivative coefficient can be 
made arbitrarily large. The sole superior feature of the 
analog differentiator is its greater simplicity vis-a-vis 
the digital system. However, for the differentiation of 
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slowly varying signals, one can use a pulsed motor in the 
digital differentiator, which significantly simplifies the 
entire system. 


Summary 


Analog differentiating devices with large working 
ranges (an amplifier at the output of an RC network, or 





jiven 


ana- 
™ 


n the 
the 


ng 


, or 





an amplifier with aperiodic feedback) have a limited mag- where T, is the length of the ith period. 


nitude of 24, and of the derivative coefficient. In an By integrating (A.6) we get 

electromagnetic differentiator, the quantity Qj, can be IB r 

made arbitrarily small, and 2,,,,/ Qmin can be made foT,— A 008 (QT; + ¢) + Dp os =i. 

less than in electronic differentiators, but the derivative 

coefficient is of the same order of magnitude. We now find the expression defining the number n of 


It is advantageous to construct discrete differentiators  2®Cumulated pulses during the interval [y, +r }: 
in the form of digital systems. A digital differentiator per- 1B 1B 
mits one to have an arbitrarily small value of Qin and foT 1 — ZF 008 (QT; + @) + Oo se=i, 
an arbitrarily large derivative coefficient. 1B 1B 
foT a — “Gy 008 (QT: + QT, + ¢) + ZF 008 (p + QT;) = 14, 


Appendix 1 (A.7) 
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For simplicity of computation, we assume that the 
motor characteristic v = F(u,,) is linear. The free play 
in the reducer manifests itself when the sign of the input 
signal's slope changes. We now determine the motor speed 1B 
at the moment when the free play has all been taken up, + Fy 008 (9 + QT, +...+O7,_,) = 1, 
given a sinusoidal input signal: 


IB 
foT , — D °° (Q7;+...+ QT, -+ @) 


02 7? Ip _1B . 1B 
vp = kh B (1 — cos Qte) = kk pt, (A.1) fo iT— _ (ean+ ®) + reegwem 


where vf is the motor speed at the moment when the free 1B 1B 

play has all been taken up. Knowing the shaft rotation fot — “qy 008(Qt + g) + Q coe =A. 
angle corresponding to taking up of the free play for a 
given reducer transfer number, we can, by integrating (A.1) 
and solving the resulting equation for te, determine this 
quantity in the form 


One may show analogously that, during the interval 
([¢+ T,¢+2r], there will, in the general case, be accumu- 
lated (n+ m) pulses: 


3 
/~ be lB lB 
tr=// TF Bo (A.2) fot — “Gy 008 (227 + ) + FH cos (gp + Qt) =n 4 m. (A.8) 
By substituting (2) in(1), we find the expression for By subtracting (A.7) from (A.8), we get 
V¢ in terms of ¢¢: 1B 
2 m =~ [2 cos (Qt + 9) — cos (22 + @) — cos 9]. (A.9) 
Pf 
a kak B-1,65 (=) (A.3) For ¢ = 32/2, the number m is a minimum: 
hell 41B — (A.10) 
We now find the ratio of the’maximum speed during Mnin= —gQ ~ sin? > sin Qt = 1Bt (Qr)*. . 
one period at a given frequency to the velocity v,: 
¢ For g=™, m is maximum. The approximate maxi- 
; Q\3 mum value of m is 
4 = “ar (Rip) ins (A.11) 
a hn m max ~ /Bt (Qr). 
The angle of rotation for free play take-up can be pre- 
sented as a product m¢k,, where mg Corresponds to the free Correspondingly, (A.7) attains a maximum forg = 1/2. 


play in the first stage of the reducer (the free play in the The approximate maximum number of pulses accumulated 
other stages has a lower value). By multiplying the num-  4Uring one interval of integration may be estimated by the 





erator and denominator in (4) by kg, we obtain finally expression 
2 *max~ * Vo + Bryay!). (A.12) 
mekns* 4 Th ill be 
4 = 1.65 ( ) =z (A.5) e counter w used optimally if, for a given n, 
Vor the quantity mppjp is maximal. It can be shown that the 


ratio My)j,/0, found from (A.10) and (A.7), will be maxi- 


Appeadiz 3 mal asl p/f,y—1. 





The length of the period between two neighboring 
pulses for a sinusoidal input signal with arbitrary initial LITERATURE CITED 
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CHOICE OF FREQUENCY RANGE FOR INDUSTRIAL PULSE- 
FREQUENCY TELEMETRY SYSTEM DEVICES 


A. M. Pshenichnikov 


(Moscow) 


Translated from Avtomatika i Telemekhanika, Vol. 21, No. 4, pp. 525-529, April, 1960 


Original article submitted November 17, 1959 


The paper considers the effect of frequency range of pulse-frequency telemetry system devices on the errors 
and noise resistance of transmission, speed of action, and receiver dynamics. 


The following factors influence the choice of frequency 
range of pulse-frequency telemetry system devices: 

1) telemetry error as a function of frequency range; 

2) noise resistance of signal transmission along the 
communications channel, and the effect of noise from the 
50 cps power supply; 

3) indicating instrument dynamics and speed of sys- 
tem action; 

4) bandwidth of the secondary channel. 


1. Telemetry Error as a Function of Fre- 
quency Range 








As a function of the ratio of the frequencies which 
correspond to zero value and to the nominal value of the 
measured parameter, the error « in transforming this param- 
eter to a proportional frequency in the transmitter gives 
rise to the telemetry error [1] 

a 
Se = TT, Tally’ (1) 
where f and fn are the frequencies corresponding to the 
zero value and to the nominal value,respectively,of the 
parameter being measured. 
Analogously, the error 8 arising in the receiver when 


the frequency is transformed to a proportional current gives 
tise to the error 


8 
Pee i— holy ( ) 
The ratio of the errors is 
Fe (3) 
a re) | —hell 


The figure shows the relationship of the ratio of errors, 
Me /a or8y_/B, to f/f. As is obvious from (3) and 
from the graphs on the figure, the telemetering error de- 
creases as the ratio fy / fy decreases. 


2. Noise Resistance of Signal Transmission 





along the Communications Channel and the 








Effect of Noise from the 50 cps Power Supply 





With weak fluctuation noise, the error is determined by 
the distortion of the pulse edges of the transmitted signals. 
This error, as a rule, is a very small quantity. 
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Dependence of the ratio of the 
telemetering errors to the ratio 
of the frequencies corresponding 
to the zero value and to the nom- 
inal value of the parameter being 
measured. 


An analysis of the noise resistance of telemetry signal 
transmission along a channel with weak fluctuation noise, 
based on the theory of potential noise resistance, was made 
by V. A. Kashirin and G. A. Shastova [2]. However, in con- 
sidering pulse-frequency modulation with a constant mark- 
space ratio PFM Il, the authors of [2] analyzec the trans- 
mission of triangular pulses with amplitude 2U,,,. The trans- 
mission of triangular pulses negates all the advantages of 
PFM II, and does not correspond to actual conditions. 

The noise resistance of PFM II to fluctuation noise can 
be determined by expanding rectangular pulses in a Fourier 
series. The pulse frequency in PFM II can be presented in 
the following manner as a function of the measured param- 
eter A: 


10) = fit Afr, (4) 


where fj is the mean pulse frequency and Af /2 = 
(fn-h)/2 is the frequency deviation. 


361 











The expression for a rectangular signal of amplitude 
2U,,, after expansion in a Fourier series, can be written 
in the form 

, 4U / 

A(h, t) = Uy + —"|sin 2 (7, -+ ca \l + 
7 3 (5) 
+ + sin 6x (/ <a ] 

The mean square error in transmitting the parameter 

Aq due to fluctuation noise is defined by the expression [2] 


S . Ss 
ome. s4/) = 7 (SSS 9 


/T 
j = ‘ «© [ @A(A, t) 
2V2Y S[SP pa 
0 * di. lai 
where o is the specific intensity of the noise in Wcps) 
and T is the time to establish the indication. 


If the filter passes the first and third harmonics then, 
after some transformation, we get, for T >» 21/fj, 








(6) 


Ps 


UAT? 





ome. sq PFM — 


Comparisons of the noise resistance with frequency 
modulation (FM) and with PFM II must necessarily be car- 
ried out for signals of identical amplitudes. With this, 
the amplitude of the first harmonic must be considered 
equal, not to 4U,/ m but ~ 1.15 U,,. In this case, 

0.22 ¢ (8) 
U AT" * 


For the same amplitude, the error with FM is (cf., 


“me. sq PFM= 


(2) 0.28 6 


eee g 
U ,,AfT” % 


eae. sq FM~ 


Consequently, if the use of filters provides transmission 
of the first and third harmonics, with PFM II, then with this 
the mean square error from fluctuation noise will be less 
than with FM with filters of the same passband Low-fre- 
quency pulse-frequency systems for transmitting first and 
third harmonics require the same (and sometimes a smaller) 
frequency range than the existing frequency telemetry sys- 
tems. 

N. V. Pozin [3] carried out a determination of the error 
with PFM due to the relatively strong noise when the basic 
signal distortions are due to the breakdown into pulses and 
pauses, formulas being derived for approximate determina - 
tions of the mean,and the mean square errors as functions 
of the probability of a unit distortion, the filter passband 
2F pb» the frequencies corresponding to the zero and nomi- 
nal values of the measured parameter, and the establish- 
ment time of the indicating instrument: 





ime. = V/ Tos [ (Pho) —(v— B)'] 


11) 
where P is the probability of unit distortion. ¢ 
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The minus sign between the parenthetical expressions 
in (10) and (11) is due to the fact that, with PFM II, the 
subcarrier frequency filter must pass at least the third har- 
monic and, consequently, Fob =3 fy: This latter deter- 
mines the choice of the frequency range for PFM II. Other- 
wise, the probability of unit distortion would not be the 
same for all pulse durations, as was assumed [3] in the 
derivation of (10) and (11). 

The table gives the results of computing 6. and 
5 me.sq for the telemetry devices of the pulse-frequency 
devices of Ts_EM Mosenergo, IAT AN SSSR and TsNIIKA 
for operation with filters of passbands of 40 cps (developed 
by IAT AN SSSR), 80 and 140 cps (telegraph filters), start- 
ing from the assumption that the formulas for determining 
5 me 40d 5 me.sq to a first approximation, reproduce the 
state of affairs when the transition is made from PFM II 
to FM [3]. 

In the table, Bo» 80s and Py denote, respectively, the 
probability of unit signal distortion for filter passbands of 
40, 80, and 140 cps, with identical signal voltages in all 


cases. 
It is clear from the tabulated data that, for the stan- 


dard telegraph filters, the very best frequency range of those 
considered is the range from 4 to 20 cps. 

In addition to the considerations already mentioned, 
the noise from the 50 cps power supply (line) has an effect 
on the choice of frequency range for the devices of fre- 
quency and pulse-frequency telemetry systems. The pre- 
sence of the frequencies of 25, 50, and 100 cps in the fre- 
quency range leads to a strengthening of the effect of noise 
from the power supply on the operation of the receiving de- 
vice in the neighborhood of the points cited. In the trans- 
mission of the telemetered quantity at a frequency of 25 
cps there appears, at the receiving side after some nonlinear 
transformations, a component of the first harmonic of the 
industrial frequency (50 cps) which is essentially influenced 
by the power supply. 

Since it is desirable to have the ratio fy/ fy as low as 
possible, so as to provide accuracy of transmission, one 
should choose, for the devices of pulse-frequency telemetry 
systems, a nominal frequency less than 25 cps. If one takes 
into account the possible oscillations of the line frequency 
(under trouble conditions ,the frequency in energy systems 
falls to 46 cps) and the fact that, in using standard telegraph 
filters with 140 cps pass bands, it is desirable to transmit the 
third harmonic, then one should take, as the very best nomi- 
nal frequency, one approximately equal to 20 cps. 

3. Dynamics of the Indicating Instrument 
and the’System's Speed of Action 








The least (zero) frequency of teletransmission is de- 
termined by the highest frequency of variation of the meas- 
ured parameter, Fin’ and, for analog reception, by the in- 
dicating instrument dynamics. Since industrial telemetry 
devices are used, basically, for the measurement of rel- 
atively slowly flowing processes, the condition 


fo > 2F m, (12) 
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flowing from Kotel'nikov's theorem, is always met in prac- 
tice. Thus, the zero (null) frequency is determined by the 
indicating instrument's dynamics. 

The practical behavior of a pointer indicating elec- 
tromagnetic galvanometer, measuring mean current of a 
condenser (or transformer) receiver, becomes completely 
satisfactory for frequencies at the receiver input of 4 to 5 
cps (current pulses of frequencies from 8 to 10 cps flow 
through the instrument). 

The speed of action of a telemetry device is defined 
as both the time to establish the frequency at the trans- 
mitter’s output and as the time necessary to transform the 
frequency to a new indication of the indicating device 
(for analog reproduction) or, again, as the time necessary 
to translate the frequency to code (for digital reproduc- 
tion). In the very fastest transmitter, the maximum time 
for establishing the output frequency will be when the sig- 
nal changes from its nominal valué to its null value. This 
time is defined by the rr 

Mh< <7: 

The speed of a telemetering system with analog re- 
ception is defined as the time to establish the indication 
of the indicating instrument and, for a frequency range of 
4 to 20 cps, is about 4 seconds. 

For digital reproduction, the time necessary to trans- 
late the frequencies into code is determined by the receiver 
error 8+. and by the frequencies corresponding to the zero 
value and the nominal value of the measured parameter. 

With simple pulse counting, 


7 
q = 
2 


(1 3) 


100 
Bie” (hy — fo) 





(14) 


With a doubled frequency and a transformation of 
PFM II to PFM I (pulse-frequency modulation with con- 
stant pulse duration), 

T= 50 


Pre % (iN fo) . 





(15) 


4. Secondary Channel Bandwidth 





As was mentioned earlier, when signals are transmitted 
by PFM II, the subcarrier frequency filter must pass at least 
the third harmonic. However, it is possible to use trans- 
mitting and receiving devices of a pulse-frequency tele- 
metry system with narrow-band filter, passing only the first 


harmonic of the nominal frequency. With this, teletrans- 
mission in a large part of the frequency range will be im- 
plemented, not by PFM II, but by FM. In this case, thanks 
to the narrow-band filters, the signal-to-noise ratio at the 
receiver input will be reduced. 

In the devices with secondary compression used today, 
the narrowest secondary channel is a subchannel of an 80 
cps width. Since the receiving device of pulse-frequency 
telemetry provides normal operation when a sinusoidal sig- 
nal whose frequency is of the order of 20 cps is applied to 
its input, telemetry devices with such stations will 
operate reliably. 


Summary 





1. Decreasing the {/ fy ratio lowers the require- 
ments placed on the accuracy of transforming de (voltage) 
to frequency and on the inverse transformation of frequency 
to current. 

2. Increasing the frequency range, fy to f , for the 
same filter passband , first decreases the mean square error 
for weak and for strong noise and, second, decreases the 
conversion (coding) time when digital reproduction is used. 

8. The highest teletransmission frequency for PFM II 
is determined by the noise from the 50 cps line end by the 
existing standard telegraph filters. For filters with a 140 
cps passband, it is desirable to transmit the signal's third 
harmonic. The highest frequency is about 20 cps. 

4. The lowest (null) frequency with analog reproduc- 
tion is determined by the indicating instrument's dynamics, 
and is about 4 to 5 cps. 

5. For a nominal teletransmission frequency of the 
order of 20 cps, a telemetering device can operate with 
the use of any telegraphic filter. 
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ON THE RELIABILITY OF SCHEMES FOR CONNECTING 
DISPERSED OBJECTS ' of 3 
, que 
V. A. Il'in cir 
(Moscow) 
Translated from Avtomatika i Telemekhanika, Vol. 21, No. 4, pp. 530-532, April, 1960 
Original article submitted July 27, 1959 
The paper considers the reliability of schemes for the parallel and series connection of a remote control of 
system's executive points in a communications line. The domain of applicability for each of these schemes is cir 
shown. tw 


In [1], the reliability of various remote control line 


structures for dispersed objects was considered. As a de- oe 

velopment of this work, we consider here the reliability of LD} jp} [DJ Tre iD a wi 
‘ : ‘ee 

schemes for connecting the executive points in communi EP | EP | EP | 


cations lines. The choice of some connection scheme or 
another determines, to a significant degree, the structure, Fig. 1 ar 
and significantly affects the structural reliability, of remote in 
control systems with dispersed objects (oil and gas industries, pr 
pipelines, irrigation systems, transportation, mines and 
factories, cities’ communal economies, etc.). To a cer- 
tain extent, this also applies to automated systems. 

Figures 1, 2, and 3 show schemes for parallel, series, 
and mixed connections of a group of m dispersed executive 
points(EP) to a common communications bus via distribu- 


tive devices D which might be parts of the EP's. Ina 
scheme of the mixed type, the individual distributive de- 
vices implement, for example, only the functions of 
switching devices (Fig. 3). Analogous schemes are also 
used for radio channels. 

The executive points on Fig. 1 are connected in paral- 
lel via the distributive devices D to a common communi- 
cations line which joins them to the central point (CP). 

In many cases, the distributive devices are executed in the 
form of matching transformers. The parallel connection 
scheme (Fig. 1) is employed in remote control systems with 
either frequency or temporal separation of the signals. 

With series connections (Fig. 2), the communications 
line is divided into a chain of segments which are connec- 
ted in series via the distributive devices D. With a series connection, each distributive device can 

In one of the simplest remote control systems with also implement the function of signal amplification, simi- 
dispersed switching [2], each distributive device contains lar to what takes place in radio relay lines with the same 
relay Contacts, connected in series in the line, and dis- connection scheme. 
connecting the remaining parts of the system in one of the We now consider the reliability of systems with paral- 
relay positions. When a series of selection pulses is trans- lel and with series connections of the EP'sin the communi- 
mitted, the “next® link is connected to the line at each cations line. Such connection schemes are the fundamen- 
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Fig. 3 
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cycle, and the line is “selected” successively in time. tal ones. 
In remote control systems with dispersed objects, if We make the following assumptions. 
these systems operate in accordance with the scheme of 1. All distributive devices in one and the same sys- 


Fig. 2 and are constructed of contactless elements with tem (scheme) are identical. For each of them there exists 
rectangular hysteresis loops, the distributive devices are a probability of a short-circuited input or output, P,., and 
executed in the form of transformers. a probability Pope of an opened circuit, for example, from 
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a winding blowout. The probabilities P,, and P for a 
given EP are mutually independent, and do not depend on 
the probabilities of short circuits or open circuits of the 
other EP's. 

The probabilities Ps. and P., of a nonoperative state 
of a given distributive device are related to the mean fre- 
quencies of failure due to short circuits fg and due to open 
circuits fone by the following expressions: 


“| 


* ~ P 
ke = Sand onc pre. 
sc opc 
Here, Tu and 7. ¢ are the mean nonoperating times 
of a given distributive device due, respectively, to short 
circuits and to open circuits during the period of time be- 
tween maintenance cycles; 


n n 
~_ _— 1 
T sc= — lim Tie, 1. = —!) T 
sc Py tse» “O56 lim a i ope 


where n is the number of faults (failures). 

2. The communications line itself is considered to be 
absolutely reliable. Such an assumption simplifies the 
analysis and can not change the character of the results 
in any essential way since, when account is taken of the 
probabilities of short circuits and open circuits in identical 
links of the line, there is only an increase in the corres- 
ponding P._ and P,,..- 

Then, for schemes with parallel and with series con- 
nections, the probability of system failure due to short 
circuiting of the distributive devices to ground or to the 
second conductor in a two-conductor line, without taking 
account of failures due to open circuits, will be expressed 


by the formula Pulte 


For mP,,. « 1, Da mPo. 


The probability of a nonoperative state of an object 
in a system due to open circuits in the distributive devices 


when parallel connections are employed (Fig. 1) isP =P... 


since the open circuiting of any distributive device disrup 
the operation only of the EP connected to it. 

For a scheme with series connections (Fig. 2), the 
mean probability of a nonoperative state of an object in the 


system due to open circuits in the distributive devices, with- 


out account being taken of faults due to short circuits, will 
be 


M3 


| 


pe) 
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where nj is the number of objects connected between the 


CP and the ith distributive device ( n; = i) and, for 


mP ope « 1, will be 


m-+- 1 
P =— 2 Popo 
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If an open circuit in a distributive device disrupts only 
the operation of the EP connected to it in a paralel-con- 
nected scheme, the same open circuit in a series-connected 
scheme will cause failure of an average of (m+1)/2 EP's, 

It also follows from the expressions obtained that the 
probabilities of an object's nonoperative state due to a 
short circuit in a distributive device in systems with paral- 
lel and with series connections are equal to each other 
when the probabilities of short-circuited distributive de- 
vices are equal. However, in a parallel-connected scheme, 
the probability of a nonoperative state of the system due to 
a short circuit in the distributive devices or in the EP's can 
be decreased to an arbitrarily small quantity by connecting 
the EP's to the line via cutouts or limiters simila:ly to what 
is done in connecting electric energy users to a power line, 
or radio relay points to a radio relay network. 

For frequency systems of remote control, and for sys- 
tems which use exponential transformers [3] for the time 
separation of signals, the EP apparatus is comparatively 
simple in parallel-connected schemes. 

If existing systems of remote control with time separa- 
tion of signals are connected in accordance with the scheme 
of Fig. 1, it is necessary to have, at each EP, a one-termi- 
nal set (incomplete in many cases) which is frequently un- 
reasonable (for a small number of signals on one EP). 

The series-connected scheme can be realized in re- 
mote control systems with time separation of signals with 
contact distributed switches [2] or in contactless systems 
with elements with rectangular hysteresis loops. 

Such systems, even when constructed from contact- 
less elements, have significantly less structural reliability, 
as our analysis shows. 

Frequency systems of remote control can also be con- 
structed in accordance with the series-connected scheme 
(Fig. 2) but this is unreasonable, in connection with the 
lower structural reliability of such systems. 

The use of series-connected schemes, from the point 
of view of reliability, may be justified principally in 
those cases when it is necessary to amplify signals at given 
points of the line, or it is necessary to switch groups of 
EP’s in communications lines with tree structures, if the 
solution of such problems gives rise to difficulties, when 
parallel-connected schemes are used. 
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( Moscow) 


Translated from Avtomatika i Telemekhanika, Vol. 21, No. 4, pp. 533-541, April, 1960 


Original article submitted February 12, 1958 


Methods are considered for increasing the reliability of operation of relay-contact devices by means of spare 
elements. A new method is suggested for sequential circuits which allows the device's reliability to be increased 


with the minimum number of spare relays. 


In recent years the questions of reliability of elec- 
trical devices have received great attention. With this, 
the problem of increasing reliability is solved by two 
routes. The first of these is the increase of strength and 
durability of the elements entering into devices, the re- 
placement of contact elements by contactless ones, for 
example, by semiconducting and magnetic elements, etc. 
The second route, independent of the first, is the wide 
use of stand-by elements, assemblies and blocks, i.e., 
the use of back-up. 

In these cases, the operating parts of the apparatus 
or of its individual blocks, or even individual elements, 
are supplied with reserves which are either switched into 
operation automatically when the basic set fails or operate 
constantly in parallel with it. Such reserve elements are 
“redundant,” since their introduction adds no new func- 
tions to the device. 

As an example of the use of the back-up method of 
increasing reliability, one may adduce the automatic 
track blocking circuit of the Moscow subway in which, as 
is well known, a type DSR relay is connected in parallel 
for every two tracks. 

In the engineering and scientific literature devoted to 
the questions of stand-by capability, one considers princi- 
pally the reliability of radioelectronic apparatus, con- 
structed from contactless elements and used in a number 
of special domains, for example, in computing technology, 
in radio relay communications, for the control of moving 
objects, etc. In correspondence with the purpose of such 
apparatus, there arises the requirement that it operate 
without failures during some interval of time, so that here 
reliability is defined as the probability that, in the course 
of a definite time interval, the device will correctly carry 
out its function. 

In industrial devices for automation and remote con- 
trol, one ordinarily considers a device to be sufficiently 
reliable if a random failure of its elements does not en- 
tail an incorrect functioning which leads to an unsafe state 
of the object. Operational failures in these cases are con- 
sidered as “protective”, and are ordinarily accompanied 
by either a blocking device and an alarm indication of 
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the fault, or a transfer of the device to a safer state. Thus, 
for example, when a railway element fails in the automatic 
blocking circuit, a red light signal automatically appears, 

requiring that the train stop. 

However, as is understandable, each such case of pro- 
tective failure disrupts the continuity of the technological 
process. With the modern tendency to intensify production, 
it becomes important to increase the reliability and unceas- 
ing operation of industrial automation and remote control 
devices. 

In the present paper we consider several questions in 
the construction of reliable relay-contact devices which 
permit a single failure to occur without disruption of opera- 
tion. The principles to be considered in the sequel may be 
extended also to contactless devices using relay -acting 
(switching) elements. 

The effect achieved when the ordinary method of 
parallel back-up [1] is used can be illustrated by the exam- 
ple of the simplest relay-contact circuit (Fig. 1a), whose 
function is to close circuit a-b when relay Al is excited. 

We denote the probability of correct closure of the 
relay contact by p(A) = p, and the probability of the event 
that circuit a-b will be closed by p(a-b). In the given 
case, p(a-b) = p. 

Connecting relay A2 in parallel to relay Al (simple 
back-up) and connecting the contacts of relay Al and A2 
in parallel (Fig. 1,b) increases the probability of correct 
closure of circuit a-b up to the value p(a-b) = 2p-p*, ie., 
for 1-p « 1, increases the reliability by almost one order. 

Conversely, a series connection of the contacts of re- 
lays Al and A2 (Fig. 1c) lowers the probability of closure 
of circuit a-b to the value p(a-b) = p®. 

Figure 1 d-f, also gives several variants of the sim- 
plest schemes for the two- and three-fold back-up of re- 
lay Al, and the corresponding values of the probability of 
closure of circuit a-b. To determine the expressions for 
the probabilities, it was assumed that all the relays A1,A2, 
A3, and A4 were monotypic, as were their contacts, that 
the probability of contact closure was the same for all, and 
that cases of incorrect operation were statistically inde- 
pendent. 
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As is obvious from these examples, by paralleling the 
corresponding number of relays (with multiple back-up) 
and by properly connecting their contacts in a circuit which 
is equivalent in its operation to the action of one contact 
(or any other more complicated structure), one can attain 
any previously given degree of reliability, even with the 
use of elements which are themselves of low reliability. 

Of course, such back-up leads to a significant increase in 
complexity and cost of the devices, and to a large expendi- 
ture of circuit elements. 

For sequential relay-contact circuits, a more economi- 
cal back-up, which may be called “combined* is possible. 
The essence of this back-up amounts to this, that the back- 
up elements, the number of which, as a rule, is less than 
the number of the basic elements, back each of the basic 
elements up at that moment (cycle) when the operation of 
the entire device is determined by the operation of this 
basic element. 

We consider this new back-up method in terms of the 
example of a centralized control relay distributor of type 
DVK-3 — DVK-3a. The functional schematic of such a 
distributor is shown in Fig. 2, and the sequence of its ac- 


tions is given in Table 1 [2]. One easily convinces one- 
self that a fault in any one of the counter relays 1-10, or 
of their contacts, causes the distributor to cease operation, 
as well as a fault in the action of the entire one-terminal 
set (and, if this one-terminal set is the central control 
point, the entire system fails). 

We isolate, in Table 1, those cycles which would be 
terminal points, completing each distributor step, i.e., 
cycles 0,3,6,9,12,15, etc., and we write the state of the 
distributor at each of these cycles in the form of a binary 
number. We adopt the convention that an excited relay 
state will be denoted by a "1", an unexcited state by a 
"0". Table 2 gives the sequence of these binary numbers, 
which completely characterizes the sequence of operations 
of the distributor relays and which mirrors, by the correspon- 
ding binary numbers, their states at each step. Thus, for 
example, the fifteenth cycle of Table 1 corresponds to the 
binary number 000001000010 (the fifth cycle of Table 2), 
while the eighteenth corresponds to 000010000000. 

By considering this table of binary numbers as a set 
defined by given combinations of relay states, and consid- 
ering all the rest of the set of binary numbers which are not 


TABLE 1 Sequence of Operations of a Relay Distributor L is the line relay, RR is the 


receiving relay. 
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TABLE 2 Sequence of Distributor Operations, Given in 
the Form of a Table of Binary Numbers 
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TABLE 3 Distances between Relay State Combinations 
at Each Cycle of Distributor Operation 
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written in the table as distortions of these combinations, one ber of nonrepeating combinations with this will equal 2” 


can try to use, for constructing self-correcting circuits, the 
same methods which were developed for constructing self- 
correcting codes, also described by binary numbers which 
are tabulated in analogous fashion [3,4]. 

It is well known that, for the construction of error-free 
codes, it is necessary to create redundant elements of the 
combinations, and thereby to achieve the result that the 
code combinations differ from each other by a sufficient 
number of bits ("distance” between combinations). Thus, 
to construct a code for the correction of unitary errors, it 
is necessary that the code combinations differ from each 
other by at least three bits. 

We now extend these propositions to the given com- 
bination of relay states. Table 3 gives the values of the 
“distances” between the combinations of relay states (the 
binary numbers of Table 2) for the distributor of the cir- 
cuit in Fig. 2, from which it is clear that, in a number of 
cases, this condition is not satisfied, so that it is required 
to introduce auxiliary elements in the combinations (re- 
lays) to provide the requisite redundancy. 

The number of auxiliary relays can be determined 
from the well-known expression of Hamming, 


a 
Pre 
SS n + 1 ; 





where np is the number of basic relays in the circuit and n 
is the total number of relays, including the auxiliary ones. 

The value of m = n-ng, i.e., the number of auxiliary 
relays, can be found from the following table, set up in 
accordance with this formula: 


n=12345 6 7 8 9 10 44 12 
n=35679 10 11 12 13 14 15 16 
m=:23334 44444444 


In the case under consideration, it turns out to be suf- 
ficient to add four auxiliary relays. Connecting these re- 
lays in the over-all operation of the distributor can be im- 
plemented in various ways with the same condition that re- 


peated combinations not occur. As is well known, the num- 


i.e., the number of arrangements which differ from each 
other in at least one place. It is most convenient to con- 
struct the table of auxiliary relay operations such that, at 
each cycle, the minimum number of auxiliary relays, i.e., 
one relay, will change state. This condition is satisfied 
by, for example, the Gray code [5], used in a somewhat 
different form for constructing Table 4, giving the opera- 
tion of auxiliary relays I-IV. One easily convinces one- 
self that, with this, the distance between the binary num- 
bers of the combinations of relay states at different cycles 
of distributor operation in no case is less than three, which 
is the requirement imposed. 

We now consider three cases of constructing self -cor- 
recting circuits for relay counters. 

1, Circuits for self-correction it there appears, in a 
binary number characterizing an operate cycle of a given 
relay, a zero instead of a one (for example, if a prior 
counter relay did not close its contact in the same operate 
circuit of the given relay). 

2. The same, but tor the appearance of a one in the 
binary number instead of a zero (for example, for the seal- 
ing of a contact of the prior relay). 

3. The same, but for the appearance in a binary num- 
ber of a zero instead of a one, or a one instead of a zero. 

In the first case, to provide for operation of the given 
next counter relay, its operate circuit, passing through 
normally open contacts of the previous counter relay, must 
be backed up by the contacts of that auxiliary relay which 
also changed its state at the previous cycle (a normally 
open contact if it were excited, and a normally closed con- 
tact if it were unblocked). 

In the second case, protection from a false operation on 
an incorrect cycle due to a random closing of a contact of 
a previous relay is attained by introducing into the operate 
circuit of the given relay the contacts of three auxiliary 
(back-up) relays which do not change their. state at the 
given cycle, and the combination of whose states at the 
given cycle differs by the state of at least one relay from 
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TABLE 4 Binary Numbers for Operation of the Distributor with Auxiliary Relays 
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Fig. 3 


the combinations of all other cycles. These combinations 
are given in the right part of Table 4. 


The third case comprehends the first two. To con- 

struct a protected self-correcting relay counter circuit, 
it is necessary to create as many back-up operate cir- 
cuits as there are series contacts of basic and auxiliary 
relays in the basic operate circuits. Certainly, the con- 
tacts of these back-up circuits can be joined by the use of 
diode elements, as shown in Fig. 3. 

It is necessary to dwell on the cut-off relay circuits. 
As is clear from the scheme of Fig. 2, the blocking and 
cut-off circuits of the counter relays of an ordinary relay 
distributor are interconnected, and pass in series through 
inherently normally open contacts and normally closed 
contacts of the relay of the next higher ordinal number. 
For the case of nonclosure of this normally closed contact, 
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it might be backed up by contacts of auxiliary relays. How- 
ever, with nonrelease of one of these contacts, the relay is 
not released, which disrupts the distributor's operation. 
Therefore, correct distributor operation can be achieved 
only by separating the relays’ locking and cut-off circuits, 
whereby the release of the relays is implemented, not by 
means of the supply circuit, but by counter flows in the 
relays’ second windings (active unlocking). With this 
method, unlocking of a given counter relay is obtained by 
means of the operate circuit of the next relay of the same 
parity (i.e., even or odd, respectively), and may therefore 
be identified with it, the differentiation being only by the 
state of the forcing element, and the cycles on which they 
are used. 

Figure 3 shows the circuit for connecting several of the 
basic counter relays of a self-correcting distributor, con- 
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Fig. 4. a is the operate circuit of relay II, b is the release circuit of relay II; 
A is the operate circuit of relay 2, B is the operate circuit of relay 8, C is the 
operate circuit of relays 6 and IV, D is the release circuit of relay 9. 


structed in accordance with Table 4 and with the Condi- Using the same principle, one may construct multiple 
tions to be given now. relay contact circuits in which are corrected double, tri- 

As was stated earlier, auxiliary relays I-IV are connected ple, etc. errors, i.e., which will provide correct operation 
analogously to the basic counter relays in that, at the for simultaneous failures of two, three, etc. relays, at the 


corresponding cycles of distributor operation, they operate expense of a corresponding increase in redundancy. 
or are released simultaneously with the basic relays. Since 


auxilfary relays I and II operate and release twice during LITERATURE CITED 
an operate cycle of the distributor, one then uses for them 
the binary number of such circuits, as is shown in Fig. 4 
for auxiliary relay II. 

The difference in the circuits for the basic and the 
auxiliary relays stems from the fact that there is used, in 
one of the back-up circuits for operating a basic relay, a 
proper contact of that auxiliary relay which must also 
operate (or release) on the same cycle. To obtain coher- 
ent operation both of the basic and of the auxiliary relay - 
by this back-up circuit, this circuit is so separated forthe °% A- A. Kharkevich, An Outline of General Communi- 


1. EE. F. Moore, and C. E. Shannon, “Reliable circuits 
using less reliable relays,” J. Franklin Inst. 261-262, 
No. 2-3 (1956). Faac. 

2. F.I. Marushko, “The DVK-3a and DVK-3 central- 
ized control systems,” Collected Papers on the New 
Signal Centralization and Blocking System and 
Communications Technology [in Russian KTranszhel - 
dorizdat, 1955), 


basic and the auxiliary relay that the winding of the auxil- cation Theory [in Russian] (Gostekhizdat, 1955). 

iary relay is supplied only after the basic relay has operated. 4: %- W. Hamming, “Error detecting and error correct- 
Such a relay distributor provides correct operation for ing codes," Collected Papers (subject not given) 

any single fault, i.e., for the failure of any one element, edited by A. M. Petrov [Russian translation) (Tl, 1956). 

either basic or auxiliary. The distributor's output circuit 5. Fundamentals of Automatic Control.2 [in Russian) 

must, of course, be constructed like an operate circuit, (edited by V. V. Solodovnikov) (Mashgiz, 1959), 

i.e., it must be backed up. Chap. 2. 


Thus, increasing by almost an order the reliability of 
distributor operation is attained most economically by adding « his Jast sentence is factually(and mathematically ) 
four relays only, instead of the ten necessary for parallel incorrect, as are the data of Table 3 on page 368, 
back-up. It should be mentioned that, according to Ham- The maximum number of basic relays for which four auxiliary 
ming'’s formula, four auxiliary relays can provide back-up _ relays will suffice is not 15, but 11 [Publisher's note]. 
for 15 basic relays. * 
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RECTIFIED VOLTAGE REGULATOR 
O. G. Malkina 
(Moscow) 


Translated from Avtomatika i Telemekhanika, Vol. 21, No. 4, pp. 542-547, April, 1960 
Original article submitted July 21, 1959 


A new variation of a circuit for regulating low rectified-voltage feeding loads on the order of hundreds of 
watts is discussed, As an example, curves and oscillograms of the output voltage change as a function of the 
load change are presented for a regulator wired according to the given circuit. 


The use of transistors for computing machines and reactors Ch,-Chg, to the three-phase transformer. The 
other devices requires power sources of low dc voltages. stepped-down ac voltage from the secondary windings of 
These voltages may be obtained by lowering and rectify- the transformer is rectified by rectifiers B,-Bg, the variable 
ing the ac power line voltage. The use of ac sources means component is filtered out by the filter (Chg, Cf, and Cp), 
that filters and voltage regulators must be built. The ac and the resultant de is fed to the load. 


source can be ether single- or three-phase. 

One of the advantages of using rectified voltages from 
three-phase ac sources is that the ripple on the output is 
less; this simplifies and eases the filter requirements of the 
voltage regulator. The voltage regulators considered in 
the article are designed for loads on the order of hundreds 
of watts, and use saturable reactors for regulation. These 
regulators permit ac line voltage fluctuations of + 5-10%, 

a temperature change of the surrounding medium from -30 
to +50°C, and continuous control of the regulated voltage 
over a range of * 20-25%, 

The fluctuations of the regulated voltage are less than 
+0.5-1.5% for this type of regulator, and the full regulation 
time for a step change of load from 100% to 50% and, con- 
versely, from 50% to 100% of the nominal value does not 
exceed 5-15 msec, 

The circuit diagram of the regulator is shown in Fig. 1. 

The regulator consists of the following basic elements: 

1. Saturable reactors Ch, - Chg, made with an external 
feedback winding We and one control winding Wy, shunted 
by the capacitance Cy in order to eliminate the effect of 
the inductive emf of even harmonics. The reactors can 
also be made with an internal feedback winding. 

2. The three-phase transformer Tr (or three single- 
phase transformers). 

3. Power rectifiers B,-Bg. nal 

4. The filter, consisting of choke Ch¢ and capacitances 
Cf, and Cf,. 

5. The meter-amplifier element which, in its turn, 
consists of: a) voltage divider made of resistances R,-Rg, 
and Ry-Ry; b) Zener reference diodes Vc, and Vq,, type 
D809 or D810; c) rectifiers By-By, type D7D-D7Zh; d) 
junction transistors Q,, Q, and Q, (type P4B); e) resistances 
Rs, Rg, and Ry, which limit the currents in the control wind- 
ing and transistors Q,; and Q,. 


In this case, the saturable reactors operate as single- 
phase magnetic amplifiers, working into a variable load 
resistance. A simplified equivalent circuit diagram of the 
load windings of such an amplifier is shown in Fig. 2a. 

Commercial transformers have an insignificant leakage 
reactance, and the quantities x, and x", can be neglected. 
If the load resistance is large, it is also possible to neglect 
the choke resistance (teh), and primary and secondary wind- 
ing resistance (r, and r’,) and the core loss resistance (f%). 
Therefore, the simplified equivalent circuit diagram of the 
load windings of the magnetic amplifier corresponds to the 
circuit shown in Fig. 2b, for no-load on the regulator. 

The amplitude of the output voltage U,,,,, will be de- 
termined, in this case, by the relation between the choke 
(magnetic amplifier) reactance xch, and the transformer 
reactance X). The relation chosen between x,}, and X¢ and 
the given value of the regulated voltage determine the 
transformation ratio of the transformer. 

When the transformer is loaded, Ry <« Xp, and the sim- 
plified equivalent circuit of the load windings of the satur- 
able reactor corresponds to the circuit shown in Fig. 2c. 

The regulation of the output voltage in this type of 
regulator is based on the fact that the reactance of the 
saturable reactor Xo depends on the magnitude of the dc 
excitation ampere-turns. 

The main excitation of the circuit is caused by the 
load current, that is, by the ampere-turns of the feedback 
winding, which has a coefficient somewhat less than unity 
(this coefficient can even be equal to unity). The load 
current is minimum at no-load, and is determined by the 
currents of the voltage dividers Ry+R, and Ry+Ry, and the 
resistance x, is maximum (in the circuit of Fig. 2b). 

At maximum load, Xap is minimum (in the circuit of 
Fig. 2c), and the change in xch due to load current com- 
pensates for the change in the output voltage due to the 

The ac of the three-phase power line, at a 400 cps change in the voltage drop across the circuit resistance 
frequency, is connected through three single-phase saturable R = Ich + Ty + Tp. 
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Fig. 2. Simplified equivalent circuit diagram of the load 
windings of a single-phase saturable reactor at different 


loads. 


The de excitation of the saturable reactors is caused 
not only by the ampere-turns of the feedback winding Wo, 
but also by the ampere-turns of the control winding Wy- 
The latter is used to get better stabilization of the output 
voltage, and to control the amplitude of the regulated vol- 
tage over a range of 4 20-25%. 


The direction and amplitude of the current in the con- 
trol winding depend on the states of transistors Q, and Q, 
in the metering element. 


Voltage reference diodes (Bo, and B,.) are wired in the 
base circuit of each transistor(Q, and Q). Therefore, the 
Circuit in the transistor base will, in practice, flow on'v 


when the voltage between the moving arm of the divider 
and the emitter of the given transistor(U,. or U,.) is greater 
than the sum voltage of the reference diode (U,, or U,) and 
the rectifiers (Uy7-s) or Uys ~y)- The greater the dif- 
ference in the voltages mentioned, the greater the base cur- 
rent of the transistor will be, and the greater the collector 
current. 

When both transistors are cut off, the circuit diagram 
of the control winding will be practically analogous to the 
circuit shown in Fig. 3a. The ampere-turns of the control 
winding in this case assist the excitation action of the feed~- 
back winding, and by the same token, cause additional de- 
crease of the reactance Xp. 
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Fig. 3. Simplified equivalent-circuit dia- 
grams of the control winding of the saturable 
reactor for different states of the transistors 
in the metering element. 


Fig. 4. Dependence of the output voltage of the regulator 
on the current in the control winding. a) for R,, = 30 ohms, 
b) 4.2 ohms, c) 1.2 ohms, d) minimum value of the regu- 
lator output voltage. 





A change in the collector current of transistor Q, is 
equivalent to a change in the resistance R, , shown in Fig. 
3b. A decrease in the resistance that is, an increase 


® 
in the collector current of this transistor, causes an increase 


in the voltage drop across resistance R, and a decrease of 
the current in the control winding. 

When transistor Q, is fully on, the change in the col- 
lector current in transistorQ, is equivalent to a change in 
resistance R,_, shown in Fig. 3e. A decrease in resistance 

will cause an increase of an opposing current in the 
control winding, and therefore these ampere-turns will 
cause a relative increase in the reactance Xqp. 

The use of feedback with a coefficient approaching 
unity ensures an abrupt change in the reactance x,), and 
consequently, in theaoutput. voltage, when the current in 
the control winding is changed. 

The change in the output voltage of the regulator as 
the current in the control winding is changed is shown in 
Fig. 4; the control winding of this reguiator has been dis- 
connected from the metering element, and connected to 
an independent source of dc. 

It can be seen from Fig. 4 that to change the output 
voltage of the given regulator from 12 to 25 v acrossa load 
resistance of 30 ohms, it is sufficient tc change the control 
current from -25 to +25 ma, and from -50 to +50 ma for 
small load resistances. There are 70 turns on the control 
winding of the reactors of the given regulator. It can also 


be seen from Fig. 4 that a negative control current exceed - 
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ing 50 ma is not permitted in this type of regulator. In or- 
der to hold the output regulated voltage at 25 v, for example, 
the control current must be altered within the limits of 22 
to 52 ma, and to hold this voltage at 20 v, the control cur- 
rent must be altered between -6 to +8 ma. 


All the foregoing determines the choice of the resis- 
tances Rs - R; in the metering element. 


When the control winding is connected to the metering 
element, the output voltage is stabilized because a change 
in the output voltage due to a variation in the load, or 
power line voltage, causes the base current of a transistor 
to change; this, in turn, alters the control current so as to 
counteract the output voltage change. 


In principle, it is not possible to compensate completely 
for the voltage deviation by changing the current in the con- 
trol winding, since there must always be some voltage dif- 
ference to cause a new value of the control current. How- 
ever, the output voltage is maintained quite accurately due 
to the large amplification factor. 

Some curves of the dependence of the output voltage 
on the load change are shown in Fig. 5 for one of the regula- 
tors of the type discussed in this paper. These dependencies 
are given for different ac source voltages. 


The amplitude of the regulated output voltage of this 
type of regulator is determined by the position of the moving 
arms of the voltage dividers, which consist of resistances 


Ry, Rg, amd Rg, Ry. 
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The currents in the dividers are significantly greater 


than the base currents of the transistors; therefore the dividers 


are only slightly loaded. 


The voltage obeys the relation 


Cou = OTe = ’ 
for the first divider, where R, is the resistance of the di- 
vider from point A to the moving arm. 

Since the voltage Ue, can be greater under regulation 
than the sum of the voltages Ur, + Uwr-s) by only a very 
small amount, an increase in the resistance Ry (that is, a 
movement of the arm of the divider) will cause a decrease 
in the output voltage amplitude 

The regulated output voltage can be lowered by mov- 
ing the arm of the first divider from point A to point B, 
but only to a value which does not require a change in the 
direction of the current flow in the control winding. Fur- 
ther reduction of the output voltage can be realized by 
moving the arms of both voltage dividers. For finer con- 
trol, regulation of the output voltage must be accomplished 
by moving the arms and the voltage dividers simultaneously, 
so that the resistance ratios of the dividers 





Ry + Rs and/*2 + Ry, 
do not differ much. ty Ry. 

The change in the collector currents of junction tran- 
sistors due to a temperature variation leads to a decrease 
in the control winding current, which in turn causes a 
change in the output voltage of the regulator. The output 
voltage change will cause a corresponding variation in the 
base currents of the transistors, which compensates for the 
collector current change and practically restores the out- 
put voltage amplitude. But when the regulator is heated, 
the change in the collector current of the transistor due to 
a temperature increase must be less than the initial col- 
lector current, so that the decrease in the transistor base 
current can compensate for the increase in the collector 
current. Transistor Q, (Fig. 1) is used to prevent Al, from 
being greater than the initial I-. With no base connection, 
the collector current of Q, is minimum, and its increase 
depends on the temperature of the environment. 

As the temperature increases, the collector current 
change of Q; leads to a decrease in the control winding 
current (Fig. 3d), but an increase in the collector current 
of Qs, conversely, causes an increase in the control winding 
current. Thus, the collector current change in Q, in this 
case compensates for the change in the collector current 
of Q;. 

The output voltage change for a regulator, wired ac- 
cording to the circuit discussed is shown on oscillograms 
(Fig. 6) for various load changes. 


Q! sec 
Saas | 
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10 msec 
1-0 | 
15 nsec 
Un2N p= \, 
Vv i=Wamp i 
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Fig. 6. a) Output voltage change for load current 
changes from zero to 15 amp and from 15 amp 
to zero, b) output voltage change for a load cur- 
rent change from 5 to 10 amp, and from i0 to 5 
amp. 


The recovery time obtained for this particular regula - 
tor is not a limit for this type of regulator. 
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STABILIZING THE TEMPERATURE OF HEATED THERMISTORS 
G. K. Nechaev 
(Kiev) 


Translated from Avtomatika i Telemekhanika, Vol. 21, No. 4, pp. 548-554, April, 1960 
Original article submitted October 12, 1959 


An analytic method of calculating circuits for stabilizing the temperature of heated thermistors 
used in time delay apparatus is presented. 


At present, a number of methods of obtaining time possibility of making such a thermostat is mentioned in the 
delay with electrical circuits are used. Circuits using the literature [2,3]. 
thermal inertia of thermistors are convenient due to their Let us suppose that there is an insulating cylinder on 
simplicity, reliability, and ability to give long time de- whose surface, or inside of which, a heater is placed. By 
lays reaching dozens of minutes. However, these circuits regulating the heater current in a suitable manner, the 
are sensitive to a temperature change of the surrounding temperature inside the cylinder can be held constant when 
medium; this limits their use. the medium temperature changes within certain limits. It 
In order to eliminate this dependence on the external should be noted that the temperature inside the thermostat 
temperature, Udalov [1] proposed that heated thermistors must be higher than the temperature of the surrounding 
be used, whose temperature was regulated in a simple man: medium. 
ner by changing the current of the heater. Two satisfactorily simple circuits for controlling the 
It must be noted that not all heated thermistors are temperature can be proposed; these are shown in Figs. 2 
suitable for this use. Two types of heated thermistors exist; and 3, in which the thermistor acts as the sensing element. 
they are shown in Fig. 1. In the first, the semiconducting In these circuits, Ry is the thermistor in the regulator cir- 
element, a, is inside the heater, b, and hence, heat trans- cuit, r is a linear resistance, rq is the heater resistance, 
fer occurs between the semiconducting element and the Rrc is the semiconducting element of the heated thermis- 
envelope, c, on which the heater is placed. By stabilizing tor, and r, is the barretter. The circuit using the barretter 
the envelope temperature, the heat transfer conditions of has the advantage over the circuit of Fig. 2, in that it is 
the thermistor and its initial temperature are held approxi independent of supply voltage fluctuations, but such a 


mately constant. circuit requires somewhat greater power than the circuit 
The second type of thermistor, in which the heater is With the linear resistance. 
placed inside the semiconducting element, which has a The circuit works in the following manner. At a low 


tubular shape, presents a somewhat different picture. In | temperature of the medium, the thermistor resistance has 
this construction, heat transfer occurs between the semicon a considerable magnitude, the current in it is small, and 
ducting element and the medium; a temperature change of 4 large amount of power is liberated in the heater. As the 
the latter cannot be compensated for by changing the heat- temperature of the medium is raised, the current in Ry in- 
er current with the aid of the proposed circuit. Only the creases, and by the same token, the current in, and power 
temperature of the thermistor can be stabilized. Analysis liberated by the heater decrease. Consequently, it is essen- 
shows that this condition is not sufficient to maintain the _ tial to choose the parameters of the circuit so that as the 
delay transient constant when a circuit with the thermistor temperature of the medium changes, the temperature inside 
is turned on, if the temperature of the surrounding medium the thermostat or heated thermistor will remain within given 


changes. Consequently, the second type of thermistor will limits. 
not do for the circuit being considered. It has been experimentally established that the heat 


transfer between the thermistor or thermostat and the sur- 
rounding medium follows Newton's formula well within the 
operating temperature range (1) 


The operation of the circuit can also be made inde- 
pendent from the temperature of the medium by means of 
thermostation of the thermistor. However, this method is 
applicable only in those cases where the thermostat is P = Fa(Ty—T ) = H(T'1 — To) = K(T — T)), 
simple in construction. A heated thermistor of the first where P is the power dissipated by the thermistor or thermo- 
type can be considered as a microthermostat , inside of stat into the surrounding medium, F is the area of the exter- 
which the temperature can be stabilized within certain nal surface of the thermistor or thermostat, a is the surface 
limits. By analogy, it is possible to make a small-sized heat emission coefficient, H is the overall dissipation factor 
thermostat, and place a direct-heated thermistor in it. The of the thermistor, T, and T, are the temperatures of the 
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Fig. 1. Types of heated thermistors. I) semiconducting 
element inside the heater, II) heater inside the semicon- 
ducting element. 


thermostat surface or the surface of the thermistor heater, 
and that of the surrounding medium, T is the temperature 
inside the thermostat or thermistor heater, and K is the ov- 
er-all heat transfer constant. 

Temperature T will be constant if the power dissipated 
by the resistance rp; is equal to the power emitted by the 
thermistor or thermostat into the surrounding medium: 

Py = forn = K(T —T,). (2) 

If the semiconducting or thermostated element liker- 
ates a certain power, it is considered to be a deduction 
from the computed value of P},; when calculations are made. 
In time relay circuits, the semiconducting element is or- 
dinarily unexcited before the main circuit is turned on, and 
thus it is not necessary to account for this power. 

Consequently, the power P;,; should be a linear function 
of the medium temperature T,. The proposed circuits per- 
mit one to obtain a dependence P}{T,) which is nearly 
linear inthe static mode. In order to have this, it is essen- 
tial that the curve P};(T,) assume values at a number of 
points in the temperature range which correspond to a lin- 
ear dependence of P(T ), that is, so that the difference 
AP = P(T,) — Py(T9) is zero at these points. 

Let us take three such points, which correspond to the 
limits of the temperature range and its center, as shown in 
Fig. 4. The curve P};(T,) must have an inflection point 
in order to intersect Pg(T,) at points 1,2, and n. This 
means that d*P};/dT% must be zero at the given point n. 
Let us consider this condition for the proposed circuit. a 

For the circuit of Fig. 2, the power liberated in 1,,, is 
determined by the expression 








Pr 
Pp, =U*r, A ee (3) 
where U is the power supply voltage of the circuit. 
Then 
d*P (4) 
H 
—— 2U? rr, 


0 





((10,)? 4- R, ROL LR, (re + ry) + ry) — 3 (r + ry) R, (RY? 
[2t.(r 4-7) + rryl* ? 








Fig. 2. Circuit for controlling the temperature of a 
heated thermistor. 


whence, 
rr, _R,(2(R,Y—R, R,] 
(R,) + RR, 


r+ry (5) 





The temperature dependence of the thermistor is de- 
fined by the expression 


B B 
Ry = Rr exp (pz) = 1h: " 


where Rry is the thermistor resistance at ternperature Ty 


and B is a constant defining the properties of the semicon- 
ductor from which the thermistor is made. 
Then 
R. =_— B R 
al ti (1) 


R= B 
1 aq (B+ 27)R,, 


Putting expression (7) into formula (5), we get 
. Rrp B — 2T, (8) 


o- = R, ¥. 
r+r 2 B+T, "pr 
The expression for P., 1s written in the form (9) 
U*Rt 3, UR +, 
i . 








P= j = ” 

= [Ry, (r +?y) +rryl? = [Ry, (r + ,) + rr,]* 
for temperatures of the medium T, and T3, corresponding 
to the limits of the working range. 

‘ : 
The ratio a must be equal to the ratio P;/P,, 
u2 

which is computed or determined in the ordinary manner. 





Consequently, 
R 4. Pry ] 
PE Wily Pa: tment, «| os 
ie ae rr, (10) 








Fig. 3. Circuit for controlling the temperature of a 
heated thermistor with a barretter. 


377 












































$0 

25 to SQ 
ov Ih V; 
015 0 i as 


Fig. 4. Dependence of the powers P, P;,, and 
Py on the temperature of the surrounding 
medium. 


Then, according to expressions (8) and (10), we get 


wou rs (Rr, - re ] = €1 (€2 + 1) lf (11) 
p Ry, (Ry, -- ™ €2 (e, + 1p) 





where 


Ry = Ry exp (7 —_ r-) = Rye. 

The unknown B, which characterizes the thermistor, 
can be determined from(11). But since this equation is 
exponential, and is not analytically solved for B, then the 
curve kf B) is plotted, and B, corresponding to a given kD 
is determined from it. Then the resistances r and ry are 
determined. However, there is only one expression, (8), 
which relates these quantities. The second equation can 
be obtained from the condition of maximum efficiency of 
the circuit at one of the points in the temperature range. 
The circuit efficiency is defined by the expression 

P,,- Rory (12) 


{= 


P ~ (RF rth Fry) + rr) 


where P is the power needed by the whole circuit. 
The quantity r, corresponding to a maximum 7, is 


determined from the equality dn/ dr, = 0: 





r= ——, (13) 





Tp is chosen.as the temperature at which the efficiency 
is maximum. Then the equation for determining r is ob- 
tained from (8) and (13): 

14 
Ce PO — 70. and 

The quantity Rrp is chosen to correspond to existing 
types of thermistors. The value of 1, is defined by the 
expression 


Ry © (15) 
— —>— ° 
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5 °K 






B(x) 


K 

P 
Fig. 5. Dependence B( ky) for finding the com- 
puted value of B. 


Then the power supply voltage of the circuit is deter- 
mined. The power, liberated in Ty; at the temperature Tp, 
is found from the formula 





P, p= | ae : Ry (16 
-P ~~ (7 +9) + rr, » 
or after transformation, 
U*r 
(pn SAS ie Thy (am 
u.p (r +- r ( _ -) : 
Next the supply voltage is determined: 
$— , 
Uaetnatny/ =BP, "7 
H 


The next step in the calculation is to determine the 
maximum power dissipated in the thermistor. The power 
liberated in the thermistor is defined by the expression 





Ur, 2 
la lz (r--+ri)+ =| Ry. (19) 


Ry is found from the condition dPy/ dT = 0; it cor- 
responds to maximum Py: 
») Te ca 
T. Fr, Try 
Then 
 RES.s. iad 
T Max ~ 4r (r +- ry) (21) 
The temperature corresponding to Prmax is found 
from (8) and (20): 
BT p 


If T,4 proves to be outside of the working temperattt] 
range, the maximum P- is taken as that for temperature 
Ts. 

To conclude the computation, the maximum static 
error of the circuit is determined. In order to do this, the 
temperature at which SP = P(T ») — Pz;(T,) has extreme 
values is found. The expression 

(28 
2U°R, Ri rr = — K [R,(r + 7,) + Pry). 


is found from the condition d(AP) /dT, = 0. 





ns 


s deter- 
ure T 


(16) 


(17) 


(18) 
ne the 
> power 
ssion 

(19) 


t cor 


(20) 


(21) 


found 


-—_7"F 





Putting the values of R, and Ry from (6) and (7) into 
(23), we get the equation necessary to find Toy. 


B B 3 
Over Th, | Fr, =~ (x * F-) (r+ ry) + ra (24) 


oe) 


H OM 





Since this equation is not solved analytically for 1 64 
one is obliged to use a graphical method. The maximum 
static error of the circuit is found from the formula 


, , Salle 
at, = 1 = 


—(T—T (28) 


om) * 

Let us examine the computation for a circuit with a 
barretter.It is characterized by constancy of the current 
flowing through the barretter. This current is defined by 
the expression 








Ri(r+r.)4+err 
} anes T H H 26 
Rota? (26) 
The power liberated in ry will be 
, Rr 2 21) 
"7 NCEE | : 


By analogy with the foregoing, an expression similar 
to (8) corresponds to the condition for obtaining an inflec- 
tion point on the curve P}¢T 9). If r = oo in the circuit, we 
get from (8) 

r, = Rr. (28) 
P 

It is possible only in rare cases to make a circuit 
which will fulfill this condition. By connecting a resis- 
tance r in the circuit, one can change the total current 
and make it equal to the nominal current of the barretter. 

The quantity B of the required thermistor is deter- 
mined from the given ratio of the powers at the edges of 
the temperature range. 

This ratio corresponds to (11) for the first circuit. 

If the temperature of the heated thermistor is regulated, 
then r, is the given quantity. Choosing a value for Rrps 
we find r: 

Ty, hy 


et 


(29) 


r 


The power voltage for the circuit is determined from 
the values of the lower limit of.barretration and from the 
nominal current of the barretter: 

ee {= (r+ r,) +7? y 


R, ly 


(30) 





|+o. B. 


In this expression, the thermistor resistance is taken 
at the lower limit of the temperature range. When the 
thermistor temperature is raised, the operating point on 
the volt-ampere characteristic of the barretter will be 
shifted to the middle of the barretration zone. 

The power liberated in the thermistor will be 


(31) 





np rr. 2 
P= ‘lRGaFR | 3 


The amount of thermistor resistance which corre- 
sponds to P-r:44y is determined by an expression similar 
to (20). The maximum power liberated ‘in the thermistor 
will be 

rr 
Sea (32) 
™/* 4(r+r,)° 


Hence, the parameters of all the elements in the cir- 
cuit can be determined. In order to illustrate the method 
presented, the calculation of the circuit shown in Fig. 2 is 
presented as an example. 

Example. Compute the circuit elements of the tem- 
perature regulator of a heated thermistor operating in the 
temperature range from 69, = -30,to 6g =+30°C. The 
temperature of the semiconducting element @ = 50°C. 
The power dissipated by the heater at -30°C is PH, * 10 
watts. 

It is assumed that the best transfer between the therm- 
istor and the medium obeys Newton's formula. 

Let us find the power dissipated by the thermistor at 
temperatures 6, = 30°C and at 6p = 0°C: 


8 — 002 





Pus = Pm p—o,, = 2° Ww, 
6—90 
Foi Pm a 
P 
= 6.25 w and ko =p =4. 
Pus 
Table 1 ; 
B,° K  |2500|3000/3500 





ky = |2,83/3,92| 4,8 

















In order to determine B corresponding to k,, = 4, the 
function k,(B) is computed from (11) in some range of 
change in B. The results of the calculation are presented 
in Table 1. 

Having plotted a curve through these points (Fig. 5), 
we find from it that B = 3025°K, corresponding to k, = 4. 

We will take a thermistor having Rg = 20 ohms at F4=298°K. 

Then the values of the thermistor resistance at three 
points in the temperature range will be 


R, =165ohms R, = 14,2 ohms R, = 42.8 ohms 
1 2 Pp 


Using (8), we find 
y = 0.376. 
The quantity r = 38.5 ohms is found from (14), written 
oe ?_ 30.5 r~ 305 =0. 

The heater resistance is determined from (15): 

Thy = 27.8 ohms. 
The power circuit voltage is computed from (18): 

U = 43.3 volts. 
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Table 2 
































s —30 —15 0 15 30: 
e 8.25 3.97 2.14 :.77 0.714 
P. w 10.0 8.13 6.25 4.38 2.9 
Py Ww 9.9 8.22 6.25 4.24 2.58 
ro w 1.65 2.86 4.05 5.13 5.0 
~ GS 49.2 50.8 50.0 48.8 50.6 





The maximum power dissipated in the thermistor will cuit can be successfully used for regulating the temperature 
be Prmax = 5.18 w, according to (21). of the semiconducting element in a heated thermistor. 

The thermistor resistance which corresponds to P+ 2x 
is found from (20): 


Rr, = 16.1 ohms. 
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USE OF HEATER RESISTANCES FOR CURRENT 


AND VOLTAGE STABILIZATION 
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Translated from Avtomatika i Telemekhanika,Vol. 21, No. 4, pp. 555-557, April, 1960 


Original article submitted October 8, 1959 


A new type of electrical stabilizer with metal heater resistance is considered and analyzed. Results of ex- 
perimental investigation are given,and technical characteristics are noted. 


Bridge stabilizers with nonlinear resistances — incan- 
descent lamps, thermistors, and barretters — are widely 
used as low-power electric current stabilizers [1]. 

Bridge stabilizers have a number of advantages — sim- 
plicity, wide frequency response, and the absence of mov- 
ing parts. On the other hand, such instruments have a com- 
paratively low stabilization coefficient; furthermore, they 
cannot be used in circuits where the load and the supply 
source are commoned. 

In the stabilizer developed by the author the basic 
advantages of bridge stabilizers are retained, and the stabi- 
lization coefficient is increased. The stabilizer is suitable 
for applications in circuits where the load and the supply 
source are commoned. 

The stabilizer consists of a shunt and a metallic resis- 
tance with indirect heating. The heater resistance is re- 
presented by a platinum heater in a glass insulating tube, 
around which a platinum spiral is wrapped (Fig. 1). The 
heater and spiral are placedinan unevacuated brass con- 
tainer serving to protect the heater resistances trom mech- 
anical damage. The diameter of the heater is usually 
within the limits of 10 to 50 y; the inner diameter of the 
insulating tube is approximately double that of the heater. 
The diameter of the spiral wire varies within the range of 
5 to 30 wp. The heater resistance is attached to a socket 
from a type 6N8S radio tube. 

The time constant rT of the heater resistances as a 
function of heater and spiral dimensions is within the lim- 
its 0.03 to 0.80 sec. 

As indicated by the investigation, the variation of the 
spiral resistance R, with concurrent action of the heating 
current I and the current flowing directly through the spiral 
I, is of the form 

Ry = Nyy + 12 + lh (1) 


where Ry, is the initial resistance of the sensing element 
(spiral), c and c, are constants depending upon the con- 
figuration of the heater resistance, a and 6 are constants 
within wide limits independent of the lay of winding and 
of the spiral diameter. 

For platinum heaters a © 2.2, where the temperature 
rise of the heater due to the heating current is within the 





View without cover 





Fig. 1. Metallic resistance 
with indirect heater. 1) Cov- 
er; 2) heaterresistance; 3) 
compensating resistance; 4) 
socket. 











Fig. 2. Circuit diagram of the stabilizer. 
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limits of 40 to 400°. For heater resistances with nichrome 
or constantan heaters, « © 1.7 under the same conditions. 

With a platinum spiral, 8 * 2.0 under the condition 
that the temperature rise of the spiral due to the current 
flowing through it does not exceed 20°. With increasing 
temperature 8 does increase somewhat. 


The arrangement of the stabilizer with heater resistance 


is shown in Fig. 2. 
The relationship between current I, in the load R, and 
the heating current I is of the form 


“por Ro 
s 


B+ 8, 4+-B ° (2) 


Differentiating (2) while taking (1) into account, we 
obtain 


a (Ry ++ Ra + Ry) —1Ry(acl®! +. 0,31 G8 
dl, 2 (ft; -- Jt, 4+- Rs) — 2\ ac “1 Cipa'g =) (3) 


dl (Ay + Re + £s)* 
From that, we find,after rearrangement , 








dl. ( eBlg 
a Ry. 


pice, 58+ Bo — col”) 
dl -Rs+ Rs! 


(Ry -+ Re+ Ra)? : 


Equating the expression obtained to zero, we obtain 
the condition where the current in the load is independent 
of the current flowing through the heater®; 

Ry + Ry + Rs =a I*. (4) 

In this manner the condition for stabilization will ap- 
pear when the resistance increase of the spiral due to the 
heating current, multiplied a times, is equal to the sum 
of the resistances of the spiral, load,and shunt. 

We will determine the conditions for the stabilization 
coefficient of the stabilizer to be a maximum. For that 
purpose we determine the conditions for a minimum of the 
second derivative d*I, /dI* in the operating sector of the 
stabilizer characteristic, that is, in the region where the 
first derivative dl,/dI is equal to zero. 

Differentiating equation (3) and setting the first de- 
rivative equal to zero, we obtain, after a number of re- 
arrangements, 











ne Re (1 — a) 5) 
di? e818 
explh 
i iderably 1 
The expression Rat Ret Pe is considerably less 


than unity. Therefore, the magnitude of the current stabil- 
ized has practically no influence on the stabilization coef- 
ficient. 

It follows from (5) that with decreasing shunt resistance 
R, the second derivative decreases and, consequently, the 
stabilization coefficient increases. In that case, according 
to (4), the operating temperature of the heater resistance 
is simultaneously decreased. The stabilization coefficient 
also increases,due to resistance increase of load Rs, de- 
crease of the nonlinear exponent of heating current a, and 


ting temperature of the heater resistance increases and its 
overload capacity decreases. 

According to (4) a higher load resistance is permissible 
for larger values of the nonlinear exponent a; the power 
output of the stabilizer increases simultaneously. There- 
fore,a platinum wire, which has a higher nonlinear exponent 
than nichrome wire, is used as a heater in the stabilizer. 

It was shown by the experimental investigation of a 
number of stabilizers devised using various metallic heater 
resistances that, with an increase of the heater current, the 
load current initially increases, and, thereafter, after reach- 
ing some maximum value, decreases. This is due to cur- 
rent redistribution between the load and the shunt as a re- 
sult of intense resistance increase of the spiral. 

The stabilization coefficient of stabilizers is,on the 
average equal to 50 fora 10% change in the heating (stabil- 
izing) current, and about 30 for a 20% change of the stabil- 
izing current. 

The resistance of the platinum heater increases signifi- 
cantly upon increase of heater current. Consequently, on 
varying the voltage applied to the heater, the resistance of 
the latter changes in such a fashion as to oppose a change of 
the heater current. As a result the stabilizer will be less 
sensitive to variations of terminal voltage U,p, than to 
heater-current variations. 

The table shows the relationship between the load cur- 
rent I; and the voltage U, transmitted to the terminals of 
a stabilizer with heater resistances of the platinum heater 
being of 504 indiameter and a platinum spiral with wire 
diameter of 24 », with R, equal to 1.0 ohm, Rg equal to 
5.0 ohms, and at an ambient temperature t of 20 and 40°. 

As can be seen from the table, the stabilization coeff: 
cient of a voltage stabilizer is about 70 for a 20% change 
in the voltage. 

The load current decreases when the ambient tempera- 
ture changes as a result of resistance increase of the spiral 
and of the heater. For current stabilization, as well as for 
voltage stabilization, this error can be almost fully com- 
pensated for by winding part of the shunt resistance R, with 
copper. 

For current stabilizers the copper resistance is 30% of 
the total in R, , and for voltage stabilizers is 50%. The 
compensating copper resistance, wound on a small plate, 
is located in the same container as the heater resistance 
(Fig. 1). 

During ambient-temperature changes the region of 
stabilization is somewhat displaced. This is caused by the 
fact that whenever the ambient temperature increases, the 
condition for stabilization (4) is disturbed — resistance R, 
increases, and the heating current within the voltage stab- 
ilizer decreases slightly. The temperature error of stabil- 


izers does not exceed 0.2% per 10°. 
To protect the stabilizers from large overloads, safety 


fuses made of platinum wire of a diameter equal to ap- 
proximately 0.8 of the heater diameter should be connect- 
ed in series with the heater. ©:8 18 





increase of the current stabilized; but in that case the opera- *In satisfying equation (4) the expression Ri + Ret Rs ¢ i. 
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t = 20° 





U. v 0.5 | 1.0 | 1.2 | 1.3 | 1.4 


1.5 | 1.6 | 1.7 | 14.8 | 1.9 | 2.0 





7,,ma 9.20 | 13.20 | 13.85 | 14.02 | 14.17 
































14.23 | 14.25 | 14.25 | 14.24) 14.20) 14.12 
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The stabilizers developed are suitable for stabilization 
of direct current and alternating current (voltage), in the 
frequency region of from a few cycles to several hundred 
kilocycles, without wave-form distortion. 

As compared with thermistors, the stabilizer with 
heater resistance has far lower inertia, higher stability, 
and insignificant temperature error. 

A significant advantage of the stabilizer is that it 
will, after an interruption in operation, reestablish, with 
a high degree of accuracy, previous values of stabilized 
current (voltage). It was shown by an investigation of 


three months’ duration that the voltage change at the 
stabilizer terminals will not exceed + 0.15% of the ini- 
tial value. 

The main advantages of the stabilizer are simplicity, 
the presence of a commoned point between the load and 
the supply source, and a high stabilization coefficient. 
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LETTER TO THE EDITOR 


Translated from Avtomatika i Telemekhanika Vol. 21, No. 4, p. 558, April, 1960 


The work of the author [1] cites new equations for 
basic parameters of the two-point regulation process with 
a new time correction method. The conclusion is drawn 
on the basis of introduction of the concept of residual 
delay At’ so that the method of calculation approaches 
the method laid out in [2]. However, for practical calcu- 
lation of a system with the proposed method of correction, 
it is expedient to express the equations for the regulation 
process parameters directly in terms of given values Qi), 
Qout C, At,and Ton. The time of starting the heat in- 
flux into the object Ton 1s rigorously given by the correc- 
tion device. In that case the equations will be as follows: 


Ton (Qin — Sout — 4#@ out 


AP m+) = —— 


(1) 


Q (2) 
= Ad —et———, 


As m(—) 


m (—) © 


The full amplitude of autooscillation 

, , , , Q; Pin C 
BP im = AP ing) + AP my —) = Ton eas 
2 in— Pout 


out 


(3) 
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T oft =7 on 


= Tos (n — 1), 


= Tog + T org = Toff 
Equations (1) to (5) given can be obtained from 
equations (1) to (5) of [1], if it is considered that 
At’ = T' 5, — At Sout, 
in a Qout 


In addition to the concept of degree of correction k, 
introduced in [1], it is also expedient to introduce a sec- 


ond degree of correction k,, which represents the relation- 

ship of full oscillation amplitudes of the process without 

correction to the amplitudes of the process with correc- 

tion, that is, 

AS 

k= > 
Ad. 

It is apparent that 
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Expressing the full period of autooscillations in the 
system with correction in terms of a new degree of cor- 
rection k, and a power of influx n, we obtain 


Atn? 


7 * hs" 


The curves of the relationship between the autoos- 
cillation period T’ and a power of the influx n for At = 1, 
and at various values of the degree of correction k, will 
be of a form coinciding with the curves in Fig. 3 of [1]. 
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